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REAL PARTY-IN-INTEREST 

The assignee of the entire right, title, and interest in the patent application 
Bechtel BWXT Idaho, LLC. 



RELATED APPEALS AND INTERFERENCES 

There are three related appeals of other United States patent applications that 
may directly affect, or be directly affected by, or have a bearing on, the Board's decision. 
The other applications are all divisional applications of the application that is the subject 
of this appeal. Application, serial no. 10/269,807, filed October 10, 2003, came before the 
Board in Appeal No. 2005-0855. The Board mailed its decision to the appellant on May 
31, 2005. United States Patent No. 7,058,153 subsequently issued on June 6, 2006. 
Application, serial no. 10/383,096, filed March 5, 2003, is currently pending on appeal in 
the briefing stage. In application, serial no. 10/788,743, filed February 25, 2004, a notice 
of appeal was filed on May 24, 2006. There are currently no related interferences known 
to appellant, appellant's legal representative, or the assignee which will directly affect, or 
be directly affected by, or have a bearing on, the Board's decision. 



STATUS OF THE CLAIMS 

Claims 20-24, 26-28, and 30-40 are pending in the application. Claims 20-24, 26- 
28, and 30-40 currently stand rejected. The rejections of claims 20-24, 26-28, and 30-40 
are appealed. Appellant agrees with the examiner that claim 8 has been cancelled. 
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STATUS OF AMENDMENTS 

No amendments were filed or entered subsequent to the final office action mailed 
February 14, 2006. 



r 



SUMMARY OF CLAIMED SUBJECT MATTER 



The present invention is directed to apparatus for performing non-destructive 
testing of materials using positron annihilation. The invention as claimed is summarized 
below with reference to the independent claims and claims separately argued. Claims 20, 
26, 31, 34, 36 and 37 are independent claims. Claims 21, 22, 24, 27, 28, 32, 33, 35 and 38 
are dependent claims argued separately. The claims contain reference numerals and 
reference to the specification and drawings. All references are shown in the application at 
least where indicated herein. 

20. Non-destructive testing apparatus (10, Figure 1, ^fl[ 19, 20, and 34-45), comprising: 
a photon source (12, Figures 1 and 5, Iflf 19, 20, and 34-37), said photon 
source (12) producing photons (16, Figures 1 and 5, fflj 19, 20, and 34-37) having a 
predetermined energy and directing the photons (16) toward a specimen (18, Figures 1 and 
5, Iflf 19-30, 34, 39-46, and 50-60) being tested, the photons (16) from said photon source 
(12) resulting in the creation of positrons within the specimen (18) being tested; 

a detector (14, Figures 1 and 5 19, 20, and 38-40) positioned adjacent the 
specimen (18) being tested, said detector (14) producing raw data (22, Figure 1, fflf 19, 24, 
27-29, 41, 50, and 55) indicative of a positron annihilation event; and 

a data processing system (24, Figure 1, 19, 20, and 41-43) operatively 
associated with said detector (14) and said photon source (12), said data processing system 
(24) operating in accordance with a normal activation/analysis process (38, Figures 2 and 
3, Iff 22-24 and 48-50) when a half-life of a selected positron emitter within the specimen 
(18) being tested is greater than a predetermined half-life, said data processing system (24) 
operating in accordance with a rapid activation/analysis process (40, Figures 2 and 4, fflf 
22, 25, 26, 48, and 51-54) when a half-life of the selected positron emitter within the 
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specimen (18) being tested is less than the predetermined half-life, said data processing 
system (24), when operated in accordance with the rapid activation/analysis process (40), 
alternatively activating said photon source (12) and detecting raw data (22) indicative of a 
positron annihilation event, said data processing system (24) including a Doppler 
broadening algorithm (62, Figure 6, 27-30, 55, 56, 58, and 60), said Doppler 
broadening algorithm (62) processing raw data (22) indicative of a positron annihilation 
event to produce output data (26, Figure 1, Tfl[ 19, 20, and 41-43) indicative of the presence 
or absence of a lattice defect in the specimen (18) being tested. 

21. The non-destructive testing apparatus (10) of claim 20, wherein said detector (14) 
produces raw data (22) indicative of a positron formation event, and wherein said data 
processing system (24) includes a positron lifetime algorithm (64, Figure 6, ]flf 27-30, 55, 
57, and 60), said positron lifetime algorithm (64) processing raw data (22) indicative of a 
positron formation event to produce output data (26) indicative of a changing presence or 
absence of a lattice defect. 

22. The non-destructive testing apparatus (10) of claim 20, further comprising a 
second detector (14) positioned adjacent the specimen being tested, said second detector 
(14) producing raw data (22) indicative of a positron formation event, wherein said data 
processing system (24) includes a positron lifetime algorithm (64, Figure 6, 1fl[ 27-30, 55, 
57, and 60), said positron lifetime algorithm (64) processing data indicative of a positron 
formation event to produce output data (26) indicative of a changing presence or absence 
of a lattice defect. 
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24. The non-destructive testing apparatus (10) of claim 20, wherein said data 
processing system (24) includes a three-dimensional imaging algorithm (66, Figure 6, Tfl[ 
27-30, 55, and 58-60), said three-dimensional imaging algorithm (66) processing raw data 
(22) indicative of a positron annihilation event to produce output data (26) indicative of a 
location of the presence or absence of a lattice defect within the specimen (18) being 
tested. 

26. Non-destructive testing apparatus (10, Figure 1, Tflf 19, 20, and 34-45), comprising: 
positron activation means (12, Figures 1 and 5, 19, 20, and 34-37) for 

activating a positron emitter within a specimen (18, Figures 1 and 5, fflj 19-30, 34, 39-46, 

and 50-60) being tested; 

detector means (14, Figures 1 and 5 Iflf 19, 20, and 38-40) for detecting a 

positron annihilation event within the specimen (18) being tested and for producing raw 

data (22, Figure 1, Iflf 19, 24 > 27-29, 41, 50, and 55) indicative of the positron annihilation 

event; 

means for alternately activating (40, Figures 2 and 4, fflf 22, 25, 26, 48, and 
51-54) the positron emitter within the specimen (18) being tested and detecting a positron 
annihilation event; and 

data processing means (24, Figure 1, fflf 19, 20, and 41-43) operatively 
associated with said detector means (14), said data processing means (24) processing raw 
data (22) indicative of the positron annihilation event in accordance with a Doppler 
broadening algorithm (62, Figure 6, Iflf 27-30, 55, 56, 58, and 60) to produce output data 
(26, Figure 1, fflf 19, 20, and 41-43) indicative of the presence or absence of a lattice defect 
in the specimen (18) being tested. 
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27. The non-destructive testing apparatus (10) of claim 26, wherein said detector 
means (14) detects a positron formation event and a positron annihilation event and 
produces raw data (22) indicative of the positron formation event and the positron 
annihilation event, and wherein said data processing means (24) processes raw data (22) 
indicative of the positron formation event in accordance with a positron lifetime algorithm 
(64) to produce output data (26) indicative of a changing presence or absence of a lattice 
defect. 

28. The non-destructive testing apparatus (10) of claim 26, further comprising second 
detector means (14) for detecting a positron formation event and for producing raw data 
(22) indicative of the positron formation event, wherein said data processing means (24) 
processes raw data (22) indicative of the positron formation event in accordance with a 
positron lifetime algorithm (64) to produce output data (26) indicative of a changing 
presence or absence of a lattice defect. 

31. Non-destructive testing apparatus (10, Figure 1, fflf 19, 20, and 34-45), comprising: 
a photon source (12, Figures 1 and 5, fflf 19, 20, and 34-37), said photon 
source (12) producing photons (16, Figures 1 and 5, Tffl 19, 20, and 34-37) having a 
predetermined energy and directing the photons (16) toward a specimen (18, Figures 1 and 
5, 1fl[ 19-30, 34, 39-46, and 50-60) being tested, the photons (16) from said photon source 
(12) resulting in the creation of positrons within the specimen (18) being tested; 

a detector (14, Figures 1 and 5 fflf 19, 20, and 38-40) positioned adjacent the 
specimen (18) being tested, said detector (14) producing raw data (22, Figure 1, fflj 19, 24, 
27-29, 41, 50, and 55) related to a positron annihilation event; and 

a Doppler broadening processor (62, Figure 6, 27-30, 55, 56, 58, and 60) 
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operatively associated with said detector (14) and responsive to the raw data (22) produced 
thereby, said Doppler broadening processor (62) producing output data (26, Figure 1, Iflf 
19, 20, and 41-43) indicative of the presence or absence of a lattice defect in the specimen 
(18) being tested. 

32. The non-destructive testing apparatus (10) of claim 31, further comprising three- 
dimensional imaging apparatus (66, Figure 6, fflf 27-30, 55, and 58-60) operatively 
associated with said detector (14) and responsive to the raw data (22) produced thereby, 
said three-dimensional imaging apparatus (66) producing output data (26) indicative of a 
location of the presence or absence of a lattice defect within the specimen (18) being 
tested. 

33. The non-destructive testing apparatus (10) of claim 31, wherein said detector (14) 
produces raw data (22) that include data indicative of a positron formation event and data 
indicative of a positron annihilation event, said non-destructive testing apparatus (10) 
further comprising a positron lifetime processor (64, Figure 6, 1fl[ 27-30, 55, 57, and 60) 
operatively associated with said detector (140 and responsive to the raw data (22) 
produced thereby, said positron lifetime processor (64) producing output data (26) 
indicative of the presence or absence of a lattice defect of the specimen (18) being tested 
and indicative of a changing presence or absence of a lattice defect. 

34. Non-destructive testing apparatus (10, Figure 1, fflf 19, 20, and 34-45), comprising: 

a photon source (12, Figures 1 and 5, ^flj 19, 20, and 34-37), said photon 
source (12) producing photons (16, Figures 1 and 5, fflf 19, 20, and 34-37) having a 
predetermined energy and directing the photons (16) toward a specimen (18, Figures 1 and 
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5, flf 19-30, 34, 39-46, and 50-60) being tested, the photons (16) from said photon source 
(12) resulting in the creation of positrons within the specimen (18) being tested; 

a detector (14, Figures 1 and 5 fflf 19, 20, and 38-40) positioned adjacent the 
specimen (18) being tested, said detector (14) producing raw data (22, Figure 1,^19, 24, 
27-29, 41, 50, and 55) indicative of a positron formation event and a positron annihilation 
event; and 

a positron lifetime processor (64, Figure 6, 1fl[ 27-30, 55, 57, and 60) 
operatively associated with said detector (14) and responsive to the raw data (22) produced 
thereby, said positron lifetime processor (64) producing output data (26, Figure 1 , fflf 1 9, 
20, and 41-43) indicative of a the presence or absence of a lattice defect in the specimen 
(18) being tested and indicative of a changing presence or absence of a lattice defect. 

35. The non-destructive testing apparatus (10) of claim 34, further comprising three- 
dimensional imaging apparatus (66, Figure 6, fflf 27-30, 55, and 58-60) operatively 
associated with said detector (14) and responsive to the raw data (22) produced thereby, 
said three-dimensional imaging apparatus (66) producing output data (26) indicative of a 
location of the presence or absence of a lattice defect within the specimen (18) being 
tested. 

36. Non-destructive testing apparatus (10, Figure 1, fflj 19, 20, and 34-45), comprising: 

a photon source (12, Figures 1 and 5, fflf 19, 20, and 34-37), said photon 
source (12) producing photons (16, Figures 1 and 5, fflj 19, 20, and 34-37) having a 
predetermined energy and directing the photons (16) toward a specimen (18, Figures 1 and 
5, 1fl[ 19-30, 34, 39-46, and 50-60) being tested, the photons (16) from said photon source 
(12) resulting in the creation of positrons within the specimen (18) being tested; 



a detector (14, Figures 1 and 5 19, 20, and 38-40) positioned adjacent the 
specimen (18) being tested, said detector (14) producing raw data (22, Figure 1, fflf 19, 24, 
27-29, 41, 50, and 55) indicative of a positron formation event and a positron annihilation 
event; and 

a data processing system (24, Figure 1, fflf 19, 20, and 41-43) operatively 
associated with said detector (14), said data processing system (24) including: 

a Doppler broadening algorithm (62, Figure 6, fflf 27-30, 55, 56, 58, 
and 60), said Doppler broadening algorithm (62) processing raw data (22) 
indicative of a positron annihilation event to produce output data (26, Figure 
1, IHf 19, 20, and 41-43) indicative of a presence or absence of a lattice defect 
in the specimen (18) being tested; 

a positron lifetime algorithm (64, Figure 6, Iffl 27-30, 55, 57, and 
60), said positron lifetime algorithm (64) processing raw data (22) indicative 
of a positron formation event to produce output data (26, Figure 1, fflf 19, 20, 
and 41-43) indicative of a changing presence or absence of a lattice defect; 
and 

a three-dimensional imaging algorithm (66, Figure 6, fflj 27-30, 55, 
and 58-60), said three-dimensional imaging algorithm (66) processing raw 
data (22) indicative of a positron annihilation event to produce output data 
(26, Figure 1, Iffl 19, 20, and 41-43) indicative of a location of the presence or 
absence of a lattice defect within the specimen (18) being tested. 

37. Non-destructive testing apparatus (10, Figure 1, 1fl[ 19, 20, and 34-45), comprising: 
a photon source (12, Figures 1 and 5, 1fl| 19, 20, and 34-37), said photon 
source (12) producing photons (16, Figures 1 and 5, Tflf 19, 20, and 34-37) having a 
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predetermined energy and directing the photons (16) toward a specimen (18, 
Figures 1 and 5, ^ 19-30, 34, 39-46, and 50-60) being tested, the photons (16) 
from said photon source (12) resulting in the creation of positrons within the 
specimen (18) being tested; 

a detector (14, Figures 1 and 5 19, 20, and 38-40) positioned adjacent 
the specimen (18) being tested, said detector (14) producing raw data (22, Figure 1, 
fflf 19, 24, 27-29, 41, 50, and 55) indicative of a positron annihilation event; and 

a data processing system (24, Figure 1, 19, 20, and 41-43) operatively 
associated with said detector (14) and said photon source (12), said data processing 
system (24) including: 

a Doppler broadening algorithm (62, Figure 6, fflf 27-30, 55, 56, 58, and 
60), said Doppler broadening algorithm (62) processing raw data (22) 
indicative of a positron annihilation event to produce output data (26, 
Figure 1, fflf 19, 20, and 41-43) indicative of a presence or absence of a 
lattice defect in the specimen (18) being tested; and 

a positron lifetime algorithm (64, Figure 6, fflf 27-30, 55, 57, and 60), said 
positron lifetime algorithm (64) processing raw data (22) indicative of a 
positron formation event to produce output data (26) indicative of a 
changing presence or absence of a lattice defect, 
said data processing system (24) operating in accordance with a normal 
activation/analysis process (38, Figures 2 and 3, ^fl|22-24 and 48-50) when a half- 
life of a selected positron emitter within the specimen (18) being tested is greater 
than a predetermined half-life, said data processing system (24) operating in 
accordance with a rapid activation/analysis process (40, Figures 2 and 4, 1fif 22 > 25, 
26, 48, and 51-54) when a half-life of the selected positron emitter within the 
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specimen (18) being tested is less than the predetermined half-life, said data 
processing system (24), when operated in accordance with the rapid 
activation/analysis process (40), alternatively activating said photon source (12) 
and detecting raw data (22) indicative of a positron annihilation event, said data 
processing system (24) using said Doppler broadening algorithm (62) and said 
positron lifetime algorithm (64) to process raw data (22) indicative of a positron 
annihilation event and produce output data (26) indicative of the presence or 
absence of a lattice defect in the specimen (18) being tested. 

38. The non-destructive testing apparatus (10) of claim 37, wherein said data 
processing system (24) further comprises a three-dimensional imaging algorithm (66, 
Figure 6, fflf 27-30, 55, and 58-60) and wherein said data processing system (24) uses said 
three-dimensional imaging algorithm (66) to produce output data (26) indicative of a 
location of the presence or absence of a lattice defect within the specimen (18) being 
tested. 
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GROUNDS OF REJECTION TO BE REVIEWED ON APPEAL 

1 . Whether claims 20-24, 26-28 and 30-40 are patentable under 35 U.S.C. § 1 1 2, 
first paragraph, as complying with the enablement requirement. 

2. Whether claims 20-24, 26-28 and 30-40 are patentable under 35 U.S.C. § 1 1 2, 
second paragraph, as distinctly claiming the subject matter of the invention. 



- 15 - 



ARGUMENT 



Opening Statement 

This appeal truly represents a situation of re-plowing ground already plowed. The 

precise issues presented here regarding the enablement of the claims and sufficiency of the 

present specification have already been decided by the Board in the appeal described 

above. Whereas issued U.S. Patent No. 7,058,153 (the subject of the previous appeal) is 

directed to method claims, comprising use of one or more of the various algorithms 

described in the present specification, the present application is directed to apparatus 

claims comprising, for example, a data processing system including one or more of the 

various algorithms described in the present specification. In the previous appeal (like this 

one), the examiner had rejected the claims for indefiniteness, and lack of enablement, 

arguing that one of ordinary skill in the art would not know how to choose from among 

the many algorithms referenced in the specification and known in the prior art. The Board 

reversed. Indeed, the Board's treatment of the examiner's rejections in the previous 

appeal applies equally well here. The Board also reversed the examiner's rejections for 

lack of enablement as follows: 

[W]hile the appellant's disclosure does not convey much detail 
as to the various algorithms described therein, these algorithms 
admittedly are generally known in the art. The examiner has 
failed to advance any cogent reasoning as to why the disclosure 
would not have enabled a person of ordinary skill in the art to 
employ these algorithms without undue experimentation. 

Ex parte Akers, Appeal No. 2005-0855, at 8 (Bd. Pat. App. & Interf. 2005) (slip op.). The 

Board reversed the examiner's rejections based on indefiniteness of the claims as follows: 

The examiner views the appealed claims to be indefinite for 
reasons essentially similar to those listed above in connection 
with the enablement rejection. In this regard, the examiner 
points to the claim limitations relation to . . . the positron 
lifetime, Doppler broadening, three-dimensional imaging and 

activation/analysis process algorithms [T]he examiner's 

criticism of the claim limitations pertaining to the algorithms . . 
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. focus[es]on the breadth of these limitations. It is well settled, 
however, that mere breadth does not equate to indefiniteness. 
In re Miller , 441 F.2d 689, 169 USPQ 597, 600 (CCPA 1971). 

Id. at 9. For the same reasons set forth above, in the instant appeal, the examiner failed to 

carry his prima facie burden to demonstrate indefiniteness and lack of enablement. And, 

so, ground already plowed must nonetheless be plowed again. 

I. Whether claims 20-24, 26-28 and 30-40 are patentable under 35 U.S.C. §112, 
first paragraph,* as complying with the enablement requirement. 

A. Legal Standard For Rejecting Claims Under 35 U.S.C. §112, First 
Paragraph 

The legal standard for determining whether the disclosure provides a sufficient 
description of the invention is whether a person reasonably skilled in the art could make or 
use the invention without undue experimentation based on the disclosure and on 
information known in the art. United States v. Telectronics, Inc., 857 F.2d 778, 8 USPQ2d 
1217 (Fed. Cir. 1988). The fact that experimentation may be complex does not necessarily 
make it undue if the art typically engages in such experimentation. In re Wands, 858 F.2d 
731, 8 USPQ2d 1400 (Fed. Cir. 1988). See Ex parte Lemelson, 2002 WL 32334419, *3 
(Bd. Pat. App. & Interf. 2002) ("The fact that an invention is sophisticated and complex 
does not, by itself, lead to the conclusion that undue experimentation would be required to 
make and use the invention.") That is, the test of enablement is not whether any 
experimentation is required, but whether, if experimentation is necessary, it is undue. In 
reAngstadt, 537 F.2d 498, 190 USPQ 214 (CCPA 1976). 

The factors to be considered when determining whether there is sufficient evidence 
to support a determination that a disclosure does not satisfy the enablement requirement 
and whether any necessary experimentation is "undue" include, but are not limited to: The 
breadth of the claims; the nature of the invention; the state of the prior art; the level of one 
of ordinary skill; the level of predictability in the art; the amount of direction provided by 
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the inventor; the existence of working examples; and the quantity of experimentation 
needed to make or use the invention based on the content of the disclosure. In re Wands, 
858 F.2d 731, 737 (Fed. Cir. 1988); MPEP 2164.01(a). It is improper to conclude that a 
disclosure is not enabling based on an analysis of only one of the above factors while 
ignoring one or more of the others. See, e.g., Ex parte Noguchi, 2002 WL 1801470 (Bd. 
Pat. App. & Interf. 2002) (supplying information directed to one factor is insufficient); Ex 
parte Kopetzki, 2001 WL 11197754 (Bd. Pat. App. & Interf. 2001) (merely reciting facts 
without reference to the factors is insufficient); MPEP 2164.01(a). Moreover, whether the 
specification is enabling must be determined with reference only to the limitations 
claimed. CFMT, Inc. v. Yieldup InVl Corp., 349 F.3d 1333, 1338 (Fed. Cir. 2003). 

With regard to the burden of proof required to support a rejection under §112, the 
Patent Office is required to assume that the specification complies with the enablement 
provision of §112 unless it has acceptable evidence or reasoning to suggest otherwise. 
See, e.g., In re Marzocchi, 439 F.2d 220, 169 USPQ 367 (CCPA 1979). The Patent Office 
thus must provide reasons, supported by the record as a whole, why the specification is not 
enabling. Then and only then does the burden shift to the applicant to show that one of 
ordinary skill in the art could have practiced the claimed invention without undue 
experimentation. Gould v. Missinghoff, 229 USPQ 1 (D.D.C. 1985), off 'd in part, vacated 
in part t and remanded sub. nom., Gould v. Quigg, 822 F.2d 1074, 3 USPQ2d 1302 (Fed. 
Cir. 1 987). Mere conclusionary statements as to the level of ordinary skill in the art are 
not a sufficient basis for a rejection under 35 U.S.C. §112. In re Brebner, 455 F.2d 1402, 
173 USPQ 169 (CCPA 1972). In addition, the law does not require, and indeed prefers, 
that a patent specification omit that which is well-known. In re Buchner, 929 F.2d 660, 18 
USPQ2d 1331 (Fed. Cir. 1991). 

B. Response to Examiner's Rejections 
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The examiner failed to carry his prima facie burden of demonstrating that claims 
20-24, 26-28 and 30-40 are not enabled as required by §112, first paragraph. First, the 
examiner failed to analyze the enablement issue with reference to the eight required 
factors. Second, the examiner's rejections are based on limitations that are not part of the 
claims. Finally, given that the Doppler broadening, positron lifetime and three- 
dimensional imaging algorithms are well-known in the art, the level of detail contained in 
the specification sufficiently enables the claims. 

1 . The eight required factors were not addressed. 

First, while the examiner sets forth a laundry list of variables that he argues one of 
ordinary skill in the art must consider when choosing an appropriate Doppler-broadening 
algorithm, positron lifetime algorithm, or three-dimensional imaging algorithm, the 
examiner has nonetheless failed to meet his prima facie burden. In rejecting claims for 
non-enablement, the examiner was required to analyze each of the eight required factors 
set forth in In re Wands, 858 F.2d 731, 737 (Fed. Cir. 1988), with reference to the claims, 
which he did not. Nor did the examiner supply facts to support each and every factor 
required by In re Wands. The Board has repeatedly reversed similar rejections, holding 
that the examiner failed to meet his prima facie burden. See, e.g., Ex parte Noguchi, 2002 
WL 1801470 (Bd. Pat. App. & Interf. 2002) (supplying information directed to one factor 
is insufficient); Ex parte Kopetzki, 2001 WL 11197754 (Bd. Pat. App. & Interf. 2001) 
(merely reciting facts without reference to the factors is insufficient). Even assuming that 
the examiner put forth any "facts" (which is denied), the "mere citing of a series of facts, 
without relating them to the criteria for determining undue experimentation set forth in In 
re Wands, is not sufficient to establish a prima facie case of nonenablement." Ex parte 
Kopetzki, 2001 WL 11197754, at *3. In addition, by reviewing the variables that one of 
ordinary skill in the art would know to consider, the examiner makes appellant's point that 
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one of ordinary skill in the art would well know how to pick a suitable algorithm, select 
appropriate constants and combine algorithms, making any experimentation routine and 
not undue. 

2. Rejections are improperly based on limitations not in the claims. 

Second, the examiner argues that undue experimentation is necessary to select 
algorithms that will (1) obtain accurate detection of lattice defects (i.e., not too few and 
not too many); (2) obtain reasonable statistical analyses; (3) determine which algorithm 
"actually provides a reasonable determination of the presence of enough lattice defects to 
be of concern;" and (4) determine "whether there are present enough lattice defects to 

cause metal fatigue, etc " Final Office Action, p. 6. However, these limitations are 

not in the claims and, therefore cannot be considered when determining whether the 
specification enables the claims. 

Claims 20, 31, 36 and 37 claim apparatus comprising a data processing system, 
including a Doppler broadening algorithm; or data processing means that process raw data 
in accordance with a Doppler broadening algorithm; or a Doppler broadening processor 
further limited by "producing output data indicative of the presence or absence of a lattice 
defect in the specimen being tested." Claims 21, 22, 27 and 28 claim apparatus 
comprising a data processing system including a positron lifetime algorithm, or data 
processing means that processes raw data in accordance with a positron lifetime algorithm, 
further limited by "producing output data indicative of a changing presence or absence of 
a lattice defect in the specimen being tested." Claims 33 and 34 claim an apparatus 
comprising "a positron lifetime processor producing output data indicative of the presence 
or absence of a lattice defect of the specimen being tested and indicative of a changing 
presence or absence of a lattice defect." Claims 24, 32, 35 and 38 claim an apparatus 
comprising a data processing system including a three-dimensional imaging algorithm, or 
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three-dimensional imaging apparatus further limited by producing "output data indicative 
of a location of the presence or absence of a lattice defect within the specimen being 
tested." None of these claims require accuracy; reasonable statistical analysis; 
"providing] a reasonable determination of the presence of enough lattice defects to be of 
concern;" or determining "whether there are present enough lattice defects to cause metal 
fatigue." As the Federal Circuit has held, "Title 35 does not require that a patent 
disclosure enable ... a perfected, commercially viable embodiment absent a claim 
limitation to that effect." CFMT, Inc. v. Yieldup Int'l Corp., 349 F.3d 1333, 1338 (Fed. 
Cir. 2003). Therefore, the examiner cannot evaluate enablement based on limitations not 
contained in the claims. Id., 349 F.3d at 1338 (since claims contained no standard for 
cleanliness, claims for cleaning methods were enabled for any level of contaminant 
removal); Ex parte Lemelson, 2002 WL 32334419, *4. 

Similarly, but without reference to any claims, the examiner objects to the 
sufficiency of the disclosure in Figure 6, regarding a feedback arrangement between the 
data processing system 60 and the algorithms 62, 64 and 66. Final Office Action, p. 9. 
However, the examiner has not tied these arguments to the elements of any claim to show 
how such elements are not enabled. CFMT, Inc., 349 F.3d at 1338 (level of disclosure 
necessary for enablement varies with the scope of the claims). Therefore, the examiner's 
arguments in this regard cannot be used to demonstrate a prima facie case of non- 
enablement. 

The case of Ex parte Lemelson, 2002 WL 32334419, presents similar facts to the 
ones here and is persuasive in its analysis of enablement. In that case, the invention was 
directed to a method and apparatus for controlling the travel of a powered vehicle, 
comprising calculating whether the vehicle and an object are on a collision course by 
measuring the distance to and relative velocity between the two, and using fuzzy logic 
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interference rules to determine a combination of steering and acceleration to avoid a 

collision. The examiner rejected the claims for lack of enablement on grounds almost 

identical to the ones here: 

The rejection also states that the disclosure is insufficient for 
teaching how to distinguish one object from another object or 
how to determine distance and relative velocity of each object. . 
. . The examiner observes that the object detection and 
recognition and the fuzzy logic rules require complex 
algorithms and extensive processing that are not sufficiently 
described in the disclosure [citation omitted]. In the final 
rejection the examiner added that the complexity of putting all 
the features into one real-time system required undue 
experimentation. 

Id., at *2. The appellants argued that fuzzy logic and image analysis were techniques 
well-known in the art. Board reversed the examiner's rejections. First, the Board held that 
the examiner improperly based his rejections on a real-time operation requirement not 
included in the claims. It also held, "Second, the examiner's findings of undue 
experimentation are mere conclusions based on the examiner's own speculations. The 
examiner essentially finds that there would be undue experimentation because the claimed 
invention requires sophisticated and complex operations. The fact that an invention is 
sophisticated and complex does not, by itself, lead to the conclusion that undue 
experimentation would be required to make and use the invention." Id. 9 at *3. 

In addition, claims 21-24, 39 and 40 depend from claim 20 and are also allowable 
at least for the reasons that claim 20 is allowable. Claims 27, 28, and 30 depend from 
claim 26 and are also allowable at least for the reasons that claim 26 is allowable. Claims 
32 and 33 depend from claim 31 and are also allowable at least for the reasons that claim 
31 is allowable. Claim 35 depends from claim 34 and is also allowable at least for the 
reasons that claim 34 is allowable. Claim 38 depends from claim 37 and is also allowable 
at least for the reasons that claim 37 is allowable. 

3. The specification is enabling. 
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Finally, given that Doppler broadening algorithms, positron lifetime algorithms 
and three-dimensional imaging algorithms are well-known, the level of detail provided in 
the present specification is sufficiently enabling, as the Board determined with reference to 
the appeal in the divisional case issued as U.S. Patent No. 7,058,153. See supra at pp. 16- 
17. The appellant specifically states in the specification that the Doppler broadening 
algorithms, positron lifetime algorithms and three-dimensional imaging algorithms are 
well-known in the art and could be readily implemented by persons having ordinary skill 
in the art after having become familiar with the teachings of the present invention. 
Specification, Iflf [0056] [0059], [0057]. Because these algorithms are well-known to 
those of ordinary skill in the art who are familiar with their applications, the specification 
need not and does not describe them further. At no time has the examiner ever disputed 
the fact that these algorithms are well-known. Since that fact is undisputed, the examiner 
was bound to establish that the amount of experimentation in using these algorithms 
familiar to those of ordinary skill in the art would be undue. As explained above, he failed 
to make such a showing because he did not analyze the undue experimentation factors 
from In re Wands. 

During prosecution, the appellant provided examples of issued patents to show the 
general level of detail required to enable similar inventions for those of ordinary skill in 
the art. See Amendment (12/14/05), pp. 16-19, 21 and 24. In his Summary of Argument 
presented in the appellant's Amendment dated December 14, 2005, the appellant stated 
that "the best measure of level of disclosure required in this particular field of endeavor 
can be obtained by reviewing issued U.S. patents in the same field." Amendment 
(12/14/05), p. 16, Ins. 15-17. The examiner argued that this statement was "misleading," 
because he believed the appellant disclosed only one U.S. patent - U.S. Patent No. 
6,178,218 to Akers (Doppler broadening). Id. at p. 18, In. 23. In fact, including the Akers 
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Patent, the appellant provided examples of three such patents showing enabling 
disclosures for inventions involving the algorithms at issue here, the two others being U.S. 
Patent No. 4,064,438 to Alex, et al. (positron lifetime algorithm; Doppler broadening), 
U.S. Patent No. 4,756,866 to Alvarez (3-D imaging algorithm). Id. at p. 21, In. 6; p. 24, 
In. 5. Further, the other references previously cited by the examiner (i.e., Derlet, 
Bandzuch, Zhu, and Shaffer) illustrate that persons having ordinary skill in the art are 
well-aware of the various issues and complexities involved in selecting, modifying and 
combining any of these algorithms, rendering any experimentation required routine and 
not undue. 

The examiner argues that "the specification indicates that there are different known 
Doppler broadening algorithms but , fails to disclose which of these known algorithms 
would be suitable for use in the present invention." Final Office Action, p. 3. That is not 
the case. The specification specifically discloses a suitable algorithm when it states that 
"in one preferred embodiment of the invention, the Doppler broadening algorithm 62 may 
comprise the Doppler broadening algorithm described in U.S. Patent No. 6,178,218 Bl . . . 
." Page 28, If [0056]. Contrary to the examiner's suggestion, there is no experimentation 
involved at all in selecting the precise algorithm the appellant identified in the 
specification. 

The examiner argues that the specification states that the "positron lifetime 
algorithm and the 3-D algorithm, generally each requires at least two detectors ." Final 
Office Action, p. 8. He next asserts that claims 21, 33, 34, 36 and 37 require the use of 
only a single detector. Thus, the examiner argues that "Appellant's disclosure is clearly 
insufficient and non-enabling as to which specific known positron lifetime algorithm can 
be used with only a single detector and, which specific Doppler broadening algorithm can 
be used with only a single detector." Id. at p. 8. The examiner's argument is incorrect for 
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two reasons: the specification is not limited to the use of two detectors and the claims are 
not limited to the use of one detector. 

Claims 20 (from which claim 21 depends), claim 31 (from which claim 33 
depends), claim 34, claim 36 and claim 37 claim "a non-destructive testing apparatus, 
comprising . . . a detector . . . (Emphasis added.) The Federal Circuit could not have 
been more clear when it cautioned, "This court has repeatedly emphasized that an 
indefinite article 4 a' or 'an' in patent parlance carries the meaning of 'one or more" in 
open-ended claims containing the transitional phrase, 'comprising.'" KCJ Corp. v. Kinetic 
Concepts, Inc., 223 F.3d 1351, 1356 (Fed. Cir. 2000). The Federal Circuit then stated, 
"Moreover, standing alone, a disclosure of a preferred or exemplary embodiment 
encompassing a singular element does not disclaim a plural embodiment." Id. at 1356. 
Thus, in accordance with the holding of the Federal Circuit, use of the article "a" in 
conjunction with the word "detector" and the transitional phrase "comprising" means that 
the claims set forth above comprise "one or more" detectors. The specification enables the 
use of "one or more" detectors. While the specification notes that two detectors may be 
preferred, this does not operate to exclude one detector either. Moreover, even if the 
examiner were correct that the specification enables only the use of two detectors (which 
is denied), the scope of claims 21, 33, 34, 36 and 37 covers two (i.e., one or more) 
detectors. Since one embodiment covered by the claim is enabled, then the enablement 
requirement is satisfied. The words of the Board in another case cannot be improved 
upon: "The examiner's reasoning is logical but not entirely consistent with the law: 
enabling the 'full scope' of a claim does not necessarily require enabling every 
embodiment within the claim." Ex parte Saito, 2005 WL 3524580, *3 (Bd. Pat. App. & 
Interf. 2005). 

The examiner also argues that the specification does not enable claims 21, 22, 27, 
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28, 33, 34, 36 and 37 with respect to the "'changing' characteristics of the system." Final 

Office Action, p. 8. These claims all include a limitation of processing data "to produce 

output data indicative of a changing presence or absence of a lattice defect." Contrary to 

the examiner's assertion, the specification is enabling because it explains: 

In addition, information obtained from the mean lifetime of 
various defect components may be used to derive information 
relating to changing characteristics of the defects present in the 
specimen. The positron lifetime algorithm 64 basically 
involves a determination of the time between positron 
formation and annihilation. In order to do so, the positron 
lifetime algorithm detects some precursor event associated with 
the formation of the positron, as well as the gamma rays 20 
produced by the positron annihilation event. The time between 
these two events is the positron lifetime. 

Specification, U[0057] (emphasis added). The specification goes on to say that because 
positron lifetime algorithms are well-known, they will not be described in further detail. 
Compare Specification, 1f[0057] with U.S. Patent No. 4,064,438, Col. 6, Ins. 60-68 
(detection of hydrogen embrittlement by positron lifetime measurements is accomplished 
according to known methods). Thus, the level of detail set forth in the specification is 
sufficiently enabling for one of ordinary skill in the art. 

Claims 21-24, 39 and 40 depend from claim 20 and are also allowable at least for 
the reasons that claim 20 is allowable. Claims 27, 28, and 30 depend from claim 26 and 
are also allowable at least for the reasons that claim 26 is allowable. Claims 32 and 33 
depend from claim 31 and are also allowable at least for the reasons that claim 31 is 
allowable. Claim 35 depends from claim 34 and is also allowable at least for the reasons 
that claim 34 is allowable. Claim 38 depends from claim 37 and is also allowable at least 
for the reasons that claim 37 is allowable. 

Having failed to analyze the In re Wands factors, the examiner failed to meet his 
prima facie burden to demonstrate that the claims are not enabled. In any event, as the 
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Board previously determined in the divisional case, the descriptions of the various 

algorithms used in detecting positron annihilation events are sufficiently enabling. 

II. Whether claims 20-24, 26-28 and 30-40 are patentable under 35 U.S.C. §112, 
second paragraph, as distinctly claiming the subject matter of the invention. 

A. Legal Standard For Rejecting Claims Under 35 U.S.C. §112, Second 
Paragraph 

Section 112, second paragraph, requires that the claims particularly point out and 
distinctly claim what the applicant regards as his invention. 35 U.S.C. § 112. The 
standard is one of reasonably particularity and completeness. MPEP 2173.02. "The 
indefiniteness inquiry focuses on whether those skilled in the art would understand the 
scope of the claim when the claim is read in light of the rest of the specification." Union 
Pacific Res. Co. v. Chesapeake Energy Corp., 236 F.3d 684, 691 (Fed. Cir. 2001); Credle 
v. Bond, 25 F.3d 1566, 1576 (Fed. Cir. 1994) (all that is required under § 112, 1f2 is that 
one of ordinary skill in the art would understand the claim language based on the 
specification and the drawings). The degree of precision necessary for adequate claims 
depends on the nature of the subject matter. Miles Laboratories, Inc., v. Shandon, Inc., 27 
USPQ2d 1123 (Fed. Cir. 1993). As such, the Patent Office does not permit per se 
rejections of certain language as being indefinite under § 112, ^| 2. MPEP 2173.02 
("Office policy is not to employ per se rules to make technical rejections.") The examiner 
bears the initial burden of presenting evidence or reasoning to explain why persons skilled 
in the art would understand the claim language based on the specification. See In re 
Oetiker, 977 F.2d 1443, 1445 (Fed. Cir. 1992). "After evidence or argument is submitted 
by the applicant in response, patentability is determined on the totality of the record with 
due consideration to persuasiveness of argument." In re Oetiker, 977 F.2d at 1445. 

B. Response to Examiner's Rejections 

The examiner never established a prima facie case of indefiniteness, as explained 
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more fully below. Moreover, even assuming that he did carry his initial burden (which is 

denied), the appellant set forth significant evidence and argument to demonstrate why one 

of ordinary skill in the art would understand the claim language based on the specification. 

Instead of addressing the appellant's argument and evidence regarding definiteness, the 

examiner dismissed them summarily as "not convincing" and "not persuasive" and merely 

incorporated by reference his previous action of September 14, 2005. Thus, the examiner 

failed to carry his ultimate burden of persuasion. 

Independent claim 20 comprises "a data processing system including a Doppler 

broadening algorithm"; independent claim 26, "data processing means processing raw data 

. . . in accordance with a Doppler broadening algorithm"; and, independent claim 31, "a 

Doppler broadening processor". Independent claim 34 comprises "a positron lifetime 

processor"; independent claim 36, "a Doppler broadening algorithm ... a positron lifetime 

algorithm . . . and a three-dimensional imaging algorithm"; and independent claim 37, "a 

Doppler broadening algorithm . . . and a positron lifetime algorithm". The examiner 

previously rejected claims 20-24, 26-28 and 30-40 under 35. U.S.C. §112, second 

paragraph, as failing to particularly point out and distinctly claim the subject matter of the 

invention, because the "claims are vague, indefinite and incomplete as to how and in what 

manner one should select a positron lifetime, Doppler broadening or three-dimensional 

algorithm for a plurality of available algorithms .... Office Action (9/14/05) at p. 8. The 

appellant responded with the following evidence and argument, explaining why the 

examiner had failed to meet his prima facie burden of demonstrating indefiniteness: 

In responding to these rejections, applicant hereby 
repeats the arguments set forth above in response to the 
Section 112, first paragraph, rejections. That is, the various 
algorithms are not only enabled when considered in light of 
the level of ordinary skill in the art, but are also not vague nor 
indefinite for the same reason, i.e., the person having ordinary 
skill in the art would understand the algorithms and possess 
the knowledge required to evaluate the many algorithms 
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available and select those that would be appropriate for the 
desired application. 

In addition, with regard to the Doppler broadening 
algorithm, paragraph 0056 of the specification states that 
several different types of Doppler broadening techniques have 
been developed and could be used. The specification also 
states that the Doppler broadening algorithm may comprise the 
algorithm disclosed in U.S. Patent No. 6,178,218. 

The positron lifetime algorithm is described in 
paragraph 0057 of the specification. That paragraph also 
makes clear that "systems for detecting positron lifetimes, as 
well as the algorithms utilized thereby, are well-known in the 
art and could be easily provided by persons having ordinary 
skill in the art after having become familiar with the teachings 
of the present invention." 

With respect to three-dimensional imaging algorithms, 
paragraph 0059 of the specification also states that several 
different types of 3-D imaging algorithms are known in the art 
and could be used in conjunction with the present invention. 

Clearly, the terms "Doppler broadening algorithm," 
"positron lifetime algorithm," and "three dimensional imaging 
algorithm" are specifically referred to in the specification, thus 
provide sufficient support for their use in the claims. In 
addition, techniques for performing both of these processes are 
well-known to persons having ordinary skill in the art, and 
could be easily provided by persons having ordinary skill in 
the art after having become familiar with the teachings of the 
present invention. For example, the publication cited by the 
applicant in the Information Disclosure statement, namely the 
paper entitled "Positron Annihilation Spectroscopy," discusses 
the Doppler broadening and positron lifetime techniques, as 
well as methods for performing these techniques. Three- 
dimensional imaging techniques are also well-known in the art 
and are commonly used in PET medical scanners. 

Therefore, because Doppler broadening, positron 
lifetime, and 3-D imaging techniques are well known to 
persons having ordinary skill in the art and could be easily 
provided by persons having ordinary skill in the art after 
having become familiar with the teachings of the present 
invention, the references in the specification and claims to 
"Doppler broadening algorithm," "positron lifetime 
algorithm," and "three dimensional imaging algorithm" meet 
the requirements set forth in Miles Laboratories , supra , that is, 
these terms have sufficient meaning to persons having 
ordinary skill in the art to allow them to understand the bounds 
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of the claims when read in light of the specification. 



With regard to the selective activation algorithm, 
applicant notes that the selective activation algorithm is 
described in paragraph 0060 and basically involves controlling 
or operating the photon source as necessary to produce 
photons having energy levels suitable for activating the 
selected positron emitter or emitters. As described in 
paragraph 0036, electron accelerators capable of producing 
photons having a range of energies are well-known in the art. 
Therefore, a person having ordinary skill in the art could 
develop a suitable selective activation algorithm based on his 
level of skill coupled with an understanding of the teachings of 
the present invention. 

Amendment (12/15/05), pp. 24-26. Instead of addressing these arguments and presenting 
evidence why persons of ordinary skill in the art would not understand the meaning of 
Doppler broadening algorithm, positron lifetime algorithm, three-dimensional imaging 
algorithm, and selective activation algorithm within the scope of the claims, the examiner 
pronounced the arguments "not persuasive." Without any further explanation, he merely 
incorporated his previous rejections by reference into his Final Office Action. Final 
Office Action (2/14/05) at p. 13. In a similar situation where, in the examiner's answer, the 
examiner responded that the appellants' arguments were "not persuasive," the Board held 
that the examiner failed to satisfy his burden. See Ex parte Voordum, 2002 WL 1801 166, 
*4 (Bd. Pat. App. & Interf. 2002). 

As demonstrated above, one of ordinary skill in the art would well understand the 
claim language, "Doppler broadening," "positron lifetime," "selective activation," and 
"three-dimensional imaging algorithm." The examiner has failed to carry his burden to 
demonstrate otherwise. Claims 21-24, 39 and 40 depend from claim 20 and are allowable 
at least for the reasons that claim 20 is allowable. Claims 27, 28, and 30 depend from 
claim 26 and are allowable at least for the reasons that claim 26 is allowable. Claims 32 
and 33 depend from claim 31 and are allowable at least for the reasons that claim 31 is 
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allowable. Claim 35 depends from claim 34 and is allowable at least for the reasons that 
claim 34 is allowable. Claim 38 depends from claim 37 and is allowable at least for the 
reasons that claim 37 is allowable. 



The Board's previous determination of the sufficiency and enablement of the 
identical specification and the definiteness of the related method claims in the divisional 
case applies equally well here. The algorithms constituting claim limitations are described 
in the specification as well-known, which, the Board has previously held, would enable 
one of ordinary skill in the art to practice the claimed invention without undue 
experimentation. The examiner cannot sustain his prima facie burden by arguing that 
limitations not contained in the claims are not enabled. Similarly, because the algorithms 
are well-known, one of ordinary skill in the art would understand the scope of the claims 
including such algorithms as a limitation. Given the degree of skill possessed by persons 
of ordinary skill in the art, the claims are sufficiently enabled and definite to be patentable. 
The examiner's rejections should be reversed. 



CONCLUSION 



Respectfully submitted, 
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CLAIMS APPENDIX 



1-3. (Canceled). 
4. (Withdrawn). 
5-19. (Canceled). 

20. Non-destructive testing apparatus, comprising: 

a photon source, said photon source producing photons having a predetermined 
energy and directing the photons toward a specimen being tested, the photons from 
said photon source resulting in the creation of positrons within the specimen being 
tested; 

a detector positioned adjacent the specimen being tested, said detector 
producing raw data indicative of a positron annihilation event; and 

a data processing system operatively associated with said detector and said 
photon source, said data processing system operating in accordance with a normal 
activation/analysis process when a half-life of a selected positron emitter within the 
specimen being tested is greater than a predetermined half-life, said data processing 
system operating in accordance with a rapid activation/analysis process when a half- 
life of the selected positron emitter within the specimen being tested is less than the 
predetermined half-life, said data processing system, when operated in accordance 
with the rapid activation/analysis process, alternatively activating said photon source 
and detecting raw data indicative of a positron annihilation event, said data processing 
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system including a Doppler broadening algorithm, said Doppler broadening 
algorithm processing raw data indicative of a positron annihilation event to produce 
output data indicative of the presence or absence of a lattice defect in the specimen 
being tested. 

21. The non-destructive testing apparatus of claim 20, wherein said detector 
produces raw data indicative of a positron formation event, and wherein said data processing 
system includes a positron lifetime algorithm, said positron lifetime algorithm processing raw 
data indicative of a positron formation event to produce output data indicative of a changing 
presence or absence of a lattice defect. 

22. The non-destructive testing apparatus of claim 20, further comprising a second 
detector positioned adjacent the specimen being tested, said second detector producing raw 
data indicative of a positron formation event, wherein said data processing system includes a 
positron lifetime algorithm, said positron lifetime algorithm processing data indicative of a 
positron formation event to produce output data indicative of a changing presence or absence 
of a lattice defect. 

23. The non-destructive testing apparatus of claim 20, wherein said data 
processing system includes a selective activation algorithm, said selective activation 
algorithm responsive to a user input, said selective activation algorithm operating said photon 
source to produce photons having the predetermined energies in response to the user input. 
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24. The non-destructive testing apparatus of claim 20, wherein said data 
processing system includes a three-dimensional imaging algorithm, said three-dimensional 
imaging algorithm processing raw data indicative of a positron annihilation event to produce 
output data indicative of a location of the presence or absence of a lattice defect within the 
specimen being tested. 

25. (Canceled). 

26. Non-destructive testing apparatus, comprising: 

positron activation means for activating a positron emitter within a specimen 
being tested; 

detector means for detecting a positron annihilation event within the specimen 
being tested and for producing raw data indicative of the positron annihilation event; 

means for alternately activating the positron emitter within the specimen being 
tested and detecting a positron annihilation event; and 

data processing means operatively associated with said detector means, said 
data processing means processing raw data indicative of the positron annihilation 
event in accordance with a Doppler broadening algorithm to produce output data 
indicative of the presence or absence of a lattice defect in the specimen being tested. 

27. The non-destructive testing apparatus of claim 26, wherein said detector means 
detects a positron formation event and a positron annihilation event and produces raw data 
indicative of the positron formation event and the positron annihilation event, and wherein 
said data processing means processes raw data indicative of the positron formation event in 
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accordance with a positron lifetime algorithm to produce output data indicative of a changing 
presence or absence of a lattice defect. 

28. The non-destructive testing apparatus of claim 26, further comprising second 
detector means for detecting a positron formation event and for producing raw data indicative 
of the positron formation event, wherein said data processing means processes raw data 
indicative of the positron formation event in accordance with a positron lifetime algorithm to 
produce output data indicative of a changing presence or absence of a lattice defect. 

29. (Canceled). 

30. The non-destructive testing apparatus of claim 26, wherein said means for 
alternately activating the positron emitter within the specimen being tested and detecting a 
positron annihilation event comprises means for moving the specimen being tested between 
an activation position and a detection position. 

3 1 . Non-destructive testing apparatus, comprising: 

a photon source, said photon source producing photons having a predetermined 
energy and directing the photons toward a specimen being tested, the photons from 
said photon source resulting in the creation of positrons within the specimen being 
tested; 

a detector positioned adjacent the specimen being tested, said detector 
producing raw data related to a positron annihilation event; and 
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a Doppler broadening processor operatively associated with said detector and 
responsive to the raw data produced thereby, said Doppler broadening processor 
producing output data indicative of the presence or absence of a lattice defect in the 
specimen being tested. 

32. The non-destructive testing apparatus of claim 31, further comprising three- 
dimensional imaging apparatus operatively associated with said detector and responsive to the 
raw data produced thereby, said three-dimensional imaging apparatus producing output data 
indicative of a location of the presence or absence of a lattice defect within the specimen 
being tested. 

33. The non-destructive testing apparatus of claim 31, wherein said detector 
produces raw data that include data indicative of a positron formation event and data 
indicative of a positron annihilation event, said non-destructive testing apparatus further 
comprising a positron lifetime processor operatively associated with said detector and 
responsive to the raw data produced thereby, said positron lifetime processor producing 
output data indicative of the presence or absence of a lattice defect of the specimen being 
tested and indicative of a changing presence or absence of a lattice defect. 

34. Non-destructive testing apparatus, comprising: 

a photon source, said photon source producing photons having a predetermined 
energy and directing the photons toward a specimen being tested, the photons from 
said photon source resulting in the creation of positrons within the specimen being 
tested; 
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a detector positioned adjacent the specimen being tested, said detector 
producing raw data indicative of a positron formation event and a positron annihilation 
event; and 

a positron lifetime processor operatively associated with said detector and 
responsive to the raw data produced thereby, said positron lifetime processor 
producing output data indicative of a the presence or absence of a lattice defect in the 
specimen being tested and indicative of a changing presence or absence of a lattice 
defect. 

35. The non-destructive testing apparatus of claim 34, further comprising three- 
dimensional imaging apparatus operatively associated with said detector and responsive to the 
raw data produced thereby, said three-dimensional imaging apparatus producing output data 
indicative of a location of the presence or absence of a lattice defect within the specimen 
being tested. 

36. Non-destructive testing apparatus, comprising: 

a photon source, said photon source producing photons having a predetermined 
energy and directing the photons toward a specimen being tested, the photons from 
said photon source resulting in the creation of positrons within the specimen being 
tested; 

a detector positioned adjacent the specimen being tested, said detector 
producing raw data indicative of a positron formation event and a positron annihilation 
event; and 
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a data processing system operatively associated with said detector, said data 
processing system including: 

a Doppler broadening algorithm, said Doppler broadening algorithm 
processing raw data indicative of a positron annihilation event to produce output data 
indicative of a presence or absence of a lattice defect in the specimen being tested; 

a positron lifetime algorithm, said positron lifetime algorithm processing raw 
data indicative of a positron formation event to produce output data indicative of a 
changing presence or absence of a lattice defect; and 

a three-dimensional imaging algorithm, said three-dimensional imaging 
algorithm processing raw data indicative of a positron annihilation event to produce 
output data indicative of a location of the presence or absence of a lattice defect within 
the specimen being tested. 

37. Non-destructive testing apparatus, comprising: 

a photon source, said photon source producing photons having a predetermined 
energy and directing the photons toward a specimen being tested, the photons from 
said photon source resulting in the creation of positrons within the specimen being 
tested; 

a detector positioned adjacent the specimen being tested, said detector 
producing raw data indicative of a positron annihilation event; and 

a data processing system operatively associated with said detector and said 
photon source, said data processing system including: 

a Doppler broadening algorithm, said Doppler broadening algorithm 

processing raw data indicative of a positron annihilation event to produce 
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output data indicative of a presence or absence of a lattice defect in the 
specimen being tested; and 

a positron lifetime algorithm, said positron lifetime algorithm processing raw 
data indicative of a positron formation event to produce output data indicative 
of a changing presence or absence of a lattice defect, 
said data processing system operating in accordance with a normal activation/analysis 
process when a half-life of a selected positron emitter within the specimen being tested 
is greater than a predetermined half-life, said data processing system operating in 
accordance with a rapid activation/analysis process when a half-life of the selected 
positron emitter within the specimen being tested is less than the predetermined half- 
life, said data processing system, when operated in accordance with the rapid 
activation/analysis process, alternatively activating said photon source and detecting 
raw data indicative of a positron annihilation event, said data processing system using 
said Doppler broadening algorithm and said positron lifetime algorithm to process raw 
data indicative of a positron annihilation event and produce output data indicative of 
the presence or absence of a lattice defect in the specimen being tested. 



38. The non-destructive testing apparatus of claim 37, wherein said data 
processing system further comprises a three-dimensional imaging algorithm and wherein said 
data processing system uses said three-dimensional imaging algorithm to produce output data 
indicative of a location of the presence or absence of a lattice defect within the specimen 
being tested. 
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39. The non-destructive testing apparatus of claim 20, further comprising means 
for alternately moving the specimen between an activation position and a detection position, 
the activation position being adjacent said photon source, the detection position being 
adjacent said detector. 

40. The non-destructive testing apparatus of claim 20, further comprising means 
for alternately moving said photon source adjacent the specimen during the activation time 
and away from the specimen during the detection time and for alternately moving said 
detector adjacent the specimen during the detection time and away from the specimen during 
the activation time. 
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The opinion in support of the decision being entered today was not written 
for publication and is not binding precedent of the Board. 
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Before PATE, MCQUADE, and BAHR, Administrative Patent Judges . 
MCQUADE, Administrative Patent Judge . 

DECISION ON APPEAL 

Douglas W. Akers originally took this appeal from the final 
rejection (mailed June 16, 2003) of claims 1 and 8 through 39, 
all of the claims pending in the application. After the 
appellant's main and reply briefs (filed October 3, 2003 and 
January 13, 2004) and an examiner's answer (mailed December 1, 
2003) had been made of record, the examiner issued an Office 
action (mailed February 25, 2004) reopening prosecution and 
entering superseding rejections of all of the claims. Pursuant 
to 37 CFR § 1.193(b) (2) (ii), the appellant then submitted a 
request that the appeal be reinstated (filed April, 14, 2004) and 
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a supplemental brief (filed August 2, 2004). The examiner, 
implicitly granting the request, entered a second answer (mailed 
September 14, 2004) and forwarded the application to this Board 
for review of the new rejections of claims 1 and 8 through 39. 

THE INVENTION 

The invention relates to "methods . . . for performing non- 
destructive testing of materials using positron annihilation* 
(specification, page 1). Representative claim 1 reads as 
follows: 

1. A method, comprising: 

determining whether a specimen to be tested includes at 
least one positron emitter therein that will be activated in 
response to photon bombardment; 

selecting a positron emitter to be activated; 

determining a threshold photon energy required to activate 
the selected positron emitter; 

determining a half-life of the selected positron emitter; 

and 

when the half-life of the selected positron emitter is less 
than a selected half -life, then performing a rapid 
activation/analysis process, said rapid activation/analysis 
process comprising: 

activating for an activation time the selected positron 
emitter by bombarding the specimen with photons having energies 
at least as great as the threshold photon energy; 

detecting for a detection time gamma rays produced by 
annihilation of positrons with electrons in the specimen; and 

repeating said steps of activating for an activation time 
and detecting for a detection time until detecting a sufficient 
number of gamma rays to determine at least one material 
characteristic of said specimen; 
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when the half -life of the selected positron emitter is 
greater than or equal to the selected half-life, then performing 
a normal activation/analysis process, said normal 
activation/analysis process comprising: 

activating the selected positron emitter by bombarding the 
specimen with photons having energies at least as great as the 
threshold photon energy; and 

detecting gamma rays produced by annihilation of positrons 
with electrons in the specimen. 

THE REJECTIONS 

Claims 1 and 8 through 39 stand rejected under 35 U.S.C. 
§ 112, first paragraph, as being based on a specification which 
fails to comply with the enablement requirement. 

Claims 1 and 8 through 39 also stand rejected under 35 
U.S.C. § 112, second paragraph, as failing to particularly point 
out and distinctly claim the subject matter the appellant regards 
as the invention. 

Attention is directed to the main, reply and supplemental 
briefs and the second answer for the respective positions of the 
appellant and the examiner regarding the merits of these 
rejections. 1 



1 In the Of f ice action which reopened prosecution, claims 1 
and 8 through 10 additionally stood rejected under 35 U.S.C. 
§ 112, first paragraph, as being based on a specification which 
fails to comply with the written description requirement. Upon 
reconsideration, the examiner has withdrawn this rejection (see 
page 2 in the second answer) . Presumably, the examiner also has 
withdrawn any reasoning relating to the two remaining rejections 
which was set forth in the Office action but not restated in the 

(continued...) 
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DISCUSSION 

I. Preliminary matter 

In the supplemental brief (see pages 3, 4 and 17) and second 
answer (see pages 11 and 12), the appellant and the examiner 
debate the propriety of the Office action which reopened 
prosecution subsequent to appeal. As this matter is not directly 
connected with the merits of any issue involving a rejection of 
claims, it is reviewable by petition to the Director rather than 
by appeal to this Board ( see In re Henaehold . 440 F.2d 1395, 
1403-04, 169 USPQ 473, 479 (CCPA 1971), and hence will not be 
further addressed in this decision. 

II. The 35 U.S.C. § 112, first paragraph, rejection 

The dispositive issue with respect to the enablement 
requirement of 35 U.S.C. § 112, first paragraph, is whether the 
appellant's disclosure, considering the level of ordinary skill 
in the art as of the date of the application, would have enabled 
a person of such skill to make and use the claimed invention 
without undue experimentation. In re Strahilevitz . 668 F.2d 
1229, 1232, 212 USPQ 561, 563-64 (CCPA 1982). In calling into 



'(...continued) 

second answer ( see Ex parte Emm , 118 USPQ 180, 181 (Bd. App. 
1957)) . 
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question the enablement of the disclosure, the examiner has the 
initial burden of advancing acceptable reasoning inconsistent 
with enablement. I£. 

The examiner considers the appellant's disclosure to be non- 
enabling in four respects: 

. . . There is neither an adequate description nor 
enabling disclosure as to what is encompassed by the 
term, ''activating a positron emitter* . A ^positron 
emitter* is inherently already activated, i.e., it is 
radioactive. 

At best, the use of the term Activating a 
positron emitter 0 is superfluous. 

At worst, the term would imply that a positron 
emitter is either being transformed to another positron 
emitter or its energy level is further raised by the 
energy of the activating photon. There is no support 
in the specification for either one of these two 
alternatives [second answer, page 3J; 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner 
potential interferences in the data are accounted for 
in the analysis. For example, Claim 1 recites the step 
of bombarding the specimen with photons at least as 
great as the threshold photon energy required to 
activate the selected positron emitter and detecting 
gamma rays produced by annihilation of positrons with 
elections in the specimen. 

. . . [TJhere is no support as to how one would 
differentiate between the signals from the selected 
positron emitter and from the non-selected ones [second 
answer, pages 6 and 7]; 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner one 
should select an algorithm from a plurality of 
available algorithms, modify /manipulate the selected 
algorithm and evaluate the constants in the selected 
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algorithm to fit Appellant's situation. Appellant 
himself admits in the specification that there is not 
only one known algorithm but rather a plurality of 
known positron lifetime algorithms (e.g. see paragraph 
0058), Doppler broadening algorithms (paragraph 0057) 
and three-dimensional imaging algorithms (e.g., see 
paragraph 0060) . 

. . . There is no support as to how and in what 
manner one selects the specific algorithm to use for 
the two activation/ analysis processes, to evaluate 
requisite constants and to modify the selected 
algorithm to Appellant's situation [second answer, page 
7]; 

and 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner one 
can determine when the half-life of the selected 
positron emitter is less than a selected half-life 
(sic), where the composition of the specimen is unknown 
and it is this composition that is desired to be 
determined. The same lack of support exists for the 
determination of when the half-life of the selected 
positron emitter is more than a selected half-life 
(sic) [second answer, page 8]. 

A fair assessment of the appellant's disclosure shows that 
the references therein to '"activating" a positron emitter through 
photon bombardment of the specimen would be understood by the 
artisan as synonymous with producing* or 11 forming* the positron 
emitter from a stable precursor. For example, the specification 
states that "photons 16 from the photon source 12 produce such 
neutron-deficient isotopes . . . referred to herein in the 



alternative as 'positron emitters'* (page 7, enqphasis added), and 
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do not require the composition of the specimen to be unknown, and 
are not directed to a method for determining the composition of 
the specimen. 

Finally, while the appellant's disclosure does not convey 
much detail as to the various algorithms described therein, these 
algorithms admittedly are generally known in the art. The 
examiner has failed to advance any cogent reasoning as to why the 
disclosure would not have enabled a person having ordinary skill 
in the art to employ these algorithms, without undue 
experimentation, to the extent required by the appealed claims. 

Hence, the examiner's position that the appellant's 
disclosure, considering the level of ordinary skill in the art as 
of the date of the application, would not have enabled a person 
of such skill to make and use the claimed invention without undue 
experimentation is not well taken. Therefore, we shall not 
sustain the standing 35 U.S.C. § 112, first paragraph, rejection 
of 1 and 8 through 39. 

III. The 35 U.S.C. S 112, second paragraph, re-iectioq 

The second paragraph of 35 U.S.C. § 112 requires claims to 
set out and circumscribe a particular area with a reasonable 
degree of precision and particularity. In re Johnson . 558 F.2d 
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1008, 1015, 194 USPQ 187, 193 (CCPA 1977). In determining 
whether this standard is met, the definiteness of the language 
employed in the claims must be analyzed, not in a vacuum, but 
always in light of the teachings of the prior art and of the 
particular application disclosure as it would be interpreted by 
one possessing the ordinary level of skill in the pertinent art. 
Id. 

The examiner views the appealed claims to be indefinite for 
reasons essentially similar to those listed above in connection 
with the enablement rejection. In this regard, the examiner 
points to the claim limitations relating to (1) the activation of 
a positron emitter, (2) the positron lifetime, Doppler 
broadening, three-dimensional imaging and activation/analysis 
process algorithms and (3) the determination of the half-life of 
the selected positron emitter (see pages 9 through 11 in the 
second answer) . For the reasons explained previously, the 
limitations pertaining to the activation of a positron emitter 
would have been readily understood by the artisan when read in 
light of the underlying specification. Furthermore, the 
examiner's criticisms of the claim limitations pertaining to the 
algorithms and the half-life of the selected positron emitter 
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focus on the breadth of these limitations. It is well settled, 
however, that mere breadth does not equate to indef initeness . 
In re Miller . 441 F.2d 689, 169 USPQ 597, 600 (CCPA 1971). 
Accordingly, the examiner's contention that the appealed claims 
are indefinite is unpersuasive. 

Thus, we shall not sustain the standing 35 U.S.C. § 112, 
second paragraph, rejection of claims 1 and 8 through 39. 

SUMMARY 

The decision of the examiner to reject claims 1 and 8 
through 39 is reversed. 
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REVERSED 



WILLIAM F. PATE, III 
Administrative Patent Judge 






JOHN P. MCQUADE 
Administrative Patent Judge 
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Administrative Patent Judge 
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The reliability of . least squares analysis of pooitron annihilation lifetime data represented a. a sum or conrorctod exoonentials i. 
investigated w,th computer-stimulated test data. The method of analysis, which is available T . c«^«T -SSLT^! 
"SPLMOD- [R.R Vogel. Technical Repon EMBUDA08, European Mo£u«J? WcU^K«£v S^XS 
avoids direct determination of the Instrument resolution function by employing the decay curved. mf^n^J^'i 
well-known single Uf etime. ^performance of „. algorithm for ^L^Z^nZSxZSll Z^o^St^ 
by using several measures of the mlonnatioa content of the decay curve descriiw w^tW wv c'w^T V^, evaluated 
Annihilation Studies of Flcid.. ed. $. Shanna (WoHd £J^l^£SS 3T .L^v^^^^^ 

C mvI" 0 ^ k ^ rtSO,udon r wwion - 8m ^ u w the zero time channel and conlaminatJon of the reS^y «™ £ 
-ddiuon.1 hfenme exponents was evaluated. The method provide, excellent estimate, of lifetime par^et^Twhen^^e^l 



1. Introduction 

The analysis of positron annihilation lifetime data is 
an extremely difficult undertaking. The decay curve is 
usually represented by: 

c(/)= £«,x f .ex P (>x,o> (i) 
i-i 

where a, is the number of positrons utinihiin^ng with 
natural lifetime, X," 1 . Analysis of sums of exponentials 
is notoriously difficult, particularly when the number of 
components, «, is unknown. In general, even for data 
with arbitrarily small noise levels, there exist a large 
number of solutions, [a /t X„ «], that fit the data to 
within the noise level. These difficulties are com- 
pounded by the fact that in a positron annihilation 
lifetime experiment one cannot observe the function 
C(r) directly, only iu convolution with the Instrument 
resolution function. R(t): 

>(')"*<')*£(/). (2) 
The difficulties associated with the analysis of eq. (2) 
are therefore two-fold: the choice of an appropriate 
representation of C(t) and the o^teirnination of the 
function", R{t). 

0168-9002/89/503.50 <0 Elsevier Science Publishers B.V. 
(North-Holland Physics Publishing Division) 



The problems associated with analysis of data with 
eq. (1) are very well known. In general, we have no prior 
knowledge of the number of components, n, in eq. (I). 
If n is too large the solutions become very unstable to 
noise in the data. Because of this difficulty, it is com- 
mon to assume only two or three components in eq. (1). 
This has the effect of stabilizing the solution to noise 
and in some simple systems may provide a physically 
meaningful representation of the decay. However, in 
more complex systems setting n too small will result in 
an inadequate representation of the system, resulting in 
a loss of valuable information and a sex of a, and X, 
pairs that may have little or no physical meaning.' 
Because of the non-linearity of C(f ) in the X /t it is not 
possible to extend linear hypothesis tests to define an 
upper value for nl However, Provcncher (1) has de- 
scribed a hypothesis test which incorporates a correc- 
tion for non-linearity that has been found to be very 
effective, II is employed in the programs DISCRETE, 
which solves sums of exponentials and SPLMOD which 
solves more general sums of one parameter functions, 
including exponentials and convoluted exponentials. The 
use of eq. (1) assumes that a discrete sum of exponen- 
tials is a physically meaningful representation of the 
data. For many simple, homogeneous systems this is 
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indeed the case, however, for complex, heterogeneous 
systems and systems in which distinct structural do- 
mains persist for times that arc long relative to positron 
lifetimes, the lifetimes spectrum Is expected to be more 
complex and may consist of a continuous distribution 
of lifetimes over some interval For such systems an 
alternative approach is to replace the sum in eq. (1) by 
an integral (2-4]: 

C(t)-j"\*(\) exp(-Xr)dX=^{Xa(X)) ( (3) 

C(f) is simply the Laplace transform of X«(X), where 
a(X) is the annihilation rate probability density func- 
tion (pdf). The fraction of positrons annihilating with 
rates between X and X + dX is then given by or(X) dX 
(assuming <x(X) is appropriately normalized). In princi- 
ple, a(\) can be recovered by inverse Laplace transfor- 
mation of C(0- The problems of solution uniqueness 
and stability to noise associated with analysis, of eq. (1) 
also plague the inversion of eq. (3). However* recent 
advances in the solution of Fredholm integral equations 
provide a number of reliable inversion methods (4- 9 J. 

Schrader [2,3] has pioneered the application Of 
Laplace transform methods in the analysis of positron 
annihilation lifetime data and the application of infor- 
mation theory approaches to the determination of the 
information content of the experimental data. A par- 
ticularly important result of this work is the expression 
of the information content (Shannon Number) of the 
data in terms of spectrometer parameters (i.e. the num- 
ber Of channels, the channel width, fWHM and the 
total number of annihilation events observed). Not only 
does this knowledge define the complexity of the func- 
tion, £(0> tAat iS required to properly represent the 
data, but it also serves as a guide in the design of 
experiments to maximize the information content of the 
data. 

The accurate determination of the instrument resolu- 
tion function, /*(')» 16 * major problem in analysis of 
positron lifetime data. Probably the most widely used 
method is to determine R(r) from the prompt spectrum 
of a 60 Co source. However, the resolution function is 
very dependent on the distribution of energy deposited 
in the scintillators. As a result, the function determined 
with the 1.17 McV and 1:33 MeV coincident v-rays 
from 60 Co will not correctly represent the resolution 
function when data are collected with a 22 Na source. 
One approach, that attempts to overcome these difficul- 
ties, is to determine the resolution function from the 
lifetime data themselves. An analytical function, such as 
a sum of Gaussians or a modified Gaussian, is assumed 
for R{<) and its parameters are treated as variables to 
be determined along with the lifetime parameters of 
interest [10,11]. However, there are several fundamental 
objections to this approach. The addition of further 
parameters to be determined by the analysis can severely 



compound the problem of solution uniqueness. In ad- 
dition, the approach requires the choice of an analytical 
function to represent the resolution properties of the 
instrument. An incorrect choice will, of course, transmit 
errors to the lifetime estimates. 

The resolution function may also be determined in a 
separate experiment from the decay curve of a standard 
sample with a known single lifetime, such as a well-an- 
nealed metal [12]. In principle, this method can yield 
very good estimates of the instrument resolution func- 
tion. However, in practice, annihilation in the source 
and positron trapping at the surface can give rise to 
additional components. An alternative approach is to 
employ a ^"Bi source. ^Bi emits two y-rays separated 
by a time interval of 186 ps, the energies of which are 
almost Identical to those employed in a positron life- 
time experiment with a 22 Na source. Although the 
method does not introduce unwanted source compo- 
nents, it is not yet clear if the energies of the y-rays 
emitted by "'Bi are sufficiently similar to those occur- 
ring with a Ni that the types of error associated with 
the use of 60 Co are entirely avoided. 

At present, the use of standard samples with 23 Na 
sources appears to offer the best hope for reliable 
deierrnination of the instrument resolution function. We 
have been interested in the use of a deconvolution 
method due to Gauduchon and V/ahl [13] and 
Wijnaendts van Resandt ct al. [14] for this purpose. A 
pretinunary report of its use in the analysis of pcsiiron 
annihilation lifetime data with numerical Laplace inver- 
sion methods has appeared [4] but no tests of its perfor- 
mance have been presented. Our purpose in this paper 
is to explore the reliability of this deconvolution method 
in the analysis of positron annihilation lifetime data. 
The efficiency of the algorithm for extracting lifetime 
information is evaluated by using several measures of 
the information content of the decay curve described by 
Schrader and Usmar [3]. In this work, C(t) is repre- 
sented by a sum of exponential terms as in eq* (1) and 
the analysis is performed with the FORTRAN program, 
SPLMOD, developed by Vogel [15]. Although eq, (1) is 
more restrictive than eq. (3), it can be a valid represen- 
tation of annihilation lifetimes in simple, homogeneous 
systems and is considerably simpler for our purpose 
because the measure of information content of the data, 
the Shannon Number, is directly related to the number 
of components in the lifetime spectrum. 



2. Deconvolution of the instrument resolution function 

Deconvolution of the instrument resolution function 
follows the method of Gauduchon and Wahl [13] and 
Wynacndts van Resandt et aL [14], which avoids the 
direct determination of R(t) by measuring, >>,(')• toe 
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decay curve of a reference material with a well-known 
single lifetime, X,: ^ 

y,(t) - A (0 •<;<<), (4) 

where 

<UO-*r*r«p(-X T /). (5) 
Convoluting eq. (2) with C r (r) we obtain: 
Q('W(0 C(r). (6) 

Here we assume that y t (i) and ^(r) have ihe same 
resolution function, R(t). This would be true if the 
reference and sample decay curves could be measured 
simultaneously, but this is not possible in a positron 
annihilation lifetime experiment. Any instrument drift 
occurring between measurement of the sample and ref- 
erence decay curves will therefore introduce error into 
the lifetime estimates. This, of course, is true of all 
methods that involve a separate measurement of the 
resolution function of the instrument and can only be 
avoided by improving the stability of the instrument. 
Eq. (6) can be rewritten as: 

^r{«p(-V))v(0->i(O't(O. (?) 

which on taking Laplace transforms gives: 

$jHp)-Mp)£{p). (8) 

where ihe tilde denotes Laplace transformation and p 
is the Laplace variable. Eq. (8) is solved for Hp)' 

Hp)-^7 , K 1 (KMp)^p)+pMp)^(p)). 

(°) 

which on inverting Laplace transforms gives: 

• c (0 + ->v(')-K'(0 + c(o)0(f)]). 

(10) 

If CO) is represented by eq. (1), eq. (10) becomes: 
•v(') 

*" £?j^(*< f > + <X f -A,)*(0.«p<-*/0}. 

<«) 

If A f is known and a r is absorbed into the cr„ eq. (11) 
must be solved for the 2n parameters a f and X/. If X 
is not known, there arc 2« 4- 1 parameters. - ' 



The current version of SPLMOD (version 3) distributed 
by Vogel [15] can perform the fonowing types of analy- 
sis of relevance to Positron Annihilation Lifetime (PAL) 
spectroscopy; 

(a) Solution of eq. (12) when the resolution function is 
known; U with /(X„ r,) - A(f A ).exp(-A,r,) and 

- 1 (to account for the background). 

(b) Solution of eq. (12) using a reference material with a 
known single annihilation lifetime* A r . i. e . with 
/(*/» '*)-<*, -A,)jr r (r A )«exp(-A,r + and 

(c) Solution of eq. (12) using a reference maieriaJ with a 
single annihilation lifetime which need not be 
known; i.e. with f(\, t t M ) -M'*)»expf>A,r,), 

*.(f A )-3and **('*) -*('*)- 
We have employed option (b) in our simulations be- 
cause it provides more stable solution* than option (c). 

SPLMOD performs a weighted least squares analysis 
in which /(A /( t k ) is approximated by a sum of normal- 
ized cubic B-spUnes. The program employs Gauss- 
Newton iniiunuzauon, with an extensive search of 
parameter space to determine the initial parameter val- 
ues. The convolution, »«*>*«9(-X,f). is calculated 
by using a recursive formula of Grinvald and Steinberg 
[161: 

/(*/. <* + i) = {/(X„ f*)+O.Mov(r*)) 

Xcxp(-*X,nr) + 05Aty r (t M ^), (13) 

where At is the channel width of the mulU-cJuuinel 
analyzer in which the data are collected. SPLMOD 
performs several analyses in which n in eq. (1) takes the 
values: 1, 2, 3,... n mw% where nfM is set by the user. 
SPLMOD supplies a number of statistical measures, 
including Fisher probabilities and random runs profcn 
abilities, to aid in the choice of an appropriate solution. 
Measures of the information content of the data may 
also be employed to define an appropriate solution [3], 
Note that it is not necessary to independently define 
the zero time point of the spectrometer. This informa- 
lion is already contained in the reference decay curve. 
However, there is an option in SPLMOD to correct for 
small changes in the zero time point due to instrument 
drift occurring between measurement of sample and 
reference decay curves. 



3. Analysis of sums of convoluted exponentials with 
"SPLMOD" 

SPLMOD analyzes data that can be represented as a 
sum of one parameter functions: 



(12) 



4. Simulation of decay curves 

Positron annihilation decay curves have been simu- 
lated by convolving eq. (l) with a Gaussian resolution 
function with FWHM equal to 2ov1nT : 



A(0--4=- exp-(r/o) 2 . 



(14) 
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xerfc(X,er/2 -</*), (\$) 
where eric is the compliment of the error function: 
erfc(x)-l -crf(x) = l - -p T exp(-/ 2 )dr. (16) 

Conversion of the continuous function. to a table 
of channel counts, was performed by numerical 
integration. A background and appropriate Poisson 
noise were added to each channel. The reference data 
were simulated in the same way except that no back- 
ground was added. The reference and sample spectra 
were calculated over the same range of channels in each 
case. 



S. Information content of the simulated data 

Schradcr and Usmar [3] have given a detailed discus- 
sion of the determination of the information content 
(Shannon Number) of decay curves and have provided 
an approximate expression for the Shannon Number in 
terms of instrument setting of the spectrometer. The 
Shannon number, in the context of PAL spectroscopy, 
is a measure of the number of resolvable lifetimes that 
can be geometrically fitted between the largest and 
smallest resolvable lifetimes. This measure is an ap- 
proximation because it depends on ones definition of 
the relative scatter in the decay curve and on the values 
that are chosen for the largest and smallest resolvable 
lifetimes. Schrader and Usmar [3] define the Shannon 
Number, F> as: 

\n(a x tf^At/ai FWHM), (17) 

where a y and a 2 are constants, FWHM is the full 
width at half the peak height of the resolution function, 
N ch is the number or channels to the right of the 
zero-time channel, /Y 0 , employed in the analysis (a 
maximum value for is defined by the number of 
channels for which the count is greater than the square- 
root of the background count) and is defined by; 

co^, = 2H- 1 ln(y^T/n), (18) 

ij is a measure of the effective relative noise in the decay 
curve and is defined as: 



E 



yV 1 



(19) 



u mftX has a special significance in the eigenfunclion 
expansion of Laplace transforms which we need not be 
concerned with here [3,5]. Eq. (18) is derived by assum- 
ing that the shortest resolvable lifetime is equal to 



a a FWHM and the longest resolvable lifetime is equal to 
a x NjH&t. A resolution factor, p, is defined as the ratio 
of resolvable lifetimes and it follows that F is the 
number of resolvable lifetimes that can be geometrically 
fitted between c 2 FWHM and 4\N<*fit> The constants 
fl| and a 2 need to be determined. In practice, the 
shortest resolvable lifetime is probably less than FWHM 
and the longest resolvable lifetime is longer than N A At, 
so a a < 1 and a, > 1. Schrader and Usmar [3] suggest 
a x «• 3 and a 2 ■» \ as reasonable values when a a FWHM 
> 34 r p but they suggest replacing tf 2 FWHM by 34/ if 
a 2 FWHM<3dr. 

The resolution factor is given by: 

/>-exp(Tr/Wnu*) (20) 

and is therefore dependent on the definition of the 
effective relative noise of the decay curve (eq. (19)). 
Estimates of a u a 2 and the correctness of cq. (19) are 
discussed below. 



6. Analysis of simulated data 

Table 1 lists the results of an analysts of simulated 
PAL data together with the information parameters for 
the decay curve and will serve as a point of reference for 
our subsequent discussions. The simulated date contain 
three components. The reference has a single lifetime of 
125 ps. A total of 800 channels were analyzed, starting 
30 channels (0.75 ns) to the left of the peak. A plot of 



Table 1 

Calculated information eon tent of simulated data and results 



of its analysis with 5PLMOD 


FWHM 
M 

Reference lifetime 
Total events 


290 ps 
25 ps 

125 ps 
10* 


p 2.6 
4 0.015 
«r»« 3.25 
maximum 442 




F 


Resolvable lifeUmes (ns) 
min. max. 


1 1 

2 0.5 

3 0.33 


3.8 
52 
6.0 


0.29 11 
0.15 22 
0.1 33 







Lifetimes (nS) 




Simulated 


0.125 


0.400 


2.000 


Fined 


0.124 


0.397 


1.997 






Alpha 




Simulated 


0.130 


0.470 


0.400 


Fitted 


0.125 


0.474 


0.401 






Background 




Simulated 




100 




Fitted 




100.5 
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(a 2 - 0.25 to 0.20). Beyond ihis value, SPLMOD is still 
able lo identify three components but the parameter 
values become very poor. Increasing the total count to 
10 7 does not improve these estimates. At a FWHM of 
290 p$, a value that is typical of modern PAL spec- 
trometers employing plastic scintillators, the shortest 
lifetime that can be resolved by SPLMOD is approxi- 
mately 50 ps (i.e. ^ = 0.17). An estimate of a y was 



of posliron-anhihUalicm UfaUme data 

determined by varying the number of channels em- 
ployed in the analysis. In all cases the start channel was 
sufficiently far to the left of the peak (-0.74 ne) to 
sample the baseline. With ^ « 1940, the longest life- 
time that can be resolved is approximately 100 ns 
(*i-4). 

We have examined the ability oi the program to 
resolve lifetime components and were particularly inter- 



Table 3 



The resolution of lifetime components with SPLMOD. Variations in the effective relative wise and resolution factor were generated 
by changing the total number of events. Simulation parameters were; FWHM -290 ps. At -123, reference lifetime -125 w and 
background — 100 



Simulated 



Fitted 



Fitted 



Fitted 



Total events 
Relative noise 

Resolution factor 
f = 4, fl a -0.25) 
Maximum 
Lifetime 1 
<ns) 2 

3 

4 

Alpha 1 
2 
3 
4 



0.125 
0.400 
1.280 
4.096 

0.250 
0.250 
0,250 
0.250 



10 fl 
0.024 
3.0 
2.9 
6.1 
959 
0.122 
0.386 
1.213 
3.947 

0.24S 
0.247 
0.248 
0.260 



10* 
0.054 
2.4 
3.6 
4.5 
451 
0.113 
0.313 
0.895 
3J38 

0.223 
0.210 
0.258 
0J10 



10 4 

0.074 

2.4 

4.1 

2.7 
64 

0.124 
0.609 
2.766 

0-341 
0.345 
0.314 



Tabic 4 



The effect of systematic errort on lifetime parameters estimated by SPLMOD. Simulation parameters were FWHM = 290 pi, 
Ai « 12.5 ps, background « 100, total events = 10 7 and reference lifetime -125 ps 





Simulated 


Incorrect reference FWHM 


Zero channel 




data 


285 ps 


295 ps 


shift 


Lifetime 1 


0.125 


0.127 


0.138 


0.134 


(ns) 2 


0.40O 


0.409 


0.410 


0.409 




2.000 


2.004 


2.004 


2-005 






0.003 


100.0 


0.002 


Alpha 1 


0.130 


0.663 


0.148 


0.247 




0.470 


0.171 


0M9 


0.402 




0.400 


0.150 


0J95 


0J51 






0.016 


0.008 


0.001 


Background 


100 


\O0A 


91.1 


99*9 



Contaminant % 



data 



Con Lamination of reference by additional 
component with lifetime ~ 450 ps 

5 10 



15 



Lifetime 
(ns) 

Alpha 



Background 



0.125 
0.400 
2.000 
0.130 
0.470 
0.400 
100 



0.122 
0.382 
2.001 
0.131 
0.476 
0.393 
100.3 



0.120 
0.365 
2.003 
0.135 
0.474 
0J91 
1001 



0.118 
0.347 
2.005 
0.137 
0.477 
0J86 
100.1 
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cstcd in determining if the value of the resolution factor 
calculated according to eq. (20) correctly predated the 
point at which SPLMOD Med to resoNe lifetime com. 
ponanis. An example of the results for such a test is 
shown in table 3. The simulated data consisted of four 
components with successive lifetimes increasing by a 
factor of 3.2, When the total number of events is 10*, 
the calculated resolution factor is 2.9 and the Shannon 
number, is 6.1 SPLMOD is able to resolve the four 
components and provides reasonably good estimates of 
the lifetime parameters. When the total number of 
events is reduced to 10 J , the effective relative noise is 
doubled, the resolution factor increases to 3.6 and the 
Shannon number decreases to 4 J. SPLMOD Is still able 
to resolve four components but the lifetime parameter 
values become poor. When the total number of events is 
reduced to 10*. the relative noise increases to 0.074, the 



resolution factor Increases to 4.1 and the Shannon 
Number decreases to 2.7. SFlMOD is unable to resolve 
the four components. Eq. (19) thus appears to provide a 
good estimate of the effective relative noise in the data 
as measured by the predictive power of the calculated 
resolution factor. 



7. Systematic errors 

We have examined the sensitivity of the SPLMOD 
algorithm to various systematic errors in the data. The 
most important sources of systematic error arise from: 
(i) variations in the resolution function for the reference 
and sample data, (ii) differences in the zero time point 
for the reference and sample, and (iii) contaminatioo of 
the reference with additional lifetime components. 
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Fig. 2. Effect of errors in the FWHM of ihe reference decay curve on the weighted residuals. The lifetime parameters are listed in 
table 4. The FWHM of the simulated simple decay curve was 290 ps. Sample data was analyzed using reference decay curves with 

FWHM of (a) 285 ps and (b) 295 pa. 
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Rg. 3. Effect of a shift in ihc zero ilroc channel on the weighted residual*. The lifetime parameter* are listed In table 4. The reference 
decay curve was shifted by - 2.5 ps with respect to the sample decay curve. 



The fitted lifetime parameters are extremely sensitive 
to variations in the FWHM of the resolution function. 
The lifetime parameters obtained from an analysis of a 
three component sample (FWHM = 290 ps) with refer- 
ence decays having incorrect FWHM values (285 and 
295 ps) arc given in tabic 4. The lifetime parameters 
determined by SPLMOD when the FWHM for the 
reference is less than thai employed for the sample are 
very poor. The intensities of all components are grossly 
in error. In addition, a spurious short-lived component 
appears. When the FWHM of the reference is made 
larger than that employed for the sample, a spurious 
long-lived component appears (100 ns). Estimates of the 
other lifetime parameters are in error, but arc not nearly 
as poor as those determined with the smaller reference 
FWHM. Plots of the residuals are shown In fig. 2. The 
fits to the data are poor about the peak. 

A shift in the rero channel of the reference with 
respect to the sample data also generates poor estimates 
of the lifetime parameters and an additional spurious 
short-lived component. The example given in table 4 
and fig. 3 was Obtained by shifting the reference decay 
curve by -2.5 ps with respect to the sample decay 
curve. This shift is a fifth .of the channel width em- 
ployed. SPLMOD can correct for small shifts in the 
zero time channel by assuming that y(r + c) has been 
used instead of v(r) and then expanding eq. (11) in a 
Taylor series about y(f) to the first order term. This 
gives an additional linear term in eq. (12). However, in 
this case the correction did little to improve the fit to 
the data. 

The analysis assumes the reference to have a single 
lifetime component, an Ideal which is unlikely to be 
exact for real systems. Annihilation in the source and 
trapping at the surface of the reference material un- 
doubtedly introduce additional lifetime components. 



The effect of introducing an additional component with 
a lifetime of 450 ps into the reference decay curve is 
shown in table 4. The estimated lifetime parameters are 
relatively insensitive to the presence of the contaminant, 
even when it represents 15% of the reference decay. The 
largest error appears in the sample component which 
has a lifetime closest to the contaminant lifetime. 



8. Conclusions 

The results presented in the previous sections suggest 
that the deconvolution and least squares analysis meth- 
ods employed by SPLMOD can provide excellent esti- 
mates of positron annihilation lifetime parameters for 
systems that can be represented by a discrexe sum of 
exponentials. As we have mentioned, this representation 
is too restrictive for many systems of interest in which 
distributions of positron lifetimes are expected. For 
such systems numerical Laplace inversion methods are 
obviously to be preferred. The results of our tests of a 
modified version of the program, CONT1N (3,5,6), which 
solves integral equations with convoluted kernels will be 
published shortly. It is nevertheless clear from the pre- 
sent work, that estimates of the Lifetime parameters are 
critically dependent on accurate knowledge of the in- 
strument resolution function. We doubt that the decon- 
volution method employed here is any more sensitive to 
errors in the resolution function than other algorithms 
that have been described in the literature or that are 
currently in use. For example, lifetimes reported using 
"Co prompt spectra to measure the instrument resolu- 
tion function are undoubtedly in error since it is known 
that this method provides a poor estimate of the true 
resolution function occurring with 22 Na sources [3J. 
There is little doubt that these problems will become 



l 
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even . more severe as we attempt to determine^continu- 
ous lifetime distribution* with numerical Laplace inver- 
sion techniques. 



Acknowledgement* 

We are grateful to Dr. Robert Vogel for making a 
copy of SPLMOD available. This work was supported 
by grants from the National Science Foundation (DMB 
85-18941), Research Corporation (10660), Kent State 
University Foundation and the National Institutes of 
Health Biomedical Research Support Grant. 



References 

|1) S.W. Provencher, J. Cham. Phys. 64 (1976) 277Z 
(2] D.M. Schroder, in: Positron Annihilation, ed*. P.O. Cole- 
man, S.C Sharma and L.M. Diana (North-Holland, 
Amsterdam, 1982) p. 91 2. 
[3] D.M Schrader and S.G. Usmar, in: Ppsitron Annihilation 



Studies of Fluids, cd. S.C Sharma (WoHd Scienuflc, 

Singapore, 1988). 
[4J R.B. Gregory and A. Procyk, in: Positron Annihilation 

Studies of FtuJds, ed. S.C Sharma (World Scientific 

Singapore. 1988). 
J5] J.O. McWhirler and E.R. Pike. J. Phys. A 11 (1978) 1729. 
[6] S.W. Provencher, Comput. Phys. Commun. 27 (1982) 213. 
|7J S.W. Provcncfaer, CompuL Phys. Commun. 27 (1982) 229. 
I8J A.K. Uvcssy and J. SHMng, AcU CrysL Am (1985) U3. 
[9] AJC Livesey and J.C. Brochon. Biophys. J. 52 (1987) 693. 
nOj P. Korkegaard, M. Eldrup, O.e. Mogensen and NJ. Peder- 

sen, CompuL Phys. Commun. 23 (1981) 307. 
[11] tt Sornuum, P. Kindl and W. Puff, NucL Instr. and Melh 

206 (1983) 203. 

fl2] L.C Smedskjaer and MJ. Fhiss, Meih. Exptl Phys. 21 
G983) 77. 

(13) P. Oauduchon and P. Warn, Biophys. Chem. 8 (1978) 87. 

(14) R.W. Wijnacndts van Resandt, R.H- VcgeJ and S.W. 
Provenchcr, Rev. ScL Instr. 53 (1982) 1392. 

[15] R.H. Vogrf, SPLMOD Users Manual, Technical Report 
EMBL-DA08, European Molecular Biology Laboratory, 
Heidelberg, FRO (1986). 

[16] A. Orinvald and \3» Steinberg, Anal. Biochem. 59 (1974) 
583. 



8605524 



Schaffer, James Peter 

DECONVOLUTED DOPPLER BROADENED POSITRON ANNIHILATION 
SPECTROSCOPY: CHARACTERIZATION OF DEFECTS IN ALUMINUM 

Duke University Ph.D. 1985 

University 
Microfilms 

international 300 N. Zeeb Road, Ann Arbor, Ml 48106 

Copyright 1985 

by 

Schaffer, James Peter 
All Rights Reserved 



i 



THIS PAGE BIANKWSTO) 



This is an authorized facsimile, made from the microfilm 
master copy of the original dissertation or master thesis 
published by UMI. 

The bibliographic information for this thesis is contained 
in UMFs Dissertation Abstracts database, the only 
central source for accessing almost every doctoral 
dissertation accepted in North America since 1861. 



300 North Zeeb Road 
P.O. Box 1346 
Ann Arbor. Michigan 48106-1346 USA 

800.521 .0600 734.761 .4700 
web www.il.proquest.com 



UMI 



Dissertation 
Services 



From:Pro£xuest 



COMPANY 



Printed in 2005 by digital xerographic process 
on acid-free paper 



THIS PAGE BUWKlMs™ 



DECONVOLUTEO DOPPLER BROADENED POSITRON ANNIHILATION 
SPECTROSCOPY: CHARACTERIZATION OF 
DEFECTS IN ALUMINUM 



Figure 14. A Comparison of the Peak Region of the Raw Data and the 

Deconvolved Data for a Quenched Aluminum Sample. 70 

Figure lb. A Plot of the Error Parameter vs. Quench Temperature. 72 

Figure 16. A Plot of R t Values for the Raw Data vs. Quench 

Temperature. 80 

Figure 17. A Plot of Values as Calculated from Average Lineshape 
Parameter Values for the Raw Data vs. Quench 
Temperature. 32 

Figure 18. A Plot of Kj Values for the Deconvoluted Data, Using All 

Retained Coefficients, vs. Quench Temperature. 83 

Figure 19. A Comparison of the Raw Data for an Aluminum Sample 
Quenched from 600°C and the Same Sample after a High 
Temperature Anneal. 84 

Figure 20. A Comparison of the Tail Region of the Raw Data for an 
Aluminum Sample Quenched from 600°C and the Same Sample 
after a High Temperature Anneal. 85 

Figure 21. A Comparison of the Deconvoluted Data for an Aluminum 
Sample Quenched from 600°C and the Same Sample After a 
High Temperature Anneal. 87 

Figure 22. A Comparison of the Tail Region of the Deconvoluted Data 
for an Aluminum Sample Quenched From 600°C and the Same 
Sample After a High Temperature Anneal. 88 

Figure 23. A Comparison of the Tail Region of the Same Samples as 
in Figure 22, but this time Both Data Sets were 
Deconvoluted Using the Same Number of Coefficients. 92 

Figure 24. A Plot of Rj Values for the Deconvoluted Data Using a 

Constant Number of Reconstruction Coefficients. 93 



LIST OF FIGURES 



Page 



Figure 1. A Schematic Diagram Showing the Effects of the Electron- 
Positron Pair Momentum on the Energy and Direction of the 

Figure 2. f^j^Ijrtn^ Function for the *p P ler 

Figure 3. A Schematic Diagram of the Quenching Apparatus. 
Figure 4. ^Schematic ojjgr. Showing the Relationship Between 

Figure b. A Plot Comparing the Power Spectrum of the Annihilation 
Energy Distribution Before and After Deconvol SJion! 

Figure 6. A Semi -Log Plot of the Number of Reconstruction 

Coefficients Retained by the FFT/PS Algorithm as a 
Function of the Number of Counts in the Peak. 

Figure 7. A Comparison of the Raw Data and the Deconvoluted Dala 

for an Irradiated Sapphire Sample. 49 

Figure 8. J Comparison of the Raw Data. h(x), with the Function 
Obtained by the Convolution of f(x) with r(x). 

Figure 9. a Comparison of the Resolution Function Before 

JnH°A^ lut i on * Af ? er Dewnvolution Using 70 Coefficients 
and After Deconvolution Using 243 Coefficients. 

Figure 10. A Plot of the Deconvoluted Lineshape Parameter xd) as a 
Function of quench Temperature. «™meier *d) as a 

Figure 11. £ Plot of th« Raw Data Lineshape Parameter x(l) as a 
(-unction of Quench Temperature. 

Figure 12. A Comparison of the Shape of Resolution Functions 
Collected at Different Times During a Series of 
Experiments. 



21 
27 

30 

46 

47 



bO 

54 
58 
62 

64 



Figure 13. A Graphical Interpretation of the Sensitivity and Error 

Lineshape Parameters. * or 6? 



F. A LISTING OF THE FORTRAN FFT/PS DECONVOLUTION PROGRAM 117 

G. RELATIONSHIP BETWEEN C v AND C d 141 

H. THE UTILITY OF "R" PARAMETER ANALYSIS 143 
REFERENCES 145 



TABLE OF CONTENTS 



ABSTRACT 

ACKNOWLEDGEMENTS 

LIST OF FIGURES 

CHAPTERS 

I. INTRODUCTION 

II. EXPERIMENTAL PROCEDURES 

The Uoppler PAS System, 17 
Specimen Preparation, 23 

III. DECONVOLUTION OF DOPPER PAS MEASUREMENTS 

The Mathematics of Decon volution, 29 
Previous Deconvolution Techniques, 33 
The FFT/Power Spectrum Deconvolution Method, 39 
Test of the Method, 44 

IV. RESULTS AND DISCUSSION 

The Utility of the FFT/PS Deconvolution Algorithm, 52 
The Separation of Multiple Defect Types, 74 

V. CONCLUSIONS AND FUTURE WORK 

APPENDICES 

A. POSITRON THERMAL I ZAT I ON TIMES 

B. THE EQUIVALENCE OF DOPPLER AND ANGULAR CORRELATION 

MEASUREMENTS 

C. CHEMICAL ANALYSIS OF THE ALUMINUM 

D. THE UTILITY OF FOURIER ANALYSIS IN DECONVOLUTION 

E. THE MATHEMATICAL DIFFICULTIES ASSOCIATED WITH THE 

INCORRECTLY POSED INTEGRAL EQUATION 



(v) 



ABSTRACT 



Positron Annihilation Spectroscopy (PAS) has become a useful tool for 
the nondestructive investigation of the defect structure of materials. Of the 
existing PAS measurement techniques Ooppler broadening appears to offer 
several practical advantages. These advantages include low radioactive source 
strength requirements, comparatively short experimental times, and minimal 
geometrical restrictions. Unfortunately, the technique suffers from some 
inherent disadvantages. These include comparatively poor instrumental energy 
resolution and a limited ability to distinguish between multiple defect types. 

In an attempt to overcome the difficulties associated with the use of 
Ooppler broadened PAS the Fast Fourier Transform/Power Spectrum (FFT/PS) de- 
convolution algorithm has been developed in this work. This algorithm is com- 
putationally efficient, nonlteratlve, and requires minimal a priori Tnforma- 
tion. The utility of the FFT/PS deconvolution algorithm is demonstrated for 
high purity aluminum samples containing known concentrations of vacancies or 
prismatic dislocation loops. It is shown that the FFT/PS procedure is capable 
of automatically correcting for accidental and/or intentional changes in the 
instrument resolution function. In addition, based on the limited data anal- 
yzed in this work, 1t appears that both the sensitivity and the accuracy of 
the energy distribution Uneshape parameters are improved by deconvolution. 
Finally, it is demonstrated that if the filtering error introduced by the 
FFT/PS procedure Is controlled, then It is possible to distinguish between 
samples containing predominately vacancy type defects and those with prismatic 
dislocation loops as the dominant defect type. 
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CHAiTER I 
INTRODUCTION 



Over 50 years ago, Dirac postulated the existence of an atomic particle 
with a mass equal to that of an electron, but which is positively charged. In 
1932 these positively charged electrons, referred to as positrons, were first 
observed by C. 0. Anderson as a by-product of cosmic ray showers. Since that 
time positron research has been pursued 1n such varied disciplines as high en- 
eryy nuclear physics, astrophysics, and materials science. It is within the 
broad context of materials science, and more specifically in the area of non- 
destructive testing (NDT) of aluminum, that the interaction of positrons with 
matter will be examined in this work. Before going on to describe the use of 
positrons as an NDT tool in the investigation of the defect structure of 
metals, a brief review of positron physics will be presented. A comprehensive 
treatment of this subject has been written by West (1974). 

As postulated by Dirac In 1935, positrons and electrons are complemen- 
tary states of matter such that when the two particles encounter one another 
they will annihilate each other. Their kinetic and rest mass energies will be 
dissipated as the energy of one, two, or three gamma rays 1n such a way as to 
conserve momentum. It can be shown that the ratio of the two-photon to three- 
photon annihilation event is 372:1 (Stevens 1980), while the ratio of the two- 
photon to one-photon events is ~2.6xl0 6 :l, (West 1974). With these ratios as 
justification the remainder of this work will consider only the annihilation 
event yielding two photons. If both the electron and positron were at rest, a 
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positrons in vacancy type lattice defects. This effect was first observed by 
MacKenzie et al_. and Dekhtyar et aK 1n 1964. Schultz etal_. (1981) describe 
the changes in the annihilation characteristics due to the localization of a 
positron In a raonovacancy type defect as follows: 1) Because the localized 
positron has a smaller overlap of its wave function with that of the high ener- 
gy core electrons than 1t does when diffusing through the nearly perfect lat- 
tice, the annihilation rate with high energy core electrons decreases. 2) The 
density of conduction electrons in the vacancy is also decreased resulting in a 
decreased annihilation rate with conduction electrons. However, the percentage 
reduction is less than in the case of core electrons so that the net result is 
an Inc^ase in the fraction of annihilation events involving low energy conduc- 
tion electrons, thus resulting in a decrease 1n the width of the Doppler broad- 
ened annihilation radiation distribution. 3) The potential for electrons in 
the vacancy is higher than normal resulting in the conduction electrons travel- 
ling more slowly. This in turn also decreases the observed Doppler broadening. 
4) A second order effect is that the overlap of the positron wave function with 
weakly bound outer-core electrons is reduced less than is the overlap with 
inner-core electrons. This too leads to less broadening of the observed 511 
KeV line, b) The uncertainty principle requires that localization of the posi- 
tron be accompanied by an increased spread in its momentum distribution (and a 
corresponding decrease in its de Broglie wavelength). As a result, the contri- 
bution of the positron to the center-of-mass velocity Increases corresponding 
to a broadening of the annihilation energy spectrum (Schultz et aU 1981). The 
first four effects lead to a narrowing of the annihilation energy spectrum 
while the fifth effect tends to broaden the spectrum. The combination of these 
effects results in more sharply peaked (narrower) Doppler spectra for metals 
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As discussed by Siegel (1980) 1t is important to note that the posi- 
tron is not equally sensitive to all lattice defects. It 1s sensitive to 
defect regions which can be characterized as having a lower-than-average elec- 
tron density, for example, vacancies, small vacancy clusters, voids, and dis- 
location lines or joys. On the other hand, positrons are relatively Insensi- 
tive to defect regions with a higher-than-average electron density such as 
interstitials. Indeed this defect specificity is one of the advantages PAS 
measurements offer as compared with other NOT techniques such as residual 
electron resistivity (Siegel 1980). 

In addition, positrons are capable of probing the entire depth of the 
valence band (MacKenzie et a1_. 1964) which is in contrast with techniques such 
as electrical and thermal conductivity, thermoelectric power, Hall coef- 
ficient, magnetic susceptibility, and Knight shift, all of which yield infor- 
mation restricted to the electrons near the Fermi surface. 

PAS measurements are an attractive NOT tool since the positron causes 
negligible perturbation of the electronic structure as it passes through a 
crystal (Johnson et al_. 1978) and the positron-electron Interaction does not 
appreciably alter the momentum distribution of the electrons (Szotek and 
Chmielewskl 1981). A further advantage of the technique is that it yields 
information regarding both the electron momentum distribution as well as the 
electron density within the sample (Siegel 1980). 

The three most common PAS techniques are lifetime measurements, angu- 
lar correlation, and Doppler broadening. The effective electron density in 
the region from which the positron annihilates can most easily be determined 
using lifetime measurement techniques. These measurements usually use 22 Na as 
a source material since it simultaneously (< 10 psec) emits a y-ray along 
^ ththe J°Jl tr .- n : - Thi s 1 * 28 acts to start the timing apparatus 



which is stopped by the arrival of the ~ 511 KeV y-ray signifying the occur- 
rence of the annihilation event. The positron lifetime depends on the con- 
centration and distribution of the electrons in the region of the material 
through which the thermalized positron drifts before its annihilation. This 
"mean- lifetime will of course be affected by any process, such as the 
trapping of positrons at defects, which alters the local concentration of 
electrons seen by the positron (Coleman 1977). These mean lifetimes, which 
are on the order of hundreds of picoseconds, can be measured to an accuracy of 
a few picoseconds with modern equipment. 

As mentioned previously, when the annihilation event does occur the 
electron-positron pair will usually have a small net momentum in the labora- 
tory frame. This net momentum affects both the energy of the resulting Y -ray 
as well as its direction of emission as shown in Figure 1. The relationships 
between the momentum of the annihilating y-ray, the measured Doppler shift In 
energy, and the measured angle H are described in Appendix B. 

\ 

\ Center of Mass 
\ Velocity V 




Figure 1. A Schematic Diagram Showing the Effects of the Electron-Positron 
Pair Momentum on the Energy and Direction of the Two Resulting 
Annihilation y-Rays in the Laboratory Reference Frame. 
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Angular correlation measurements, which measure the deviation from co- 
linearity of the two photons emitted from the annihilation process, were first 
reported by Debenedetti and his coworkers in I960, The so called long-slit 
geometry (DeBenedetti 1950) is commonly chosen since it provides both good 
angular resolution and adequate counting rates. However, other geometries can 
be employed for specific applications. The angular distribution of a sample 
is determined by measuring the coincidence counting rate as a function of 
angular deviation from colinearity. In well annealed metals this data gener- 
ally takes the form of a central inverted parabola with Gaussian wings. The 
instrument resolution is good, typically ~ 0.5 mrad, but the counting rate is 
comparatively slow and the required source strengths are relatively high. 

Alternatively, Ooppler broadening measurements could be used to deter- 
mine the net momentum distribution of the annihilating electron-positron pairs 
by measuring the energies of the resulting Y-rays. The center of mass momen- 
tum of the annihilating pair along the line of emission of the Y-ray gives 
rise to Ooppler shifts of the Y-ray energies. The distribution of these 
momentum components then results in a Ooppler broadening to give the 511 KeV 
line a width of approximately 511 (V F /2c) KeV where V F (~ 108 cm /s) is the Fermi 
velocity of the conduction electrons and c is the speed of light (Hughes 
1982). Details of the experimental procedure for collecting Doppler-broadened 
spectra are given in Chapter II. 

Although lifetime measurements, angular correlation, and Doppler- 
broadening measurements are by far the most common PAS techniques there are 
two others which should be mentioned. One of the drawbacks of using conven- 
tional sources such as 22 Na or 68 Ge is that only certain depths of the mater- 
ial under investigation are sampled extensively by these positrons. Hughes et 
al ♦._ (]9«2) give the eq uation: 
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« = * j M Q C 2 (68) = ± 255.5 ( fifl ) e V (i. 2 ) 
where 6E is the Doppler Energy Resolution in eV 
and 6B is the Angular Correlation Resolution in mrad. 

The best Doppler systems have a resolution of about 1.2 KeV which 
corresponds to an angular correlation resolution of 5 mrad. Typical angular 
correlation systems have a resolution approximately an order of magnitude 
better than this. For Doppler broadening techniques to effectively compete 
with angular correlation techniques in high precision applications it is 
necessary to improve the effective resolution of Doppler PAS. Schultz et a].. 
(mi) point out that theoretically a 10-fold increase in effective resolving 
power of Doppler measurements can be achieved by deconvolution. 

Deconvolution is the mathematical technique used to separate the 
instrumental broadening effects from the true Doppler broadened positron 
annihilation Y -ray spectrum. Instrumental broadening is a complicated problem 
caused by many different factors. For example, physical properties of the 
detector crystal such as frontal dead layers, geometry of the active region, 
and regions of poor charge collection, can all lead to line broadening. The 
primary contribution to instrumental broadening In a solid state semiconductor 
detector is due to statistical fluctuations in the division of absorbed energy 
between ionization and heating of the crystal lattice, i.e. the Fanno effect, 
(McNelles and Campbell 197b). 

The shape of the instrument resolution function, as well as the mag- 
nitude of the line broadening, is most often determined from the 514 KeV Y -ray 
emitted from the isotope 8& Sr . Tne difficulty lies in the faa ^ ^ 

Doppler-broadening of the 511 KeV Y -ray line in the absence of defects in a 
metallic crystal is of the same order (~ KeV) as the detector resolution, and 



14 

changes Induced by defects are significantly smaller. Fortunately, the defect 
induced changes in the deconvoluted Y-ray energy distribution are much more 
dramatic than in the instrumental ly smeared raw data (McNelles and Campbell 
197b). In addition, removing instrumental effects offers another advantage. 
Such a procedure allows direct comparison of experimental results obtained in 
different laboratories. If results are reported without deconvolutlon they 
must be nondiroensionallzed in some way or else the instrument resolution func- 
tion roust also be reported. In either case the interpretation of results is 
complicated in the absence of a deconvolution procedure. 

In light of the above discussion it was decided that the results pre- 
sented in this work would be for deconvoluted Ooppler spectra. Unfortunately, 
the mathematics involved are somewhat complicated. Because of this, most of 
the existing deconvolution schemes are Iterative in nature and hence computa- 
tionally inefficient. One of the major goals of this work was to develop a 
deconvolution scheme for Doppler broadened positron annihilation measurements 
which is model Independent, does not require iteration, and which yields 
reproducible results. 

In addition to the need for deconvolution, DBPAS also currently suf- 
fers from an Inability to distinguish various defect types. For example, 
vacancies and dislocations tend to alter the annihilation lineshape in quali- 
tatively the same manner. That is, in both cases the spectrum becomes more 
narrow as the defect concentration increases. Most of the lineshape param- 
eters currently used to characterize Doppler broadened spectra are unable to 
distinguish the effects of different defect types because they describe only 
gross changes in lineshape. For example, MacKenzie's S-parameter, which is 
the most commonly used lineshape parameter, essentially measures the width of 
l h _ e . . d J str1 . but1on (MacKenzle et a1. 1970). Hence, for any given spectrum such 
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an analysis yields a single value. This is useful if there is only one degree 
of freedom in the series of experiments. If two variables are allowed to 
change, however, it is impossible to separate their effects. Clearly, 
additional parameters are required. The other parameters selected must be 
uncorrected with the S-parameter. That is, if the second parameter chosen is 
a function of the first parameter only, then no new information is gained, 
Campbell (1977) and Schultz et al_. (1981) have discussed the use of other 
possible parameters to describe lineshapes. Although these parameters may be 
more sensitive to variations in specific defect concentrations they are 
correlated with S, and therefore of little help in differentiating defect 
types . 

In light of the above discussion it can be seen that 1t would be de- 
sirable to find a family of lineshape parameters which is capable of disting- 
uishing various types of metallurgical defects. For example, one of the 
groups of parameters to be examined in this work 1s the central moments of the 
Y-ray energy distribution. It 1s well known that the values of these moments 
give information regarding the shape of the distribution. In fact, it can be 
shown that if all of the moments are known one can construct the distribution. 
The even moments give information about the width of the distribution while 
the odd moments vanish for symmetrical distributions such as the ones to be 
examined in this work (the symmetry question 1s addressed in Chapter III). 
The second central moment is called the variance and its square root 1s the 
familiar standard deviation. The fourth central moment is called the Kurtosis 
or flatness. The advantage of moment analysis is that the information ob- 
tained is not redundant. For example, it is very possible to have two dis- 
tributions with equal standard deviations (second moments) and yet widely 
d i f f ere iit f qu irt h moment s . Such a fam ily of lineshape p ara meters suggests the 
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possibility of separating the effects of multiple metallurgical defect types 
using OBPAS. 

The above brief review of positron annihilation spectroscopy as an NO! 
tool for the study of metallurgical defects indicates that the technique with 
the greatest potential for practical use outside the laboratory 1s Uoppler- 
broadenlng. This technique offers the advantages of low source strength re- 
quirements, the ability to use 68 Ge as a source material 1n order to reduce 
surface effects, and few geometrical experimental constraints. Ooppler broad- 
ening currently suffers from relatively poor energy resolution and an Inabil- 
ity to distinguish different metallurgical defect types. Based on these ob- 
servations, the purpose of this work is twofold: 1) to develop an efficient, 
model independent, noniterative deconvolution method, and 2) to develop a 
family of llneshape parameters capable of distinguishing various defect 
types. 



CHAPTER II 
EXPERIMENTAL PROCEDURES 



This chapter is divided into two sections. The first section des- 
cribes the Ooppler broadened PAS instrumentation and the methodology employed 
to collect Doppler spectra. The second section outlines the thermomechanical 
treatments used to create the desired defect concentrations in the metallic 
samples which were studied in this work. 

The Popper PAS System 

The heart of the Doppler broadened PAS instrumentation system is the 
single crystal high purity Germanium (HPGe) detector. This detector is mount- 
ed at the end of a cold finger in a cryostat and, to reduce thermal noise, is 
operated at liquid nitrogen temperatures. The cylindrical detector crystal 
used in this work is 46.7 m in diameter and 46.2 m 1n length. When the 
annihilation y-rays impinge upon the crystal their energy is absorbed by two 
mechanisms. The first is the creation of electron-hole pairs within the 
crystal. These electron-hole pairs are subsequently swept apart by the 
applied detector bias (4000 V) and flow through an external circuit creating a 
current pulse whose integral with respect to time has an amplitude propor- 
tional to the incident r-ray energy. The second absorption mechanism is heat- 
ing of the crystal. The ratio of these two absorption mechanisms is known as 



the Fanno factor and the statistical fluctuations in the division of absorbed 
eneruy between ionization and heating Is the primary contribution to instru- 
mental line broadening In solid state semiconductor detectors. For further 
details of the absorption process see Moore (1983). 

The resulting current pulse from the detector 1s rather weak and un- 
suitable for cable transmission. Because of this a preamplifier (Ortec Model 
120-b) is mounted directly adjacent to the detector output. The detector and 
preamp are O.C. coupled, as is the entire system, and powered from the same 
A. C. source to prevent ground loops. The preamp is a charge sensitive device 
capable of accepting signals corresponding to the Y-ray energy range of 0-10 
MeV. Its output is a voltage pulse with an amplitude proportional to the 
integrated charge supplied by the detector and hence proportional to the inci- 
dent y-ray energy. The magnitude of the voltage pulse amplitude 1s related to 
the incident Y-ray energy by the approximate conversion factor of 280 mV per 
MeV. Therefore for a typical 511 KeV Y-ray the output of the preamp would be 
approximately 143 mV. In addition to converting the detector charge pulse to 
a voltage pulse with an intensity suitable for cable transmission the preamp 
also shapes the resulting voltage pulse so that it is acceptable to the linear 
amplifier. 

The linear amplifier used in this work is an Ortec Model 571. It 
accepts the shaped pulses from the preamp and amplifies them to suitable volt- 
ages for transmission to the multichannel analyzer. The output pulses are 
active filter shaped (skewed Gaussian approximations) and have amplitudes in 
the 0-10 V range. Unipolar outputs are used to optimize the signal-to-noise 
ratio (resolution). Since the main amplifier is operated within its linear 
range the amplitude of the output voltage pulse is proportional to the inci- 
dent Y-ray energy as measured by the HPGe detector. 
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The output from the linear amplifier Is routed to an Ortec Model 408A 
biased amplifier. This device accepts the output pulses from the main ampli- 
fier and allows one to selectively view a:.y region of the initial 0-10 MeV en- 
ergy range. In this work the biased amplifier was used to select the energy 
window of 483 to 548 KeV. The output pulses from the biased amplifier serve 
as the Input to the analog-to-digital (ADC) converter located within the Nor- 
land IT-5300 multi-channel analyzer (MCA). In this device the amplitude of 
the voltage pulse Is converted to a digital representation of the energy of 
the detected T-ray and subsequently dynamically displayed on the MCA's CRT in 
the form of a histogram of energy versus count intensity. The MCA can be pro- 
grammed to accept input signals until a predetermined number of counts in the 
peak channel is reached. At this time the experiment is terminated and a com- 
plete Uoppler broadened spectrum 1s displayed on the CRT. The energy distri- 
bution is then ready for transmission to the computer for subsequent data man- 
ipulation. 

The two most important parameters which characterize the quality of a 
Doppler PAS system are the width, in units of energy, of the instrument reso- 
lution function and the energy resolution, in units of eV per channel, of the 
MCA. The width of the instrument resolution function (IRF) is, in essence, a 
measure of the quality of the detector crystal, although It Is also affected 
by the associated amplification electronics as well as the ADC. The IRF is 
most often determined by examining the response of the system to the monoener- 
getic (E * 514 KeV) photons emitted from the Isotope 85 Sr. If the instrumen- 
tation was perfect the energy distribution, as recorded by the MCA, would con- 
sist of a single bar in the histogram in the channel corresponding to 514 KeV. 
However, in any real instrumentation this line will be broadened for reasons 
P rt»Jjou»1jr d1scu ssed jL_JM. am P, u J?JL Pi 1 nst rum g ntal. bim d^aiag.is. .m os 1. of t eo-re.-- 
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ported in terras of the energy width of the 85 Sr line as measured across the 
two channels which contain half the number of counts in the peak channel. 
This energy measure is commonly named the full width half maximum (FWHM). 
Figure 2 shows typical IRF for the Doppler PAS system used in this work. As 
can be seen in this figure the FHHM of the IRF Is about 1.28 KeV. 

The second parameter which describes the quality of the Doppler PAS 
instrumentation is the energy resolution of the MCA which is determined by the 
following procedure. First, a Doppler broadened spectrum is collected and the 
channel corresponding to its centroid is determined. The centroid channel is 
assigned the energy value 511 KeV. Next, an energy spectrum Is obtained for 
the isotope 85 Sr. The channel corresponding to the centroid of this distribu- 
tion is calculated and assigned the energy value 514 KeV. The MCA energy 
resolution is then easily calculated by dividing the energy difference (3 KeV) 
by the number of channels between the centroids. In this work the 511 KeV 
channel was number 439 while channel number 486 corresponded to 514 KeV. 
Hence, the MCA energy resolution is 3000 eV/(486 -439) or " 64 eV per channel. 
Clearly this value can be adjusted using the biased amplifier. The reason for 
selecting the 64 eV/channel is to obtain the best possible resolution while 
maintaining the ability to display the entire region of interest (483 to 548 
KeV) on the 1024 channel MCA. The reasons for requiring the above energy 
window will be discussed in Chapter III. 

While the MCA energy resolution is a parameter which, to a certain ex- 
tent, can be varied, the IRF is essentially predetermined by the instrumenta- 
tion. The one factor contributing to the width of the IRF which can be con- 
trolled is the data acquisition, or count rate. It is well known that in- 
creasing the count rate decreases the instrument resolution. Therefore, one 
must select a count raje .which results. J n^M^cceptab Inline broadening, wJhile 
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keeping the time required to collect a spectrum reasonably short. Moore 
(1983) demonstrated that with the Ooppler instrumentation used In this work 
very reproducible results can be obtained with count rates of < 2 kHz. The 
count rate Is easily adjusted by changing the sample to detector distance. 
All of the results reported in this work were obtained using count rates of 
< 2 kHz with an associated experiment time of - 30 minutes resulting in a peak 
height of 40K counts and " 3 x 10 6 total integrated counts. 

The positron source isotope in all experiments was 68 Ge which was 
chosen to avoid surface effects. The mean depth of penetration in aluminum 
for such positrons is ~ 0.50 mm. The activity of the source was initially 50 
uC. The half-life of 68 Ge is - 275 days and hence the activity of the source 
did not change significantly for any single set of experiments since a typical 
set of measurements requires less than a week to accumulate. It should be 
pointed out that the activity of the source does not measurably affect the 
observed annihilation spectra, provided the count rate remains constant. 

Although the count rate can be easily adjusted by changing the source 
to sample distance such a change introduces an additional error source. This 
error source is related to the source/sample geometry. It has been demon- 
strated that the observed spectra may depend on the angle with which the gamma 
rays impinge upon the detector crystal (Sroedskjaer 1983). It has also been 
shown that in high precision experiments among samples of the same material, 
one should use samples of nearly identical thickness to reduce the error asso- 
ciated with small angle scattering (Smedskjaer 1983). 

The 68 Ge source was deposited on a thin (~ 0.05 nm) nickel foil for 
handling convenience. This results In what 1s known as a source contribution 
to the annihilation lineshape. That is, in each experiment some of the posl- 
- trons . w1 H. ann .^ 1 l at . e 1n the n1 cke1 fon and the v-ray s produ c ed by,, th ese 
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events will contribute to the overall annihilation spectrum. However, the 
ratio of source annihilation events to sample annihilation events is small (< 
1%), and the source contribution 1s neglected in the remainder of this work. 

In all measurements the source and sample were arranged in what is 
commonly called the sandwich configuration. In this configuration the posi- 
tron source is placed between two identically treated metallic samples. The 
samples are thick enough, typically a few millimeters, to guarantee that 
virtually all of the positrons annihilate within the material being examined. 



Specimen Preparation 

The material examined in this work is high purity, large grain size 
aluminum purchased from the Materials Research Corporation. The aluminum was 
zone refined to a nominal purity of 99.9995* before being vacuum melted and 
cold fabricated to size. The chemical analysis, in PPM of metallic inter- 
stitlals, is listed in Appendix C. The aluminum was supplied as a 13 mm 
diameter rod from which - 1.5 mm thick disks were machined. Each of these 
disks had an ~ 1.5 mm hole drilled through it near its edge. The purpose of 
this hole was to allow the disks to be suspended from a titanium wire while 
within the vertical tube furnace. After the disks were machined to size, they 
were given a full anneal heat treatment consisting of 2 hours at 400°C fol- 
lowed by furnace cooling overnight. The as annealed aluminum disks were the 
startiny material for all subsequent thermal treatments. 

The thermal treatments used in this work are very similar to those 
outlined by Cotterill et aK (1972). In this work the authors describe a two 
step time/temperature cycle which allows one to produce, in a controlled 
f ash1 on , ^rjrtjampj^s -typj^_.and_t,tieii., 
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samples with dislocation loops as the major defect type. The actual thermal 
treatment employed in this work is as follows: 

1) Heat the samples to the desired quenching temperature and hold 
them at temperature for one hour to establish thermal equilibrium. 

2) Quench the samples from the furnace into ~ 293°K water. 

3) Dry the samples within 5 to li seconds. 

4) Quench from ~ 293°K to a liquid nitrogen bath at 77°K. This 
procedure has been shown to produce primarily vacancies, the concentration of 
which is a direct function of the quenching temperature. 

b) After DBPAS characterization, the aluminum disks are then annealed 
at 3S3°K for thirty minutes, which results in the production of prismatic 
dislocation loops by vacancy coalescence. 

6) The samples are cooled to 77°K and the PAS measurements corres- 
ponding to the dislocation loops are made. 

7) Finally the disks are annealed at b23°K for two hours to remove 
the dislocations and return the material to the annealed state. 

Several of the steps in this procedure require further explanation. 
The quenching procedure is taken from Cotterill et aK (1972) and their 
results together with the following discussion forms the basis for the 
assumption used in this work that the PAS measurements made after step 4 are 
for a fairly uniform distribution of vacancies of concentration approximately 
equal to the equilibrium concentration at the quenching temperature. 

Quenching experiments are usually performed on thin metallic foils so 
that the quench rate remains constant across the sample. However, as dis- 
cussed previously, constraints imposed by the positron annihilation experi- 
ments put a lower limit on the sample thickness of ~ 1.5 mm. For a complete 
review of th e complications which can eyolye.J^^ 
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sample thickness see, for example, Cotterlll et art. (1965). The main problems 
to be addressed here are the formation of vacancy clusters and the loss of 
vacancies at the sample surfaces. The second problem was investigated in a 
paper by Mikes and Harrison (1972) in which they observed the quenching of 3 
mm thick disks of aluminum and reported a very low vacancy loss. Since the 
samples used in this work are only ~ 1.5 mm thick and since the quench pro- 
cedure is similar to that employed by Wilkes and Harrison, it will be assumed 
that the vacancy loss during the quench can be ignored. As for the formation 
of vacancy clusters, the positron lifetime measurements of Cotterill et al. 
(1972) demonstrate that for this quenching procedure the probability of form- 
iny bivacancies, or larger clusters, is small. 

It is very important to dry the samples thoroughly between the water 
quench and the liquid nitrogen quench. This is necessary to avoid the forma- 
tion of a thin layer of ice on the specimen surface which would contribute an 
undesirable component to the measured PAS spectrum. As pointed out by 
Cotterlll et aK (1972) the annealing of the vacancies during this short dry- 
ing procedure may be neglected since the mean lifetime of vacancies (at room 
temperature) in a specimen quenched from 873°K is ~ 100s and increases for 
lower quenching temperatures. It is, of course, desirable to keep the drying 
time as short as possible. 

Step 5 of the procedure, the formation of prismatic dislocation loops 
from the collapse of disk shaped vacancy clusters during a low temperature 
anneal of aluminum, was reported by Hirsch et aU in 1958. Cotterill et al. 
(1972) demonstrate that the appropriate thermal treatment consisted of a 30 
minute anneal at 353°K. The theoretical relationship between the vacancy 
concentration after the quench, C v , and the concentration of dislocation 
c ores, C d (gi ven as an a tomic fraction), is presented in Appendix G. 
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The apparatus used to perform the thermal treatments is shown m 
Figure 3. A quartz tube passes through a vertical tube furnace and is 
attached to the room temperature water reservoir by a "quick disconnect" 
coupling. Just above the water level in the reservoir is the argon supply 
port. All thermal treatments were performed in an argon atmosphere to avoid 
oxidation. The other end of the tube is connected to a flow meter which is 
used to monitor the argon flow rate. The other side of the flow meter is 
vented to atmosphere. The temperature is measured by a chrome! -a lumel thermo- 
couple positioned inside the tube furnace and adjacent to the aluminum disks. 
These disks are supported by a titanium wire, the other end of which is con- 
nected to an iron slug. The slug, and hence the specimens, are held in posi- 
tion in the furnace by a magnet located outside of the tube. The actual 
quench 1s performed simply by pulling the magnet away from the tube and allow- 
ing the entire specimen-wire-slug arrangement to fall by gravity into the 
quench bath. It is, of course, important for the system to be designed so 
that the specimens fall freely through the tube and quickly reach the quench 
bath. If the specimen-wire-slug arrangement does not fall freely the tempera- 
ture of the samples when they enter the water will be significantly lower than 
the temperature measured by the thermocouple. 
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CHAPTER III 
DECONVOLUTION OF DOPPLER PAS MEASUREMENTS 



As discussed In the introduction, compared to angular correlation (AC) 
measurements, the major disadvantage of the Doppler PAS technique is Its in- 
herently poor energy resolution. In order for floppier methods to be competi- 
tive with AC methods in high resolution applications the instrumental line 
broadening must be removed from the experimentally obtained Doppler spectra. 
Mathematically, the technique which accomplishes this separation is deconvolu- 
tion. In this chapter the deconvolution of floppier PAS measurements will be 
discussed. 

The chapter Is divided into four sections. In the first section the 
mathematics of deconvolution are presented. This section justifies the use of 
deconvolution techniques to separate instrumental effects from the desired 
true floppier-broadened energy spectrum. In addition, the inherent advantages 
and mathematical difficulties associated with the use of Fourier analysis to 
solve the problem are addressed. 

In the second section of this chapter the methods Introduced by other 
investigators to solve the incorrectly posed deconvolution problem, along with 
their inherent advantages and disadvantages, are discussed. The third section 

presents the FFT/power spectrum deconvolution method (Schaffer et al. 1984); 

and describes its use in detail. In the final section of this chapter the 

FFT/power spectrum method is demonstrated using a neutron irradiated high 

purity single crystal sapphire sample. 
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The Mathematics of Deconvolution 

The first objective of this section is to demonstrate the validity of 
usiny the convolution equation to describe the relationship between h(x), the 
experimentally observed energy distribution, f(x) t the true Doppler broadened 
energy distribution, and r(x), the instrument resolution function (IRF). In 
Figure 4(a) the instrument resolution function r(x) is shown. Recall that this 
function represents the instrumental broadening of what should be a single 
energy line at the centroid of r(x). The function r(x) is assumed to be norm- 
alized to unit area. In order to understand the relationship between these 
three functions, f(x) should be viewed as a histogram comprised of segments of 
width Ax. If it can be assumed that the instrumental response is linear (the 
detector and all associated electronics are operated within their linear range) 
then each segment of the f(x) histogram will be broadened by the unit area IRF 
r(x). Hence, the value of the observed function h(x) at x' will have contribu- 
tions from many of the f(x) histogram segments in the neighborhood of x 1 . Fig- 
ure 4(b) is a schematic representation of the contribution to h(x*) from the 
histogram section of f(x) at a distance X) from the origin. This contribution 
Ah^x') will be the product of the area f (xj )Ax and rfx'-xj), that is: 

Ahi(x') - fU^rU' - X])Ax (3.1) 
The total value of h(x') will be the sum of the contributions from all pos- 
sible histogram segments in f(x). Thus, 

00 CD 

h(x') = I Ah^x 1 ) = I f(x 1 )r(x l - Xi )A X (3.2) 
i =-« 

or as the width of the segments approaches zero the summation in Eq. 3.2 is 
replaced by an integral to yield, 

oo 

h(x') = / f(x)r(x' - x)dx (3.3) 
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Clearly the same procedure holds for any value of x In h(x) and the final form 
of the relationship, with a change of variables, 1s: 

h(x) 8 L f M r < x - y)<iy (3.4) 

Eq. 3.4 is the convolution equation which serves as the starting point for the 
remainder of the mathematics used In this work. 

The advantage of using Fourier analysis when dealing with convolution 
problems can be seen clearly when Eq. 3.4 is rewritten as: 

H(w) = F(w).R(w) (3<5) 
where the functions In Eq. 3.5 are the Fourier transforms of the respective 
functions in Eq. 3.4 according to the general relationship: 

G(W) Vii L e ' 1WX9(x)dX (3.6) 
Note that in Fourier space convolution is simply a complex multiplication. 
The equivalence of Eqs. 3.4 and 3.5 is demonstrated in Appendix 0. 
Rearranging Eq. 3.5 yields 

F(w) = H(w)/R(w) ( 3 4?) 
from which one can find f(x) using the Inverse Fourier transform 

Eq. 3.8 is an expression for the true Ooppler broadened energy spectrum In 
terms of the Fourier transforms of the experimentally observed function h(x) 
and the Instrument resolution function r(x). However, the uniqueness of the 
solution given in Eq. 3.8 requires that h(x) be known exactly. This condition 
is never satisfied in a real experiment and this fact significantly compli- 
cates the mathematics of the problem. The observed function h*( x ) will be 

related to ^ j ™^^ for an u nknoMn 
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function f(x) by: 

max |h*(x)-h(x)|<e (3.9) 
-*><y<+» 

where e is a measure of the noise associated with the stochastic nature of the 
experiment. 

It can then be shown (see Appendix E) that there are two solutions 
f](x) and f2(x) = f](x) + C sin(2imx/(b-a)) with corresponding function hj and 
h£ (as defined by Eq. 3.4) both of which satisfy the inequality (3.9). Hence, 
both f] and f£ must be considered as solutions to Eq. 3.4 with a left hand 
side h*(x) taken from an experiment to within an error e (Turchin et al . 
1971). In fact, there is an infinite set of such solutions a subset of which 
contain indefinitely large rapidly oscillating components. 

In order to understand the mathematical source of these difficulties 
Eq. 3.8 was investigated further. If in Eq. 3.8 H(w) is replaced by H*(w) 

with °° 

H*( W )=J-/ e iwx h*(x)dx (3.10) 
¥ r 2T -« 

and h*(x) = h(x) + 6( x ) (3.11) 

where 6 is a stationary random process with zero mean, then it has been shown 
(Kautian 19S8) that the dispersion of f(x) is 

CO 

0(f) « <f 2 > - <f> 2 = / G ( w) 9 dw (3.12) 

fit |R(w)| 

where G(w) is the power spectrum of 5(x) 

and |R(w)| + 0 for |w| ♦ • 
Therefore the dispersion of the solution f(x) can remain finite only if the 
power spectrum of the noise approaches zero rapidly as w approaches infinity. 
In practice this will not be the case since the noise will inevitably contain a 
"white noise" component such that G(w) remains finite as |w| + * and the dis- 
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persion of the solution f(x) therefore becomes infinite (Turchin et al_. 1971). 

Alternatively one can understand the source of the mathematical diffi- 
culties as follows. The observed function h*(x) contains high frequency noise 
components which are not physically realizable in the true function h(x). In 
Eq. 3.8 these high frequency components are divided by small eigenvalues, R(w), 
which results in larye oscillations in the solution f(x) (Turchin et aK 1971). 

Recognizing these difficulties, previous investigators have proposed 
several deconvolution techniques which appear to successfully solve the prob- 
lem. Since the accuracy of the deconvoluted spectrum can never really be de- 
termined, except in certain carefully constructed numerical experiments, the 
criticisms of a specific method tend to relate to the amount of computer time 
required, the complexity of the programming effort, the degree to which the 
analysis may be automated, or to the need for some a priori information con- 
cerning the function f(x). As a preface to the introduction of the FFT/power 
spectrum deconvolution technique, the following section is a review of the 
prior techniques that have been developed to deconvolute Doppler PAS data. 



Previous Deconvolution Techniques 

One of the early Doppler PAS deconvolution schemes introduced by Hotz 
et al_. (1968) is basically a matrix approach. The convolution in its discrete 
form is 

a 

n(x) « I r(y - x)f(y)Ay ( 3 j 3) 
y a -a ' 

where Ay is the incremental width (resolution) 

and the other functions have their usual meaninys 

This can be rewritten in matrix form as 
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b i s A j1 c j (3.14) 
where b is the discrete experimental vector (spectrum) 

c is the unknown discrete vector before instrument broadening 
and A is a rectangular matrix in which each row is the detector re- 
sponse function appropriate to the corresponding component of c. 
In the method developed by Hotz et al_. (1968) the matrix A is determined using 
a5 Sr and the vector b is the output of the Doppler system MCA for a given 
experiment. In this technique an Initial estimate of the unknown vector c is 
made, say c(D, and a corresponding vector b(0 is calculated from Eq. 3.14. 
The quality of cO) is determined by finding X 2 for the fit between b and bO) 
as defined by 

X 2 * J (b (l) . bi)/bi (3-16) 

Next cO) is modified in a manner described in the original work to create a 
second estimate, c(*), of the desired function c. Then b( 2 ) 1 S found from Eq. 
3.14 and a new X 2 value determined. The process is repeated until X 2 is mini- 
mized. This process is extremely involved and inefficient. In addition the 
initial estimate cO) as well as the modification criteria used to generate 
subsequent c( 1 ) estimates can influence the final result. 

In more recent years the matrix approach has been examined in greater 
detail (Aulick and Gallie 1981). Although this work is not related to Doppler 
PAS measurements, the mathematical problem is the same. It Is demonstrated 
that 3.14 should be rewritten as: 

b = Ac (3 #16 ) 

where b = b + e 

and e is the measurement error 
The least s quares solut i on of E q« 3.V6is then 



A*b » c 
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(3.17) 



where A* is the pseudo-inverse (Penrose 1956) of A. 
However, due to the fact that b contains noise Eq. 3.16 is incorrectly posed 
and A will be ill-conditioned. Hence, the elements of At are large and the 
error in b will be magnified, obscuring the true solution which is 

c = Atb (3.18) 
Aulick and Gallie (1981) point out that to avoid this problem one needs to 
calculate solutions to Eq. 3.16 in a way which will damp out the error in b. 
One way of accomplishing this goal is to supply some additional a priori in- 
formation about the shape of the desired vector c. The additional Information 
recasts the problem as a constrained optimization problem. This technique is 
called regularization. 

An example of the type of a priori information commonly supplied 
appears in the deconvolutlon technique used by Reddy and Carrigan (1970). In 
their approach the resolution function is assumed to be the sum of two 
Gaussians with the same centroid but differing in their amplitudes and widths. 
The shape of f(x) (returning to the more common notation) is assumed to be a 
truncated parabola based on the free electron model (Reddy and Carrigan 1970). 

Introducing these constraints on the shapes of r(x) and f(x) into the 
least square solution problem yields a physically plausible result. That is, 
the function f(x) which is obtained by this method does not contain indefi- 
nitely large and rapidly oscillating components typical of the results obtain- 
ed by evaluating Eq. 3.8 directly. The difficulties associated with this 
method are related to the validity of the a priori assumptions. For example, 
the resolution function r( x ) is constrained to be symmetric. In addition as- 
suming the shape of the true Doppler energy distribution, f(x), to be a simple 
..— . * r H"cated .parabo1.a,,is„p.cobab1y .an ■■O^lnpimcatJy OthA^n^c^ 
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(MacKenzle et al_. 1976) suggest the modeling of f(x) using the superposition 
of Gaussian and inverted parabolic functions as representing contributions 
from core and valence electrons respectively. While it 1s not the purpose of 
this work to discuss the validity of the various electron distribution models, 
the very fact that there exists more than one analytic prediction of the shape 
of f(x) Is strong reason to develop a deconvolution technique which places 
minimal constraints on the shape of f(x). 

One of the most widely used deconvolution techniques in Doppler PAS 
studies is an iteratively based method introduced by Jackman et ajL (1974). 
In this method an initial estimate of the desired result f(x), say fi(x), is 
selected. The estimate f](x) is then convoluted with r(x) to obtain hjfx), 
the first estimate of the experimentally obtained energy distribution h(x). 
The function hi(x) is then compared point-by-point with h(x) and if satisfac- 
tory agreement 1s obtained fj(x) is assumed to be a reasonable approximation 
to the true solution. If unsatisfactory agreement is obtained then f](x) is 
modified to produce f2(x) and the process is repeated. 

The initial estimate of f(x), the modification procedure for obtaining 
fg(x) from f](x), and the criteria for stopping the iteration procedure are 
all current topics of investigation 1n the literature. In the original method 
the experimentally obtained function h(x) served as the initial estimate 
f|(x). The modification algorithm was 

i+2 i+2 
= f l(xi) x ( I h(xj)/ I h]^)) (3.19) 
j=i-2 j=1-2 J 

Subsequent work demonstrated that the convergence rate of the Iteration pro- 
cess could be improved if the summations in Eq. 3.19 were taken over nine 
channels instead of the original five and if the summation ratio was squared 
t)efore . .'"".I VlP.lj ^?.PJl..^-- f .L( x i.) ,(^ njLef ^er_a nd Kerr 19 75J^_ _Ihe iteraH nn p r _o= 
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cedure is stopped when the comparison of h(x) and h N (x) yields a x2 value less 
than some predetermined minimum (Jackman et al.. 1974) or when the sum of the 
squares of the deviations falls below some predetermined level (Dannefaer and 
Kerr 197b). 

The main advantage of these iteratively based techniques is that a so- 
lution is obtained without any restriction on the shape of f(x). While this 
is indeed an attractive feature of the method it also suffers from some inher- 
ent disadvantages. These disadvantages include the amount of computer time 
required, the possible dependence of the result on the initial estimate of 
f(x), and the possibility of introducing large oscillations into the solution 
if the stopping criterion is selected with too small an error bound. This 
last difficulty Is of course related to the mathematical nature of the incor- 
rectly posed problem. 

The iterative deconvolution technique essentially Ignores the statis- 
tical nature of the problem. An alternate approach is to treat Eq. 3.4 as a 
problem of mathematical statistics and introduce a regularizer (a priori In- 
formation) in the form of some probability distribution for f(x) (Turchin et 
ai- 1971). This 1s the method of statistical regularizatlon which offers the 
advantage of error analysis. The result of this process is the most probable 
f(x) for a given h(x) measured with mean square errors S(x) which are indepen- 
dent of x and normally distributed with zero mean (Goldanskii et al_. 1978). 
The regularizer is a form of a priori information concerning the shape of f(x) 
and hence somewhat restrictive. However, such an approach appears to yield 
reasonable results if a solution can be obtained using minimal a priori infor- 
mation concerning the smoothness of f(x). 

In 1981 the Fourier cut-integral method was introduced by Szotek and 

.Chmi el ewski..-. In -thi s— techniaup . f n.„_a_A 
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««>-=k/V te) 4$* (3 - 20) 

The difference in the equations is the Introduction of finite limits of inte- 
gration. If p 0 is then chosen small enough the solutions obtained are physi- 
cally plausible. The mathematical explanation for this result is that the in- 
sertion of finite limits of Integration forces the dispersion of f(x) to 
remain finite. (Refer to Eq. 3.12 and the associated discussion). The prob- 
lem associated with this method is that the optimum value of p 0 is not readily 
predictable. Approximations for f(x) are obtained using different values for 
p 0 , and one must then select the "best" approximation to the true f(x). This 
is typically done by convoluting the function with r(x) and then comparing the 
result with h(x). As p 0 becomes too large the solution begins to exhibit 
large oscillations as would be expected from the mathematical theory. 

In 1982 Verkerk published an extensive mathematical treatment of the 
deconvolutlon problem. In this work he discusses the relative strengths and 
weaknesses of many of the currently employed deconvolution techniques. Ver- 
kerk then proposes a method which in many ways 1s very similar to the Fourier 
cut-integral method. The crux of this method is that the high frequency 
Fourier coefficients of h(x) do not contain significant information and conse- 
quently only a limited number (say K) of Fourier coefficients H(w) contain the 
essential information from measurement of h(x). It should be pointed out that 
if H(w) contains N Fourier coefficients then deleting the N-K high frequency 
coefficients in the reconstruction of f(x) is equivalent to inserting finite 
limits of integration in Eq. 3.4 as in the Fourier cut-integral method. 

The deletion of the high frequency coefficients is equivalent to a 
rectangular low-pass filter and it is well known that the application of such 
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a filter to F(w) will cause nonphyslcal oscillations in the approximation to 
f(x). 

In summary, there are two different types of errors associated with 
the estimate of f(x): 1) statistical errors, which are amplified by the de- 
convolution operation and 2) systematic errors due to filtering, which are 
largest in the peak region and may cause nonphyslcal oscillations 1n the tails 
of the deconvoluted signal. Error type 1 Increases as the number of Fourier 
coefficients used in reconstruction increases, while error type 2 decreases 
(Verkerk, 1982). Thus, by using such an approach the problem is reduced to 
selecting the optimum number of coefficients. One disadvantage of Verkerk 's 
method, as with the Fourier cut-integral method, 1s that the procedure must be 
repeated several times for different numbers of coefficients until the optimum 
number is determined. It would be desirable to develop a simple one-step 
method to determine the appropriate number of coefficients to use in recon- 
structing f(x) from F(w), and indeed this Is one of the primary advantages of 
the FFT/power spectrum deconvolution method which is described in the follow- 
ing section. 



The FFT/Power Spectrum Deconvolution Method 

The deconvolution method developed 1n this work is a fast Fourier 
transform (FFT) based method which is simple to execute, fast, easily auto- 
mated, and assumes little or no a priori information concerning the shape of 
f(x). The distinctive feature of this method Is that the number of Fourier 
coefficients retained in the reconstruction of f(x) Is determined by examining 
the power spectrum of the experimental data. Using the observation that the 
power spectrum falls steeply to a minimum before rising due to the white noise 
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coefficients at this first minimum. The logic behind this approach is rather 
simple. Retaining fewer coefficients than this implies a smoothing of the 
deconvolved spectrum and the loss of potentially valuable structural informa- 
tion. Adding more coefficients gains nothing because the signal-to-noise ratio 
is of order one at the minimum. In addition, truncating the coefficients at 
the point where the white noise begins to dominate forces the dispersion of 
f(x) to remain finite. 

With this brief introduction it is now possible to list the six steps 
which together constitute the FFT/power spectrum deconvolution method: 

1) Both the raw data h(x) and the resolution function r(x) are 
background corrected. 

2) The functions h(x) and r(x) are Fourier transformed using an 
appropriate FFT algorithm. 

3) The Fourier coefficients of h(x) are examined and those up to and 
including the first minimum in the power spectrum of h are 
retained, the rest are replaced by zeros. 

4) The Fourier transform of f(x) is obtained by the point by point 
division of the Fourier transform of h(x) by the transform of r(x). 

5) f(x) is obtained by the inverse FFT. 

6) The solution is checked by forming the convolution of f(x) with 
r(x) and comparing the result with h(x). 

Each of these basic steps will now be described in further detail. 
Semiconductor detectors usually produce asymmetric Y-ray peaks with a non- 
linear underlying background. This is caused primarily by three processes: 
1) incomplete charge collection, 2) electron escape from the active volume, 
and 3) scattering of Y-rays (Chaglar et aK 1981). This background component 
^}^.^L^ms^Ll£fmJit>B exp e rimenta ] 1 v^bta4 ned-j^nw^d^ — — _ ; __-f~= 
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the deconvolution procedure 1s performed. The background correction routine 
employed in this work Involves a least squares fitting of a linear function in 
both the low and high energy tail regions. The region on the low energy side 
of the peak spans from 486 to 497 KeV, while the high energy region extends 
from 534 to 542 KeV. The straight line segments fitted in these regions are 
extrapolated to 5U3 and 518 KeV respectively where they are joined by a third 
line segment to produce a simple background approximation (Jackman et al, 
1974), The linear least squares analysis algorithm used in this work is a 
simplification of a higher order least squares routine by Gerald (1978). As 
in the work of Jackman et aK (1974), use is made of only the high energy half 
of the background corrected spectrum since it is known that incomplete charge 
collection of solid state detectors results 1n low energy tailing. In theory 
the spectrum should be symmetric about 511 KeV since statistically one should 
observe equal numbers of red and blue Ooppler shifted y-rays. 

Before proceeding to the details of the FFT algorithm one important 
point requires explanation. In order to correlate MCA channels with a cor- 
responding energy range, one must know the MCA energy resolution per channel 
as well as which channel corresponds to the center of the energy distribution 
and hence 51 1 KeV. The former value is determined as described in the experi- 
mental procedures. However, the center of the distribution must be calculated 
for each experimental spectrum. A simple centroid calculation is unacceptable 
since the higher background on the low energy side would bias the result. 
Similarly, simply selecting the channel with the maximum number of counts is 
undesirable since, due tn the statistical nature of the experiment, the 
centroid channel and the peak channel are not necessarily the same. The tech- 
nique used in this work is to determine the eleven channel region which con- 
tain s^ the lar ge st number _o jL_cgunts and call t he m iddle channel of t his re gion 
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the center of the distribution. In most cases, the channel corresponding to 
the center as calculated above was also the channel with the maximum number of 
counts. 

Proper arrangement of the data string is required before proceeding to 
the FFT algorithm. Since this algorithm requires periodic data, the raw data 
are rearranged in the following way. The first datum point of the rearranged 
set corresponds to the center of the raw data string. The remaining rear- 
ranged data points numbered from 2 to ((N/2) + l) correspond to the remaining 
values of the background corrected high energy half of the raw spectrum. It 
should be pointed out that M must satisfy the relationship N = 2™* where m is 
an integer. Since the real data string rarely has exactly ((N/2) + l) data 
values it is usually necessary to insert a string of zeros at the end of the 
data set to fulfill this requirement. In order to achieve periodicity, data 
point N should correspond to data point 2. In general the points from ((N/2) 
+ 2) to N are defined in the following way: 

g((N/2) + P) = g((N/2) + 2 - P) (3.21) 
where P varies from 2 to (N/2). In the experiment performed in this work m 
was equal to eleven which corresponds to a total of 2048 data points (includ- 
ing the zerofill). 

Once the functions h(x) and r(x) have been background corrected and 
properly arranged they are passed to the FFT algorithms. The FFT algorithms 
employed in this work are those developed by Cooley, Lewis, and Welch (1979). 
There are a number of other such algorithms available which should work 
equally well. In every case it is recommended that the Fourier coefficients 
be calculated with double precision accuracy. The advantage of using FFT's in 
place of discrete Fourier transforms (OFT) can be seen from the following 
example, For..20A8.jdata^ - 
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Implications while the FFT requires 2N Log 2 N or ~ 4.5 x 10* complex multipli- 
cations. In this example the number of calculations required is reduced by 
about two orders of magnitude. 

The justification for using only the Fourier coefficients up to and 
including the one corresponding to the first minimum in the power spectrum of 
h(x) has been briefly discussed. Similar techniques are commonly employed 1n 
electronic signal processing and it can be shown that the number of coeffi- 
cients retained is a measure of the frequency resolution. In the present case 
the resolution refers to the width of the lineshape components which can be 
resolved. The procedure by which the computer calculates the number of re- 
tained Fourier coefficients from H(w) is as follows. Values associated with 
the power spectrum are obtained by multiplying each Fourier coefficient H(K) 
by its complex conjugate. However, since the function h(x) is even, the 
Fourier coefficients H(w) are all real. Truncating H(w) at the K correspond- 
ing to the first minimum in the power spectrum is equivalent to the require- 
ment that the set of retained coefficients be monotonically decreasing, which 
is easily programmed by testing 

H(w + 1) < H(w) (3> 22) 
In this analysis only the first (N/2) + 1 Fourier coefficients are considered. 
The first coefficient corresponds to the D.C. component of zeroth frequency 
and the N/2 + 1 coefficient corresponds to the highest frequency available. 
It is well known that the remaining coefficients correspond to the negative 
frequencies and are in fact the complex conjugates of the first (N/2 + 1) 
coefficients following the relationship g(N/2 * P) = g*(N/2 + 2 -P) for P = 2 
to N/2. However since all the coefficients are real, no new information is to 
be gained by examining the negative frequencies. 
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Once the appropriate number of coefficients, say I, has been deter- 
mined the remaining N-I coefficients of H(w) are set equal to zero before 
dividing H(w) point -by-point by R(w) and proceeding to obtain f(x) by the in- 
verse Fourier transform. The consequences of this truncation, or filtering 
operation are well known. See, for example, the discussion by Verkerk (1982). 
In general the process results in some systematic errors in which the peak 
values are underestimated and nonphyslcal oscillations may appear In the de- 
convoluted function f(x). For Doppler-broadened positron annihilation spectra 
taken with our instrumentation this filtering error is negligible for Doppler 
spectra having more than 40K peak counts. 

To check the accuracy of the function f(x) requires the convolution of 
this function with r(x) and the comparison of the resulting function with 
h(x). The details of practical convolution techniques have been discussed by 
Lam et al_^ (1981), and the procedure used in this work is essentially the 
same. 

A complete listing of the Fortran FFT/power spectrum deconvolution 
program is given in Appendix F. 

Test of the Method 

In order to explore the validity of the FFT/power spectrum method the 
followiny test was performed. The resolution function r(x) was obtained from 
the 514 KeV gamma ray of 85 Sr, while the function h(x) was measured for a neu- 
tron irradiated high purity single crystal sapphire sample. This sample had 
been irradiated in EBRII at 740°C to a fluence of 3 x 10*1 neutrons/cm* with E 
< U.l MeV. Both r(x) and h(x) were then processed as described 1n the prev- 
ious section. For this test of the method f(x) was estimated from F(w) usiny 
r t.hj^ii^ 
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tral analysis. 

Figure 5 is a plot of the power spectrum for both F(w) and H(w). The 
filled circles represent the spectrum of the experimentally obtained function 
h(x). Note the rapid decrease of about six orders of magnitude in power with 
increasing frequency. This power spectrum is about what would be expected 
from a wide-band-width stochastic process contaminated with white noise. This 
white noise component is associated with the band of frequencies to the right 
of the first minimum in the power spectrum of h(x). Although the overall 
signal -to-noise ratio is very high, for frequencies in this white noise band 
the signal-to-noise ratio is less than one and therefore it is impossible to 
retrieve any structural Information about the function f(x) in this high fre- 
quency range. This suggests that only those coefficients to the left of the 
white noise band, that is to the left of the first minimum, be retained. To 
use less than this number of coefficients in reconstruction is undesirable 
since it is known that the errors associated with truncation, or filtering, 
increase with a decrease in the number of coefficients used. 

The amplification of the statistical errors in h(x) by the deconvolu- 
tion procedure can be clearly seen 1n Figure 5. The open circles represent 
the spectrum of the deconvoluted function f(x). Note the drastic increase in 
the power associated with the high frequencies of f(x) as compared with the 
corresponding frequencies in h(x). However, the low frequency coefficients to 
the left of the first minimum are well behaved. These observations are the 
foundation on which the FFT/power spectrum method is based. 

It is not surprising to find that the number of retained Fourier re- 
construction coefficients, as determined by power spectrum analysis, is a 
function of the number of peak counts in the raw data h(x). Figure 6 is a 
^emi^-log plot of jthis reUtio ^shlg Jojlirradja ted sapphire. JThis. figure. sug- 
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Figure 5. A Plot of the Power Spectrum of the Raw Data, h(x), (Filled 
Circles), vs. the Deconvoluted Data, f(x), (Open Circles). 
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gests that the relationship over this range of peak counts is an exponential 
of the form: 

y <» m £n x (3.23) 
where y is the number of coefficients, x 1s the number of counts in the peak, 
and m is approximately equal to 7 for this data set. This curve demonstrates 
the high price one must pay in terms of data collection time to increase sig- 
nificantly the number of retained coefficients. Fortunately, our experience 
has shown that with the ~100 coefficients which are usually obtainable with 
40K counts in the peak (30 minutes of data collection) the results obtained 
are very satisfactory. 

Figure 7 is a comparison of the deconvoluted function f(x) and the raw 
data h(x) after background correction. The curves have been normalized to 
equal area. As expected, the deconvoluted curve 1s both sharper and smoother 
than the raw data. There appears to be a slight oscillation in f(x) in the 
tail region. The quality of this estimate of f(x) can be examined in the 
usual way by convolution with r(x) to obtain an approximation h'(x) to the 
original experimental function h(x) as is shown in Figure 8. The overall 
quality of the fit is so good that on this scale it Is difficult to see any 
differences between the raw spectrum and the convoluted spectrum. The one 
exception occurs in the peak region where several of the raw data points of 
h(x) lie well below their counterparts in h'(x). Careful examination of the 
shape of h(x) in this region reveals some rather spurious changes in the slope 
of the spectrum. Physical constraints on the annihilation process suggest 
that, assuming no posltronium formation, the annihilation spectrum should be 
relatively smooth in this region. This suggests that these slope changes in 
h(x) are related to statistical counting errors. To test this idea two other 



53 

One of the theoretical advantages of deconvolution is a significant 
improvement in the resolution of DBPAS measurements. Figure 9 is a plot of 
the high energy half of a typical 85 Sr peak before and after deconvolution. 
The closed circles represent the raw data after background correction, while 
the open circles show the shape of the resolution curve after deconvolution 
using 243 coefficients. This is the number of coefficients selected by the 
FFT/PS method corresponding to the first minimum in the power spectrum for 
this data set. Note that the full width half maximum (FWHM) value for the raw 
data is -1280 eV while the deconvoluted curve has a FWHM of -320 eV. This 
corresponds to a factor of four improvement in resolution. However, this 
result is somewhat misleading in that for aluminum the number of retained 
coefficients is approximately 60-70. The triangles in Figure 9 show the 
resolution function after deconvolution using 70 retained coefficients. In 
this case the FWHM is -576 eV corresponding to an improvement in resolution by 
a factor of -2.2. Clearly it can be seen that the improvement in resolution 
obtained by deconvolution using the FFT/PS method is a function of the number 
of retained coefficients, which in turn is directly related to the statistical 
quality of the raw data set (Schaffer et al. 1984). This implies that the 
data collection process should be carried out for sufficiently long times in 
order to obtain good counting statistics. (The counting error is proportional 
to 1/VN, where N is the number of counts in the channel). 

Unfortunately, counting for long times introduces another error 
source. This error source is related to instrumental electronic drift due, at 
least in part, to temperature fluctuations. There are two major ways in which 
other researchers have dealt with this problem. First, the accuracy of re- 
sults can be improved by operating the entire DBPAS system within a temper- 
- — - — ature-contFolled (*~U^2°C ) - atmospftet^-^ used in 
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Figure 9. A Comparison of t he Shape of the Hi gh Ener gy Half of _the 

After Deconvolution using 70 Coefficients (a), and c) After 
Deconvolution using 243 Coefficients, (o). 



conjunction with the first, is to digitally stabilize the electronics. This 
digital stabilization is accomplished by simultaneously monitering the 511 KeV 
annihilation line along with a suitably located monoenergetic gamma ray source 
such as the 478 KeV gamma ray from 7 Be. In this procedure a spectrum is 
collected for -100S and then the centroid of the 478 KeV peak is calculated. 
Subsequently, several other spectra are collected and their corresponding 
centroids also calculated. Next, each spectrum is linearly shifted such that 
all of the 478 KeV centroids coincide. Finally the 511 KeV peaks are summed 
to produce a DBPAS spectrum with as many total integrated counts as is 
necessary to obtain suitable counting statistics. 

While temperature and/or digital stabilization should improve the re- 
producibility of uBPAS measurements, it is still possible to obtain meaningful 
results using unstabilized equipment. In this work neither of the stabiliza- 
tion techniques were employed, and the major consequence of such a procedure 
is an observable spread in the values of the lineshape parameters calculated 
for repeated DBPAS measurements of the same metallurgical sample. This prob- 
lem will be discussed in greater detail in the sections of this chapter deal- 
ing with the separation of multiple defect types. 

In order to minimize the change in the metallurgical defect structure 
of the samples as a function of time, all of the DBPAS measurements in this 
work were made at 77°K. Performing the measurements at liquid nitrogen 
temperatures raises an important theoretical question: is the assumption that 
the positrons are essentially at rest, with respect to the electrons, at the 
time of annihilation still valid? Previous experiments by Kubica and Stewart 
(197b) demonstrate that even at temperatures as low as 4.2°K the positrons 
reach thermal velocities prior to annihilation. 
■ In order to^deterjjiine-the effect of-neasuring the annihilation 
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lineshape at 77°K an experiment was performed in which six DBPAS spectra were 
obtained for a well annealed pair of aluminum samples. Three of the measure- 
ments were made at 77°K and three at room temperature. The small deviation 
(of order 1/Vn) between the six spectra showed no correlation with measurement 
temperature. It was therefore concluded that decreasing the sample tempera- 
ture from room temperature to 77°K would not affect the positron-electron 
annihilation procedure in any experimentally observable manner. 

As mentioned in the experimental procedures the samples examined in 
this work were high purity aluminum disks heat treated to control defect type 
and concentration. It Is assumed in the remainder of this work, without 
direct TEM evidence but based on the arguments presented in the experimental 
section, that one group of samples contained predominately mono vacancies while 
a second group of samples possessed a defect structure dominated by prismatic 
dislocation loops. For each pair of disks three or four annihilation energy 
spectra were obtained, from which an array of lineshape parameters were sub- 
sequently extracted numerically. The lineshape parameters examined in this 
work include the first five even central moments of the energy distributions 
defined by the relationship: 

i N i 
ra(I) = W? J < x 1 " *> f < x i> 

where N is the number of channels in the spectrum 

x is the channel corresponding to the spectrum centroid 
f(x-j) is the number of counts in channel i 

and I is an even integer (2,4,6,8,10) corresponding to the first 
through the fifth even central moments 

A second group of parameters is defined by the relationship: 
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j+4 x+59 
"(I) • 11 f(x,)]/ I f( Xl , (4il) 



with j defined as j = x + 5(1-1) 

where f(x,) is the number of counts in channel i 

x 1s the channel corresponding to the spectrum centroid 
and I i s an integer from 1 through 12 

For example, x(l) is the ratio of the number of counts in the first five 
channel regions on the high energy side of the centroid (actually it includes 
the centroid) to the total number of counts in the high energy half of the 
spectrum. Similarly, x(12) is the ratio of the number of counts in the 
twelfth five channel interval (to the right of the centroid) to the total 
number of counts in the high energy half of the spectrum. 

The third, and final, group of l.ineshape parameters examined are the 
conventional "S» parameters (MacKenzie et al. 1970). These parameters are 
defined as the ratio of the number of counts in some central energy region to 
the total number of counts in the energy spectrum. The width of the energy 
window used to define the S parameter is selected so that the value of S is in 
the range 0.4 - 0.6 (Smedskjaer and Fluss 1983). I n this study three 
S-parameters (S(1344), S(1472), and 3(1600)] were examined with respective 
central energy windows of 1.344, 1.472, and 1.600 KeV. 

These parameters were calculated for both the background corrected raw 
data and for the deconvoluted data. Hence a total of forty lineshape 
parameters were calculated for each UBPAS measurement. 

Figure 10 is a plot of the parameter x(l) for the deconvoluted data. 
The open circles represent the values of this parameter for the "as quenched" 
samples (predominately monovacancies). The filled circles correspond to the 
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low temperature anneal samples (prismatic dislocation loops), and the tri- 
angles represent the samples after the high temperature anneal (the essen- 
tially defect free state). Note that the numerical value of x(l) increases 
with quenching temperature for both defect types. This trend 1s in agreement 
with theoretical predictions since an increasing defect concentration results 
in a higher probability of annihilation with lower kinetic energy conduction 
electrons, corresponding to a more narrow energy spectrum, and hence a larger 
value for x(l). The equations used to predict the theoretical defect concen- 
trations for both the vacancies and the prismatic dislocation loops as a func- 
tion of quenching temperature are discussed in Appendix 6. 

Figure 10 demonstrates that DBPAS measurements are sensitive to the 
presence of both vacancies and dislocations in high purity aluminum. This 
result is, of course, a necessary prerequisite for the possibility of DBPAS 
measurements being able to distinguish between the defect types. From this 
plot it is possible to perform some semiquantitative analysis and estimate 
some physical properties of aluminum. For example, one can estimate values of 
the "cross section* 1 for the trapping of a positron by both vacancies and 
dislocations in aluminum. The accuracy of these estimates is limited by the 
reproducibility of the DBPAS measurements as well as the need for some 
independent information such as the lifetime of positrons as a function of 
defect type. 

In 1972 Cotterill et al_. published a paper in which they examined the 
average positron lifetime, r avg , on a set of aluminum samples very similar to 
the ones studied in this work. Their work clearly demonstrated that positron 
lifetime experiments were sensitive to both vacancies and dislocations in 
aluminum. Before discussing the similarities between their work and the 
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the parameters T avg and x(l). In a review article by Sieger (1980) it is 
shown that: 

T avg s (1 - ) T b + Vt (4.3) 
where x b is the lifetime in the bulk material 
T t is the lifetime in the trapped state 
and A t is the probability that the positron annihilates from the 

trapped state 
Similarly, Sieger demonstrates that: 

F Ml - A t )F„ + AtF t (4 . 4) 
where F is any characteristic shape parameter obtained from the analysis 
of DBPAS measurements 
F D is the value of the parameter in the bulk state 
and F t is the value of the parameter in the trapped state 

Since the parameter x(l) can be shown to be equivalent to an "S" type 
parameter, and since an "S» type parameter is a common characteristic line- 
shape parameter, Eq. (4.4) can be rewritten as: 

X(l) - (1 - A t )x(l) b + A t x(l) t ( 4 . 5) 
where x(l) b is the value of the parameter in the bulk state 
and x(l) t is the value of the parameter in the trapped state 
Clearly a comparison of Eq. (4.3) and Eq. (4.5) demonstrates the similarity 
between x avg and x(l). 

With this in mind it is not surprising to find that plots of x avg vs. 
quench temperature (T Q ) and x(l) vs. T Q share many common characteristics. 
(The T avg plot used for comparison Is Figure Z in the work of Cotterill et 
aU. 1972). In both figures the Tq at which the vacancy concentration is 
large enough to affect the measured parameter is ~200°C. The saturation value 
°! b0th l av Li n ^ Vf_600 o C, , The, curve. through 
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the points corresponding to vacancies is sigmoidal in shape for both the 
lifetime and deconvolved DBPAS measurements. As for the points corresponding 
to samples with prismatic dislocation loops, the value of both parameters 
initially deviates from their bulk values for samples annealed after a quench 
from ~4UU°C. Neither parameter reaches a saturation value for dislocations 
within the temperature range examined in this work. 

Figure 10, therefore, shows not only the ability of deconvoluted DBPAS 
measurements to detect the trapping of positrons by both vacancies and dislo- 
cations in aluminum, but also qualitatively verifies the theoretically pre- 
dicted equivalence of the lifetime parameter T avg and the Doppler parameter 
x(D. 

Figure 11 is a plot of the same lineshape parameter, x(l), vs. Tq for 
the background corrected raw data, i.e. the same data as shown in Fig. 10 but 
before the deconvolution procedure. Note the general similarity between these 
two figures. Careful examination of the raw data shows that the points 
corresponding to the 4O0°C quench are systematically displaced. That is, the 
entire group of five data points seems to be "too high" as compared with the 
rest of the data. To understand the cause of this shift it is necessary to 
examine the shape of the instrument resolution function R(x). 

During the course of these experiments the Instrument resolution was 
periodically checked by collecting an energy spectrum for 85 Sr. The monoener- 
getic line is broadened due to the physical characteristics of the DBPAS sys- 
tem as previously described. This broadened energy line serves as the resolu- 
tion function, K(x), in the FFT/PS deconvolution procedure. It also describes 
the amount by which the undistorted energy spectrum is broadened by the in- 
strumentation. 
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Figure 12 shows two such resolution spectra. The open circles repre- 
sent the instrument resolution corresponding to the DBPAS measurements for the 
samples quenched from 400°C. The filled circles represent a resolution func- 
tion which is characteristic (within statistical error) of the resolution 
functions for all of the remaining data points. This comparatively large 
shape change in the function R(x) was apparently caused by a power failure In 
the laboratory. Note that the Function R(x) corresponding to the 400°C 
samples (henceforth termed SR8S0LD) is more sharply peaked, and therefore more 
narrow, than the other resolution function (SR85NEW). 

The physical significance of the change in resolution function is as 
follows. An observed energy spectrum will be broadened in two ways. The 
first way is by a decrease in the concentration of vacancy type defects, in 
the sample being examined. The second broadening mechanism is, of course, the 
broadening due to the instrumentation. Now, if the second source of broaden- 
ing remains essentially constant, then changes in the width of the observed 
spectra may be attributed to changes in defect concentration. On the other 
hand, if the broadening due to the instrumentation decreases, and this error 
source is not properly removed, the observed spectrum will be artificially 
narrow and its apparent defect concentration will be too great. Thus, if the 
data for the 4Ul)°C samples are collected with an artificially narrow resolu- 
tion function, i.e. SR850LD, their apparent defect concentration will be 
greater than they should be. 

Of particular significance is, that while a change in experimental 
conditions can alter the shape of the resolution function and thereby distort 
the values of lineshape parameters calculated for the raw data, the deconvo- 
lution routine, when supplied with the appropriate R(x) function, auto- 
■ m i t l5Jjly_ylel d A the,deiireijj^ that for reproduca- - 
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Mlity the deconvolution of DBPAS measurements should be considered. In 
addition, as mentioned previously, since deconvolution removes the instru- 
mental effects it should yield lineshape parameters which are directly 
comparable between different laboratories. 

Despite the advantages of deconvolution it 1s important to remember 
that this technique does, of course, have its limitations. In DBPAS 
measurements there are several sources of experimental error, each of which 
has a corresponding time scale. The magnitude of this time scale is the 
overriding factor in determining the applicability of deconvolution techniques 
in minimizing the specific error in question. For example, any change in 
experimental conditions which can be treated as a step function, that is, by 
changing from one constant value to another constant value in a very short 
period of time, can be easily treated by deconvolution. Such "error" sources 
include the accidental (e.g. power failure) or intentional changes in the 
electronic amplification of the energy signal. To correct for such a change 
one simply needs to collect a resolution function before and after the event 
and use the appropriate R(x) in the deconvolution algorithm. 

On the other hand, error sources having a time scale approximately 
equal to the time required to collect a data set (tens of minutes) can not be 
corrected by deconvolution alone. Examples of such possible error sources 
include line voltage fluctuations, temperature variations, and electronic 
drift. These error sources may cause the centroid of the resolution function 
to vary during the collection period, resulting in a composite R(x) which is 
artificially broad. In addition, even if the raw data set to be deconvoluted 
and the function R(x) are collected in rapid succession the amount of centroid 
drift for the two energy spectra may well be uncorrected. This type of error 
SQuxce^rei}uXc.e^^ of correct i on .-The- DBPAS measurements - 
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should be performed in a temperature controlled laboratory. Also, digital 
stabilization, as described previously, should be used to control cent ro id 
drift. 

The final error type is the small time scale or hiyh frequency class. 
An error source such as sixty-cycle electrical noise, if it were present, 
should average to zero over the time required to collect a data set. In 
addition, a high frequency error source should have a similar effect on both 
the resolution function and the DBPAS energy distribution. 

In order to judge the effectiveness of deconvolution in eliminating 
experimental error it is necessary to define some quantitative measure of this 
error. The predominant manifestation of experimental error appears as a 
spread in lineshape parameter values calculated for multiple DBPAS measure- 
ments on the same samples. Figure 13 is a plot of the parameter x(l) vs. Tg 
for the deconvoluted data, note the spread in values of x(l) for a given 
defect configuration. 

In an attempt to evaluate the performance of the deconvolution algor- 
ithm the following parameters are introduced (see Figure 13 for a physical 
interpretation). 

1) *Pv/ p o 1s a measure of the sensitivity of this parameter to the 
existence of vacancies. 

2) 6P V /Ap y is a measure of the error associated with the spread in 
values of this parameter for vacancies. 

3) APj/P 0 is a measure of the sensitivity of this parameter to the 
existence of dislocations. 

4) &P i /i^P 1 is a measure of the error associated with the spread in 
values of this parameter for dislocations. 
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Clearly, it Is desirable for the sensitivity (AP/P 0 ) to be as large as 
possible and the error (6P/AP) as small as possible for both defect types. A 
comparison of these quantities calculated from the deconvoluted and the raw 
data is shown in Table 1. Note that for all four of the lineshape parameters 
listed the sensitivity ratio is higher and the error ratio is lower for the 
deconvoluted data than for the raw data. This analysis encourages the further 
use of the deconvolution procedure to improve the quality of the results ob- 
tained from D8PAS measurements. 

A tentative explanation for the apparent improvement in sensitivity 
and accuracy of the deconvoluted results is based upon trends in the shape of 
the energy spectrum before and after the application of the FFT/PS analysis. 
Figure 14 shows the high energy half of the peak region for an M sample 
quenched from ~600°C. The deconvoluted spectrum is, of course, more narrow 
than the raw data but note that it is also smoother (due to the filtering of 
the high frequency components). One explanation for the comparatively large 
error ratios of the raw data is that they are related to the statistical 
counting errors in the peak region of the spectrum. If this is indeed the 
case it is not surprising to observe a decrease in the error ratio after the 
filtering of the high frequency components. However, the most probable 
explanation for the increased sensitivity resulting from deconvolution is 
related to the removal of the "broadening" action of the experimental instru- 
mentation. Recall that the major effect of the instrumentation is to spread 
counts from one channel into several adjacent channels. Therefore, the true 
lineshape changes in the annihilation energy distribution, caused by the 
trapping of positrons at defects, are spread out over a much larger number of 
channels in the raw data. The deconvolution procedure corrects for this 
ef_f e.ct- and in. so ..doing, restorer some :of. the details-of:; the—true"- annihilation — 
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energy distribution. The net result is a „ i„ cre ase in defect sensitivity. It 
should be mentioned that the magnitudes of the enhancement in sensitivity and 
accuracy due to deconvolution noted in this work are consistent, but small. 
However, from this work alone it is not possible to estimate the effect of 
temperature and electronic drift upon the results. It would be extremely 
enlightening to repeat this analysis employing a temperature controlled and 
digitally stabilized system. 

Before moving on from the topic of parameter sensitivity and accuracy, 
one more pertinent observation should be discussed. Figure 15 is a plot of 
the error ratio, 6P/ AP , calculated separately at each quenching temperature 
(the parameter being examined is x(l)). Several trends are apparent. First, 
it is clear that the error associated with the dislocations is larger than the 
error associated with vacancies. This trend may be a. direct consequence of 
the size of the denominator, *P, which is smaller, at a given quench tempera- 
ture, for samples with dislocations than for the corresponding samples with 
vacancy type defects. This observation is however somewhat deceptive in that 
the horizontal axis of Figure lb is T Q and not defect concentration. A mean- 
ingful comparison of parameter sensitivity between defect types should be made 
using equivalent defect densities. Unfortunately, the dislocation densities 
which can be introduced by this thermal treatment are not large enough for 
direct comparison. In addition, it should be mentioned that the trapping of 
positrons by dislocations may be somewhat more complicated than trapping by 
vacancies. This idea is discussed in detail by Martin and Paetsch (1972). 
Very briefly, it may be that the depth of the potential well associated with a 
prismatic dislocation loop is such that detrapping of the positron, due to 
local thermal excitations, becomes significant. If this is indeed the case, 
the interpretat^ v: 



o 
- o 



o 
o 

LO 



It) 



to 

f- 

o 
c 
m 
u 

> 

c 

to 





to 




c 


o 


o 


o 






at 


0) 


u 


&. 


o ♦ 




t/1 


4-> 




<0 




O 






°i 


*t a. 


-C <D 


B 


■P CO 


0) 


o 


*— 




jC 


U C 


o 


o 


c 




0) 




3 




cr 


3P/AP 

for A 



o 
o 

CO 



o 
o 

CM 



a> ^ 

CU to 
B U 

to a> 
I— 

u 
c 

<_> 
V c 
LxJ 0) 

a* cr 



o c 
o 
+-> .1- 

o 

f— u 

o. c 



o 
o 

LO 



o 
o 



o 
o 

CO 



o 
o 

CM 



o 
o 



(SltUf) AUwnqjy) dV/de 



73 

becomes more complicated. Clearly, additional experimental data Is required 
to determine the exact reason for the observation that the error parameters 
are larger for the samples containing primarily dislocation loops than for the 
samples dominated by vacancy type defects. 

Another interesting feature of Figure 15 is the relatively constant 
value of 6P/AP for vacancies. This implies that the spread 1n the data de- 
creases with the quenching temperature, since clearly AP decreases with T Qi 
One possible explanation for this observation is related to the shift in the 
centrold of the experimentally obtained spectrum due to temperature and/or 
electronic drift. The explanation is that a centrold shift changes the peak 
shape of a broad and flat spectrum (low defect concentration) to a lesser ex- 
tent than that of a narrow and sharp spectrum (high defect concentration). 
This is, of course, an extremely qualitative argument but it does predict that 
in a temperature/digitally controlled series of experiments the relative size 
of the error bars associated with the high defect concentration samples should 
decrease. 

The preceding discussion has demonstrated the utility of deconvolution 
in DBPAS measurements. In summary, the use of the FFT/PS method improves the 
sensitivity and the relative error of the calculated Uneshape parameters. In 
addition, deconvolution automatically corrects for a wide range of error 
sources which result in changes in the instrument resolution function. Decon- 
voluted lineshape parameters should be more directly comparable between dif- 
ferent laboratories and, in addition, offer the potential ability to resolve 
fine structure in the energy spectrum which Is unavailable from direct obser- 
vation of the raw data. There is strong evidence (see for example Hood and 
Schultz, 1980) suggesting that the use of both temperature stabilization and 
.. d . 1 altal stajbi ] i zaUqn,. tp._.c.ontr.ol centrold drift*, should slgnlf^rant-^rrrlaprffyi' -„„ 
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It is Interesting to note the conditions under which the undeconvo- 
luted raw data has yielded reasonable results in the past. The main require- 
ment is, of course, experimental conditions which yield a nearly constant 
resolution function, R(x). As one would expect, temperature and digital con- 
trol as well as performing the series of experiments as close together in time 
as possible, both tend to help yield reasonable results for the raw data. In 
addition, the use of "S" type parameters, which are essentially area integrals 
over large energy windows, will be affected only slightly by the "spreading" 
action of the instrumentation. The major advantages one gains by deconvolu- 
tion are: 1) the ability to compare more accurately results for a wide var- 
iety of experimental conditions and, 2) the possibility of observing a more 
detailed structure of the positron -electron annihilation energy spectrum. 

The Separation of Multiple Defect Types 

The previous discussion has demonstrated that deconvoluted DBPAS 
measurements are able to detect the presence of both vacancies and disloca- 
tions 1n high purity aluminum. The evidence for this statement comes from 
controlled laboratory experiments using samples with reasonably well charac- 
terized defect distributions. In order for DBPAS measurements to become use- 
ful as an NOT tool in industrial applications it 1s important for the tech- 
nique to be able to distinguish between multiple defect types. The reason for 
this requirement is that in a practical application one may require informa- 
tion about a particular defect type while being comparatively unconcerned with 
other defect types. For example, if one is interested in hydrogen embrlttle- 
ment of high strength steels the "masking" of vacancies by nascent hydrogen 
may be of primary importance while slight variations in dislocation density 
may have relaUve^m " : 
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Before discussing the results of this investigation to determine the 
ability of DBPAS to separate defect types a brief review of previous work in 
this area will be presented. The basic concept Is that no two defect types 
will have identical electron configurations associated with them. That is, 
the distribution of electrons (both in terms of electron density and momentum) 
in the neighborhood of a single vacancy will be different than in the neigh- 
borhood of a dislocation core. Qualitatively the atmosphere will be the same, 
i.e. a decrease 1n electron density and a lower average electron momentum, but 
quantitatively, the measured lifetimes and fine structure of the momentum dis- 
tribution should be different. 

This variation in electron distribution has been verified repeatedly 
in the case of lifetime measurements. For example, Fluss et aj_. (1978) report 
a lifetime of -244 ps for positrons trapped at a vacancy in aluminum while 
Hautojarvi et aU (1970) give the value of -229 ps for a dislocation in the 
same material. It should be mentioned that the accurate determination of 
multiple lifetimes from a single experiment requires some rather involved cal- 
culations. The complexity of the procedure Increases rapidly with the number 
of unconstrained components one attempts to resolve and, for crystalline 
materials, the results are questionable if more than two defect types are 
considered. 

An alternate method of analyzing lifetime data is to calculate a 
single lifetime component, known as the "average" lifetime, for a sample con- 
taining multiple defect types. This approach simplifies the calculations and 
finds its greatest utility in experiments in which the concentration of one 
defect type is varying over a wide range of values while the concentrations of 
the other defect types remain relatively constant. Of course, such an 
. .-- -gMroach results in the loss of the ability to separate defect - types. _ 
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In the case of OBPAS measurements there have been fewer attempts to 
distinguish between different defect types. The most notable exception is the 
work of Nantl and Triftshauser (1978) which will be discussed in a following 
section. There are at least two major reasons for the lack of work 1n this 
area. The first reason is related to the comparatively poor energy resolution 
of the DBPAS technique. Because of the smoothing action of the instrumenta- 
tion the fine structure of the energy distribution may have been hidden in 
previous experiments. In addition, even a slight change in the instrument 
resolution function during a series of experiments could mask the subtle 
changes In lineshape required to separate defect types. Although temperature 
and digital stabilization of the electronics would Improve the situation, such 
a procedure does not address the issue of the smoothing action of the instru- 
mentation. It is the use of deconvoluted data which offers the possibility of 
examining the fine structure of the electron momentum distribution. 

As mentioned In the introduction, the most conmonly reported Doppler 
lineshape parameter is the MacKenzie "S" parameter (1970). This lineshape 
parameter is defined as the ratio of counts in some central region of the 
DBPAS energy spectrum to the total number of counts in the spectrum. Since it 
is known that the energy spectrum becomes more sharply peaked with increasing 
vacancy type defect concentration 1t follows that the numerical value of "S" 
should be proportional to the total vacancy type defect concentration. The 
disadvantage of using this lineshape parameter is that it is unable to 
distinguish between multiple defect types. For example, if both monovacancies 
and edge dislocations are present in a sample both defect types will tend to 
increase the value of "S\ however, it will be impossible to separate the 
relative defect concentrations solely from the calculated value of tt S\ 
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Mantl and Trifthauser (1978) developed a method for separating the 
effects of multiple defect types from experimentally obtained OBPAS data. 
To accomplish this they introduced a new lineshape parameter, I c , defined 
as the ratio of the sum of the number of counts in two regions, (one on the 
low energy side and one on the high energy side of the peak) divided by the 
total number of counts in the spectrum* This lineshape parameter, which in 
this work will be referred to as the winy or "W" parameter, can be shown to 
decrease with increasing vacancy type defect concentration. When analyzed 
separately the "W" parameter suffers from the same disadvantages as the "S" 
parameter in that it is incapable of separating multiple defect types. How- 
ever, Mantl and Trifthauser recognized that the combined use of "S" and "W" 
could yield useful results. 

In Appendix H the derivation of their defect specific "R" parameter is 
recreated. To summarize, it can be shown that in theory the W R M parameter, as 
defined below, is concentration independent, but characteristic of the type of 
trapping site involved (Mantl and Trifthauser 1978). The "R" parameter is 
defined as: 

R - |(S - S B )/(W - W B )| (4*6) 
where S is the "S* parameter for a given specimen, 

Sg is the value of "S" corresponding to a "defect free" specimen, 

W is the "W" parameter for a given specimen, 

W B is the value of "W" corresponding to a "defect free" specimen. 
While such a parameter is very useful, because of its ability to distinguish 
between multiple defect types, it still does not solve the problem completely 
due to some inherent theoretical limitations. The main restriction is that 
the analysis is valid if only one type of trapping center is present, or if 
one defect type is essentially., dominant (Mantl and Trifthauser 197? )_. In 
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practice this means the experimenter has the ability to calculate defect type 
and concentration provided each sample has only one predominant type of 
positron trap. The difficulties arise if a sample contains two defect types, 
say vacancies and dislocation cores, in roughly equal concentrations. 

In terms of the present experiments on high purity large grain size 
aluminum, because care has been taken to produce samples with only a single 
dominant defect type, one would expect the "R w parameter analysis to yield 
useful results. A family of R-type parameters can be calculated for the 
experiments in this work by examining ratios of the form 

R N " l(Ux(i)/Ax(j))| (4. 7) 
where R N is an R-type parameter 

Ax(i) is the difference between x(i) for the current sample and the 
value of x(i) in a defect free sample 
and i is less than j 

In the original work of Mantl and Trifthauser (1978) the energy windows 
selected to define "S" and M W M were chosen so that the annihilation events 
contributing to the M S" parameter corresponded to conduction (valence) 
electrons while core electron events resulted in contributions to the "W M 
parameter. This relationship restricts the possible values of i and j as 
defined in Equation (4.7) above. 

The specific rt R" parameter used in this study is defined as 
0 . 1 Ax(l) I 

K l ~ Kx(10) + x(ll) + x(12)J| < 4 - 8 > 
and is defined such that the energy windows are fairly close to those used by 
Mantl and Trifthauser (1978). For comparison, Rj has a low energy window of 
-511 - bll.32 KeV and a high energy window of -513.88 - 514.84 KeV compared to 
windows of -511 - 511.5 KeV and -513.6 - 516.1 KeV in the referenced work. 
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Figure 16 1s a plot of Kj as a function of quenching temperature for 
the background corrected raw data. Note that although samples with primarily 
dislocation type defects have, in general, a higher R x value than the samples 
with primarily vacancy type defects, there is a significant amount of scatter 
in the data. The source of this scatter is not known exactly but it is prob- 
ably related to the fact that the data was collected without temperature or 
digital stabilization. In fact, it is plausible that the spread in Hi values 
is a numerical amplification of the spread in values of x(l), x(10), x(ll) and 
x(12). In addition to this error source, any error associated with changes in 
the instrument resolution function, (IRF), should also be present in Figure 
16. To test the utility of deconvolution in R-parameter analysis one would 
like to have a data set in which the second error source is present but the 
first error source, related to electronic drift, is negligible. Unfortunate- 
ly, such a data set can not be obtained with the instrumentation used in this 
study. 

In an attempt to isolate the instrument resolution effects the fol- 
lowing procedure was employed. Using the fact that all DBPAS measurements for 
a given metallurgical sample were performed In a comparatively short time 
interval (-8 hrs.), it seems reasonable to assume that the spread in calcu- 
lated lineshape values is a result of the lack of temperature and digital 
stabilization rather than of changes in instrument resolution. This hypoth- 
esis is supported by the observation that for a given set of consecutive 
experiments, the IRF's measured before and after the experiments showed only 
very slight changes. If this assumption is true, one may be able to arti- 
ficially select a subset of the original data which qualitatively corresponds 
to a series of experiments with reduced electronic drift error. In this work 
each, .group of_-BB£A5-.meASjAra^^ subjected,- to„„ .-^ 
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a given thermal treatment, was Independently examined and an average value for 
each relevant Uneshape parameter was calculated. The corresponding Rj values 
calculated from the average Uneshape parameter values are shown in Figure 17. 
These values are, as 1n Figure 16, for the background corrected raw data. 

A similar procedure was then employed for the deconvoluted data. That 
is, average x(i) values were calculated for each specimen and then used to 
calculate R\ values. These Rj values for the deconvoluted data are shown in 
Figure 18. It is quite apparent that the spread in R-parameter values has 
been increased by the FFT/PS deconvolution algorithm. The explanation for the 
increased scatter in the data is not immediately obvious. In a previous 
section of this work It was demonstrated that deconvolution improved the 
quality of the calculated x(l) lineshape parameter values. That is, the 
spread in x(l) values for DBPAS measurements performed on a given sample was 
decreased (refer to Table 1). One could argue that decreasing the scatter of 
x(l) for a given sample should improve, or at least not adversely affect, the 
amount of scatter in the calculated Rj values. The difficulty is, of course, 
that Rj depends not only on x(l) but also on the values of the parameters 
x(10), x(ll), and x(12). These parameters are defined in the tail regions of 
the energy distribution and it is in this region that one must look to under- 
stand the difficulty associated with the calculation of R-type parameters for 
FFT/PS deconvoluted data. 

Figures 19 and 20 are a comparison of the energy distribution (raw 
data) for a sample quenched from 600°C and then subsequently annealed at 
523°K. Figure 19 shows the differences in the peak region of the spectra. 
Recall that x(l) is defined over the energy window of 511 - 511.32 KeV. 
In Figure 20 the vertical scale has been expanded so that the differences 
th e tai 1 re gion axe xtMW e_. fe.ajtjj.f^s^oJLdto 
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distribution in the tall region are apparent. First, it must be noted that in 
general the data points for the quenched sample He below their counterparts 
for the annealed sample as predicted by theory. Secondly, the statistical, or 
counting, error from channel to channel can be recognized. The former obser- 
vation is believed to account for the qualitative agreement between the- 
oretical predictions and the calculated R x values for the raw data. The 
latter observation, statistical noise, is the most probable explanation (along 
with electronic noise) for the spread in R 2 values as shown in Figure 16. 

For comparison. Figures 21 and 22 show the change in energy distribu- 
tion for the same two samples after deconvolution using the FFT/PS algorithm. 
In Figure 21, which shows the peak region, some of the desirable features of 
deconvolution are evident. Most notably, the statistical channel to channel 
noise has been removed (compare with Figure 19). A careful comparison of the 
raw and deconvoluted data shows that deconvolution does indeed decrease the 
full width half maximum as discussed previously. 

Figure 22 shows, with an expanded vertical scale, an undesirable 
feature of the FFT/PS deconvolution method. It was mentioned in Chapter III 
that truncating the Fourier coefficients is equivalent to the insertion of a 
low pass filter and it is well known that low pass filtering Introduces 
certain errors. Recalling that the Fourier analysis of the convolution 
equation, neglecting the statistical counting error, yields the following 
equation for the desired function f(x): 

f(x) = / u(w)F(w)e- iw x dw (4. 10) 

where u(w) = 1 for w < w c 
- 0 for w > w c 
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and w c corresponds to the frequency which is the first minimum of the 

power spectrum of h(x) 
It can be shown (Bendat and Piersol 1971) that error introduced by such a 
filter is proportional to (l/w c )sin(2irxw c ). The significance of this result 
is that the amplitude of the error is inversely proportional to w c (the number 
of retained coefficients), and the frequency of the resulting error is 
proportional to w c . 

This effect is clearly evident in Figure 22. The open circles 
represent the annealed sample for which the power spectrum analysis yielded a 
w c of fifty six (56) coefficients. For comparison, the filled circles, 
corresponding to the quenched specimen are calculated using sixty seven (67) 
retained coefficients. It is quite obvious that the error amplitude is 
decreased for the sample with the larger number of retained coefficients. A 
careful analysis of the data shows that the frequency of the error is 
proportional to w c . 

With this information, it is now possible to explain the increased 
scatter in Rj values displayed by the FFT/PS deconvoluted data. Because the 
energy window used to define the wing parameter (513.88 - 514.84 KeV) is in 
the range where the amplitude of the filtering error is of the same order of 
magnitude as the amplitude of the true signal, this error source is 
significant. If the wavelength of the error source was small compared to the 
size of the energy window, then the error contribution would average to nearly 
zero and this difficulty would not exist. The problem arises because the 
energy window spans less than a period of the sinusoidal error source and it 
spans a different fraction of the period depending on the number of retained 
Fourier coefficients. For example, in Figure 22 the energy window of interest 
.f of, the -4/inea .1 i ed sampj e spans ,apprpximatelj^ ^Ate-poiltJve-^ajf. eye] & &f _the_ _ „ 
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sinusoidal filtering error. This leads to a calculated wing parameter with an 
artificially high value, which 1n turn distorts the value of the calculated Rj 
parameter. Since multiple DBPAS measurements on even the same metallurgical 
sample may yield slightly different numbers of retained coefficients, this 
error source can lead to a scatter of Rj parameters for a single sample. 
Because the number of retained coefficients varies more widely between samples 
of differing defect density (or defect structure), the spread in Rj values 
becomes larger for such samples. 

Having related the scatter In values for the deconvoluted data to 
the phase and amplitude of the filtering error, it is now possible to try to 
correct for this error source. There are at least three ways to proceed. 
First, one could simply extend the experimental collection times to get a 
higher number of total integrated counts in the energy spectrum. This would 
result in better counting statistics, and a higher number of retained 
coefficients (Schaf f er et aJL 1984). Since an increased number of retained 
coefficients decreases the amplitude of the filtering error, such a procedure 
should yield better Ri values. Unfortunately one must pay a very high price 
in terms of collection time to significantly increase w c . 

Another approach is to use the same number of retained coefficients 
for each data set. That is, one could calculate a w c value for each data set 
to be analyzed and then decon volute all of the data sets using the minimum 
value of w c . Using this procedure will not eliminate the error source, but it 
will be controlled so that its contribution to the calculated values of R\ 
does not vary between experiments. Clearly such an approach has several 
inherent disadvantages. The primary disadvantage 1s that the absolute values 
of Ri will be a function of the particular value of w c used in the deconvolu- 
t i on procedure Qne_.wou_ld. prefer _an;: approach^whi ch-el;i minates , or : at -least 
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greatly reduces, the filtering error rather than simply controlling it. A 
method for accomplishing this goal will be discussed in the final chapter of 
this work. What follows here is a discussion of controlling the error source 
in the tail region by using the same number of retained coefficients for each 
data set. 

Figure 23 shows the energy window used to calculate the wing 
parameter. The vertical scale has been greatly expanded to show the fine 
detail of the structure. This plot is for the same samples as in Figure 22, 
but in this case both data sets have been deconvoluted using fifty seven (57) 
coefficients. As expected, the filtering error Is 1n phase. In addition, as 
theory predicts, the data points corresponding to the quenched sample are 
below the points representing the annealed sample. One can extend this 
procedure to all of the other data sets and then calculate the corresponding 
Rl parameters in the manner previously described. The results are presented 
in Figure 24. When this figure is compared with Figure 18 it can be seen that 
the filtering error does influence the scatter in Rj values. By simply 
controlling this error source, I.e. by using the same number of retained 
coefficients in the deconvolution of all data sets, the quality of the 
calculated values can be improved. 

To summarize this section, the average values of R x , along with the 
associated standard deviations for each of the two defect types, have been 
calculated for four different types of data. The four data types are: 1) the 
entire raw data set (Fig. 16), 2) the average raw data (Fig. 17), 3) the 
averages of the deconvoluted data as obtained directly from the FFT/PS 
algorithm (Fig. 18), and 4) the averages of the deconvoluted data using the 
same number of retained coefficients for each data set (Fig. 24). These 
values are jhpwn in_Table_ II.- ~. _ 
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o(D1s.) 



1 


.261 


.458 


.0677 


.1751 


2 


.240 


.440 


.0347 


.0492 


3 


.234 


.364 


.1515 


.1876 


4 


.253 


.424 


.0216 


.0290 



Note that the largest standard deviations occur for the deconvoluted 
data when the maximum number of coefficients are used in the reconstruction of 
each f(x). This result is explainable in view of the previously discussed re- 
lationship between w c and the amplitude and phase of the filtering error. 
When this error source is controlled by using a single value of w c for the de- 
convolution of all data sets (data type 4) the calculated standard deviations 
are improved. With this modification to the FFT/PS deconvolution algorithm it 
appears that if one defect type dominates, one can distinguish between 
vacancies and prismatic dislocation loops, in high purity large grain si2e 
aluminum. 

To conclude the results and discussion chapter of this work the even 
central moments of the annihilation energy distributions will be discussed. 
As mentioned in the introduction, because this technique has the potential to 
distinguish between multiple defect types, this family of lineshape parameters 
is interesting from a theoretical point of view. Briefly, this possibility 
arises because the higher order central moments are increasingly dominated by 
the tail regions of the energy distribution. Thus, the low moments give 
information about the central region (S-type information) while the higher 
order moments yield W-type information. Unfortunately, this analysis is 
extremely sensitive are-greatly--- ™ 
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amplified by moment analysis. Therefore the statistical error ruled out 
moment analysis in the raw data while filtering error precluded its use in the 
deconvoluted data. It is possible, however, that an improved filtering 
technique, discussed in the following chapter, will allow moment analysis to 
become a useful tool in future work. 



CHAPTER V 
CONCLUSIONS AND FUTURE WORK 



Over the past several decades positron annihilation spectroscopy has 
evolved Into an accepted tool for the study of the defect structure of mater- 
ials. Of the three most common measurement techniques Doppler broadening of- 
fers the advantages of reasonably high data collection rates, comparatively 
low source strength requirements, and few geometrical experimental con- 
straints. It is therefore the technique with the best chance of making the 
transition from a laboratory research device to an industrial NOT tool. Un- 
fortunately, DBPAS measurements suffer from some inherent limitations. The 
two most significant limitations are comparatively poor instrumental energy 
resolution and limited ability to distinguish between multiple defect types. 
Based on these observations, the two primary goals of this work were: 1) to 
develop an efficient, model independent, noniterative deconvolution method, 
and 2) to improve the ability of DBPAS to distinguish between multiple defect 
types. 

In order to improve the effective energy resolution of the DBPAS in- 
strumentation the fast Fourier transform/power spectrum analysis deconvolution 
algorithm was developed (Schaf fer et 1984). The convolution integral 
equation was shown in Chapter III to describe the relationship between the ex- 
perimentally observed energy distribution, the instrument resolution function, 
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and the "true" or desired Ooppler broadened eneryy distribution. Fourier 
analysis was then used to obtain what appeared to be a closed form solution 
for the desired function. Unfortunately, it was demonstrated that this 
solution is not unique. The statistical nature of the experiment places the 
equation into the mathematical class of incorrectly posed problems. In an 
attempt to find a reasonable approximation to the desired function the 
dispersion of the apparent closed form solution was investigated. It was 
observed that the dispersion became unbounded if the power spectrum of the 
statistical experimental "noise" remained finite for high frequencies. That 
is, the mathematical difficulties are related to the "white noise" component 
present in the experimental data. 

Since the shape of the power spectrum of the raw data appeared to be 
of critical importance, particularly if it contained a "white" component, it 
was investigated In greater detail. From this analysis it was observed that 
the overall signal-to-noise ratio for the experiment was very high, being of 
order ~10 6 . Using the observation that the power spectrum falls steeply to a 
minimum before rising due to the white noise component in the experimental 
data, the FFT/PS method developed in this work truncates the Fourier coeffic- 
ients of the raw data at this first minimum. The logic behind this approach 
is rather simple. Retaining fewer coefficients than this implies a smoothing 
of the deconvoluted spectrum and the loss of potentially valuable structural 
information. Adding more coefficients gains nothing because the signal- 
to-noise ratio 1s of order one at the minimum. In addition, truncating the 
coefficients at the point where the white noise begins to dominate forces the 
dispersion of the solution to remain finite. 

The truncated Fourier transform of the raw data 1s then used to cal- 
culate, a reasonable .estimate, of .the. desired, .function..™ it- was..s.howjiUn..Chapter- . 
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III that this deconvolution algorithm does indeed yield reasonable results. 
It is well known that truncating the coefficients in the manner described 
results in the introduction of a filtering error. The consequences of the 
introduction of this error source were discussed in detail in Chapter IV and 
will be reviewed below. 

In order to test the utility of the FFT/PS deconvolution algorithm a 
series of high purity, large grain size aluminum samples, containing known de- 
fect types and concentrations, were prepared. The thermal treatments employed 
to prepare the samples containing various concentrations of monovacancy or 
dislocation type defects were described in detail in Chapter II. Basically 
the procedure consisted of a two stage quench from elevated temperatures fol- 
lowed by a low temperature anneal and finally a high temperature anneal. The 
DBPAS spectrum was measured after each of these three thermal treatments. It 
was assumed, without transmission electron microscopy but based on arguments 
presented in Chapter II, that each of the three treatments introduced a single 
dominant defect type. These defect types are: 1) vacancies, after the 
quench, 2) prismatic dislocation loops, after the low temperature anneal, and 
3) the bulk or "defect free" material, after the high temperature anneal. 

In Chapter IV the lineshape parameters of the DBPAS measurements were 
calculated (for both the raw and the deconvoluted data) as a function of 
quenching temperature, A comparison of these lineshape parameters revealed 
several interesting trends. First and foremost, it appears that the FFT/PS 
deconvolution method is capable of automatically correcting for accidental 
and/or intentional changes in the instrument energy resolution function (IRF). 
It must be mentioned that there are some restrictions on the types of changes 
in the IRF which deconvolution can correct and these restrictions are dis- 
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data is that results obtained in different laboratories should be more direct- 
ly comparable if the instrumental broadening effects have been accurately re- 
moved. This ability to directly compare numerical values for standard line- 
shape parameters should help to clear up some of the discrepancies which 
currently exist in the literature. 

While it is fairly clear that accurate deconvolution should correct 
for changes in the IRF and improve the agreement of results obtained using 
different instrumentation, some of the other trends observed from the compar- 
ison of lineshape parameter calculated before and after deconvolution were not 
entirely predictable. For example, based on the limited results presented in 
this work, both the sensitivity and the accuracy of the lineshape parameters 
are improved by deconvolution. One explanation for the improved accuracy, 
i.e. the decrease in the error ratio, of the deconvoluted lineshape parameters 
is that the filtering operation has essentially removed the statistical count- 
ing errors present in the raw data. A comparison of multiple DBPAS measure- 
ments performed on the same sample show that before deconvolution, differences 
in the statistical counting error 1n the peak region can lead to a spread in 
calculated x(l) values. However, after deconvolution, the energy distribu- 
tions show less variation from experiment to experiments, and hence yield x(l) 
values with correspondingly less scatter. 

The increased sensitivity of the deconvoluted parameters may be 
related to the removal of the "broadening" action of the experimental instru- 
mentation. Recall that the major effect of the instrumentation 1s to spread 
counts from one channel into several adjacent channels. Therefore, the true 
lineshape changes in the annihilation energy distribution, which are caused 
by the trapping of positrons at defects, are spread out over a much larger 
JUSfefiT- . of . . cjhan.nfls.,1 n Jth^^aw jLala.^- . The- deconvolution procedure -r^r^r*^-f n r - 
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this effect and, in so doing, restores some of the details of the "true" 
annihilation energy distribution. This results in an Increase in observed 
defect sensitivity. 

Even with the improvements in lineshape sensitivity and accuracy due 
to the FFT/PS algorithm, the DBPAS technique must have the additional ability 
to distinguish between multiple defect types If it is to become a useful tool 
in the study of industrial materials. In order to investigate the usefulness 
of deconvolution in resolving multiple defect types, a comparison was made 
between the values of the Rl parameters calculated before and after deconvolu- 
tion. If the number of reconstruction Fourier coefficients is allowed to vary 
between samples, as would be the case if the FFT/PS method is used in its 
standard form, then the calculated Rj parameter value yields little useful 
information. 

In order to explain this somewhat negative result, the consequences of 
truncating the Fourier coefficient, also known as low pass filtering, were ex- 
amined in greater detail. It is known that the truncation of the coefficients 
introduces a filtering error with an amplitude Inversely proportional to the 
number of retained coefficients, w c , and a frequency which is proportional to 
w c . Ideally, one would like to greatly increase the number of retained coef- 
ficients in order to reduce the amplitude of this filtering error. Unfortu- 
nately, as w c is increased beyond the first minimum of the power spectrum the 
calculated solution will contain large, rapidly oscillating components which 
are not physically plausible. The mathematics of this argument were presented 
in detail in Chapter III. 

Since it is not possible to reduce the amplitude of the filtering 
error by arbitrarily increasing w c an alternative method was developed to 
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the procedure 1s to calculate a w c value for each of the data sets to be 
analyzed. Next the minimum value of w c is selected and subsequently employed 
in the deconvolution of all the data sets. Using this procedure does not 
eliminate the filtering error, but rather controls the source such that its 
contribution to the calculated Rj values remains approximately constant from 
experiment to experiment. Although this solution 1s not Ideal, the data 
presented In Chapter IV clearly Indicate that by controlling the error source 
the reproducibility of the Rj parameter, and its usefulness in the separation 
of multiple defect types, is improved. 

To summarize this section it can be stated that the FFT/PS deconvolu- 
tion algorithm appears to improve the quality of the data as obtained from 
UBPAS measurements in several ways. These improvements, as presented in this 
work include: 1) automatic correction for accidental and/or intentional 
changes in the instrument resolution function, 2) the potential of direct 
comparison of quantitative results between laboratories, 3) an improvement in 
both lineshape sensitivity and accuracy, and 4) the ability to distinguish 
between vacancy and dislocation type defects, as demonstrated for large 
grained high purity aluminum. 

To generalize these conclusions, further studies of this kind should 
be performed using a variety of different alloys with known defect type and 
concentration. In addition, it would be extremely useful to repeat the entire 
series of experiments using a DBPAS system which is both digitally and temper- 
ature stabilized. Such a data set would allow one to separate the effects of 
instrument broadening, which are theoretically correctable by deconvolution, 
from the nonsystematic effects of instrumental drift. In addition, spectra 
collected on a stabilized system could conceivably possess better counting 
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increase the number of coefficients retained in the deconvolution procedure. 

Additional experiments should be performed on the FFT/PS algorithm. 
These experiments could be designed to investigate alternative filtering 
schemes with the hope of reducing the associated filtering error. Most 
notably, the use of cosine tapering appears to be a viable approach. Once the 
filtering error has been reduced or eliminated the utility of moment analysis 
in the separation of multiple defect types should be reinvestigated. It Is 
not known at this time if the current FFT/PS algorithm is capable of disting- 
uishing the formation of positronium in a nonmetallic material. Positronium 
would show Itself as a relatively high frequency characteristic of the annihi- 
lation energy distribution. The key factor is determining whether the fre- 
quencies required to observe positronium are to the left of the first minimum 
in the power spectrum. This question could be investigated using both numeri- 
cally simulated data and data obtained from materials in which positronium 1s 
known to form. 

In order for Ooppler broadened positron annihilation spectroscopy to 
become a widely used industrial NDT tool it must continue to evolve in several 
ways. Its energy resolution must improve, the ability to separate multiple 
defect types must be further developed, and the reproducibility and reliabil- 
ity of the experimental results must increase. Although the fast Fourier 
transform/power spectrum analysis deconvolution technique has clearly not 
solved all of these problems, the results presented in this work indicate that 
this method measurably improves the quality and reproducibility of the infor- 
mation obtainable from DBPAS measurements, and serves as a sound basis for 
future refinements. 
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APPENDIX A 



In this appendix the equations used to theoretically determine the 
stopping, or thermal izati on, times of positrons in metals will be discussed. 
For a more complete treatment of this subject see Brandt and Arista (1982). As 
briefly mentioned in the text, positrons slow down mainly by two mechanisms: 
1) ionization and excitation of electrons which is characterized by the elec- 
tronic stopping power, S e , of the material, and 2) phonon interactions charac- 
terized by the phonon stopping power, S p n* The second mechanism dominates at 
epithermal energies below some critical value E c characteristic of the mater- 
ial. In the range from the maximum source energy, E s , to E c the electronic 
stopping power dominates. This energy range can be further divided into two 
subranges separated by the photon energy for which the material has its maximum 
stopping power, E m . In the range ^ < E < E s an accurate estimate for the 
stopping times is given by 



t] = 7.7p-l(E i /mc 2 )l- 2 



(A.l) 



where t] is the stopping time in psec 

p is the density of the material in g/cm 3 
mc 2 is O.bll MeV 



and 



£f is the initial positron energy 



since 



Ej/mc 2 < E s /mc 2 * 1 



(A.2) 



we can rewrite Eq. (A.l) as 



105 

*' ' 7 - 7,H (A.3) 

whicli for most TO tals of Merest results In st 0 ppt„ 3 tt«s of - , p Sec for the 
high energy region. 

For E c < E < th e electron gas theory of conduction electrons with 
Fermi energy E F is used to produce the equation 

S = JL. m2 f2i3/2 E 5/2 

e 105 h l?J (A.4) 
in the case of a continuous slow down approximation 

t = - iHl 1 / 2 f £f d£ 
s l 2 J J E i "pTO (A.5) 

Substituting Eq. (A.4) into Eq. (A.5) and integrating yields 

t . 105 h£ F 
2 " WJ~2 (A. 6) 

fc f 

where t 2 is in picoseconds. 
In the low energy region 3/2 KT < E < E c phonon interactions dominate with 

.V 

where w is the lattice deformation potential. 
Eq. (A.4) can be combined with Eq. (A.7) to define E c yielding 

E 3/2 _ 10Sb 2 m 5/2 E F 3 

c "71727-^-— (A.8) 

Inserting Eq. (A.7) i„ t0 Eq. (A . 5 ) and Integrating from 3/2 KT to E c as defined 
by Eq. (A.8) yields 

t = *h 4 p , /- 3K B T -,l/2i 

3 3^7^m 1 ' 'i (A.9) 



S = 2 m w F 
P" TUT* (A.7) 
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usiny typical values for common metals 1t can be shown that 

t3 > *2 > *1 (A.10) 
Thus, phonon thermal izati on times are larger than electronic stopping times. 
The total stopping time is given by 

t s = t } + t 2 + t 3 (A.ll) 
with typical theoretical values of 13 psec in A* and 45 psec in Pb (Brandt and 
Arista 1982). 



APPENDIX B 



The purpose of this appendix is to describe the relationships between 
the momentum distribution of the electrons in the sample and the measured 
Doppler shifts in energy as well as the angular deviations from colinearity of 
the Y-rays. In the center of mass rest reference frame the annihilation Y-rays 
will be emitted in directions separated by 180° and each with an energy of 611 
KeV which is equal to the rest mass of an electron (or positron). However, in 
the laboratory frame the momentum of the electron (the positron is assumed to 
be at rest as described in the text) produces changes in both the observed y- 
ray energy and the angle between the Y-rays. The geometry of the situation is 
shown in the diagram below. 

^ Center of Mass 

\ Velocity V 

\ 
\ 

\ 




Straightforward application of ttfe Lorentz transformations of special 
relativity combined with some trigonometric manipulation yield the result that 

2 

«4> - (%) sin * o + e (JL.) (8.1) 

c 

where is the angular deviation from 180° 

v is the electron velocity. 

(107) 
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c is the speed of light 

♦ 0 is the angle between the electron velocity vector and the 
direction of Y-ray emission 

2. .2 

2 wiih.ii is negngiDle for v « c 



and e(^.) is a term of order ±- wn1ch 1s negligible for 



(For a detailed derivation of this result see Rindler 1960). The relativistic 
definition of momentum 1s 

p = am 0 V (B.2) 

where p is the momentum 

„=(i-4-r 1/2 

and m 0 is the particle mass 
For v « c, o ♦ 1 and p - mv 
Substituting this result into Eq. (B.l) yields 

6 * " s1n *o (B.3) 

Defining the transverse momentum component as 

p T - p sin * 0 (B>4) 
and substituting Eq. (B.4) into Eq. (B.3) yields the final result 

2P, 



6$ - — I 



(B.5) 



demonstrates that measuring the distribution of the angular deviations, ««, is 
in essence equivalent to measuring the distribution of the momentum of the 
electrons within the sample. 

In the case of Doppler measurements one has 

6P Y = P e cos ♦„ (B _ 6) 
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where 6P Y Is the change in momentum of the Y-ray 

P e is the electron momentum 
and 4> 0 is as defined above 
Using Eq. (B.2) with v « c Eq. (B.6) takes the form 

6P y = nfeVg cos * 0 (B . 7) 
It can be shown that in relativistic dynamics 

E 2 = P 2 c 2 + moc2 (B>8) 
where E is the relativistic energy 
P is the relativistic momentum 
c is the speed of light 
and m 0 is the particle rest mass 

Now for a photon (Y-ray) m,, = 0 so Eq. (B.8) reduces to 

E Y * Py c (B.9) 
Combining Eqs. (B.6), (B.7), and (B.9) yields 

6E Y * 6p T c = p e (cos <. 0 )c - meCv e (cos * 0 ) (B.10) 
where 6E y is the Doppler energy shift of the observed Y-ray 
6p T is the change in momentum of the Y-ray 
and p e is the momentum of the annihilating electron 

From Eq. (B.1U) it can be seen that the momentum distribution of the annihi- 
lating electrons can be deduced from the measured Doppler energy shifts. 

Combining Eq. (B.5) with (B.10) yields the following relationship 
between 6$ and 6E: 

6E = 255.56$ (g ] j j 

where 5E is in eV 
and 6$ is in mrad 
Finally, by combining Eq. (B.10) and the fact that for an electron m 0 c2 = 6 n 
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KeV it can be shown that 

v 

6E T = 511 (cos 4 0 ){-f) (B.I 2) 

where 6E Y is in KeV 

For the case of a 1 eV electron (v e = 5.93 x 10? cm/sec) moving in the direc- 
tion of Y-ray emission, Eq. B.12 predicts a 1.01 KeV shift in observed energy. 



APPENDIX C 



The material examined in this work was high purity (99.9995%) alumi 
supplied by the Materials Research Corporation. The chemical analysis, in 
terms of PPM of metallic interstitials (and some gases) is listed below. 



t lenient 


PPM 


Element 


PPM 


Element 


PPM 


An 

Ag 


<U.l 


L1 


<0.1 


U 




A 1 


major 


Mg 


1.5 


V 


<0.1 


A c 
MS 


<U.l 


Mn 


<0.1 


W 


<0.2 


An 


<u.l 


Mo 


<0,1 


Y 




R 
D 




Na 


0.2 


Zn 


0.4 


Ba 




Nb 


<0.1 


Zr 


<0.1 


Be 




Ni 


0.1 






Bi 


<0.1 


Pb 


<0.1 


CI 


0.6 


Ca 


<0.1 


Pd 


<0.1 


F 




Cd 


<0.1 


Pt 


<0.14 


H 


3 


Co 


<0.1 


Kh 


<0.1 


N 


5 


Cr 


<0.1 


Sb 


<0.1 


0 


20 


Cu 


0.1 


Si 


0.9 


P 


0.3 


Fc 


<0.1 


Sn 


<0.3 


S 




Ga 


<0.1 


Ta 








Ge 


<0.1 


Th 








In 


<0.1 


Ti 


<0.1 






K 


<0.1 


Tl 









APPENDIX D 



In this Appendix the utility of Fourier analysis in the solution of the 

convolution integral equation will be demonstrated. In the first section of 
Chapter III it was shown that 

h(y) = / f(x)r(y - x)dx (D #1 ) 
By definition of the Fourier transform 

F(w) = f e1wx f(x)dx (D>2) 



Therefore 



F(w)K(w) = [/ e1wx f(x)dx]r / e^rUJdt] (D.3) 

OO OD 

" L L eiWXf t x ) e1Wt ''(t)dxdt (D.4) 

= / / eiw(xn) f(x)r(t)dxdt (D 5) 



Now let x + t = y 
so t = y - x 
and dt = dy 
Hence Eq. (D.b) can be rewritten as 



F(w)R(w) = eiwy f(x)r(y . x)dxdy (D 6) 

00 GO 

■ / elwyr/ f(x)r ( y . x)<Jx3dy (D 7) 
(112) 
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But the Integral within the brackets is by Eq. (D.l) equal to h(y). Therefore 
Eq. (U.7) can be rewritten as 

GO 

F(w)R(w) = / e™y n (y)dy ( 0 .8) 

and using Eq. (0.2) on the right hand side 

F(w)R(w) = H(w) (0.9) 
This demonstrates that in Fourier space convolution is simply a complex multi- 
plication. 



APPENDIX E 



The purpose of this Appendix is to demonstrate some of the mathematical 
difficulties associated with the incorrectly posed convolution integral 
equation. The primary source of the problem is that the left hand side of the 
equation 

* b 

Mx) = / f(y)r(x - y)dy (E.i) 

d 

is only known to within some error e. That is, if h(x) is the "true" function 
then 

max |h*(x) - h(x)| < e (E.2) 
a<x<b 

Any function f}(y) which when convoluted with the known resolution function 

* 

r(x) yields a function hj(x) that satisfies Eq. (E.2) must be considered a 
solution of Eq. (E.I). 

Assume there are two solutions to Eq. (E.I) which are f](x) and f 2 (x) = 
f](x) + c sin(2unw/(b-a). By definition 

* b 

Mx) = / f](y)r(x - y)dy (E.3) 
a 

and 

* b 

M x ) » / [f^y) + C sin (2imy/(b - a))]r(x - y)dy (E.4) 
a 

Subtracting Eq. (E.3) from Eq. (E.4) yields 
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* * b 

Mx) - Mx) = / r ( x " y)C sin(2irny/(b - a))dy (E.5) 
a 

Now integrating by parts with 

u = r(x - y) dv = C sin (|^)dy 

yields 

/b r(x . y)c ...(gflL^ = r(x . rifHjjpl cos 



Now integrating the last term on the right hand side by parts with 

dy v = r(x - y) 

Yields 

/J »s(|l^) *» dy = [r(x - y) cos (gflL)]* 

- £ - - y)(|^) st„ (**-, dy (E.7) 

Combining Eqs. (E.6) and (E.7) yields 

/J '<* " y) C sin (|^) dx = - C ^ N - a >l> (x - y) cos (|^)] a b (E.8) 

and combining Eqs. (E.8) and (E.5) yields the result 

l£(x) - h*(x) = £ U±£±) r(x - y) cos (E.9) 



So it can be seen that for any value of C one can select a value of N large 
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enough so that the difference between h2(x) and hj(x) is smaller than any 
previously determined number. Therefore, both hi(x) and h 2 (x) satisfy Eq. 
(£•2) and f](x) as well as f 2 (x) must be solutions to Eq. (E.l). The important 
point is that one may select C and N in such a way as to allow f](x) and f 2 (x) 
to differ from each other by as much as one likes. Clearly the problem as 
posed does not have a unique solution and a different solution technique is 
required. 



APPENDIX F 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C Q 

C THIS IS A LISTING OF THE FORTRAN 77 VERSION OF THE PROGRAM C 

C DECONFFT WHICH IS DESIGNED TO DEOONVOLUTE DOPPLER BROADENED C 

C POSITRON ANNIHILATION SPECTRA. FOR A DESCRIPTION OF THE C 

C MATHEMATICS INVOLVED SEE J.P. SCHAFFER, E.J. SHADGHNESSY, C 

C AND P.L. JONES, JOURNAL OF NUCLEAR INSTRDMENTS AND METHODS C 

C IN PHYSICS RESEARCH, B5, 74, 1984. - 

C C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

INTEGER PCSP 
CHARACTER*2 PCSZ,PCH 
CHARACTER* 80 XCB , YCH , TCE 
CHARACTER* 7 NAME, NAME E 

DIMENSION IA(2048) , IB(2048) , D0PE(12) ,XP(74) ,YP(74) 

DOUBLE PRECISION A(2048) ,B (2048) ,C(2048) ,FR(2048) ,HC (2048) 

DOUBLE PRECISION D(2048) ,E(2048) ,G(2) ,S (2048) , SI (2048) 

DOUBLE PRECISION Fl (2048) , PI, V,W,W1 

COMPLEX*16 R(4096),H(4096) ,F(4096) ,T,A1,H1(4096) 

COMPLEX*16 A2,R1(4096),T1 

PI-3.141592653589793D0 

PRINT*,' THIS PROGRAM DEALS WITH THE CONVOLUTION 1 
PRINT*, 'B(X) a F(X) *R(X) , WHERE R IS THE* 
PRINT* , • RESOLUTION FUNCTION, H IS THE ' 
PRINT*, 'BROADENED AND OBSERVED FUNCTION, AND' 
PRINT*,' F IS THE TRUE AND DESIRED PUNCTION. ' 
PRINT*,' ' 

C THE RESOLUTION FUNCTION IS USUALLY OBTAINED FROM 

C THE Sr(85) GAMMA RAY (514 KeV) . 

C 

2 PRINT*, 'THE PROGRAM REQUIRES N DATA POINTS' 
PRINT*, 'WHERE N=2**M, AND M IS AN INTEGER. * 
PRINT*, * ' 

C THESE DATA POINTS WILL BE H(X) AND R(X) AS COLLECTED 
C BY THE MCA. 



C 



PRINT*, 'ENTER THE VALUE OF M TO BE USED IN THIS' 

PRINT*, 'PROGRAM. (M=ll YIELDS N=2048) ' 

PRINT*, ' ' 

READ(*,' (12) ') M 

N=2**M 

ND2=N/2 

ND4-N/4 

PRINT* , ' ENTER THE RESOLUTION FUNCTION FILENAME' 

PRINT*,' • 

READ (*, ' (A7) ') NAME 

OPEN ^OfFILEsHAMEfPAD-'YES') 

READ (90,1) (IA(I), 1=1,2048) 

FORMAT^.(5X*5X7,) - -- : ■- - 

CLOSE (90) 

(117) 
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PRINT*, •ENTER THE EXPERIMENTAL FUNCTION FILENAME* 

PRINT*, • » 

READ(*,'(A7) •) NAME 

OPEN (90,FILE=NAME,PAD= ' YES*) 

READ (90,1) (IB(I), 1=1,2048) 

CLOSE (90) 

NIN-2048 

C THE SUBROUTINE DELSKIP DELETES THE LINE 
C SKIP IN THE RAW DATA SUPPLIED BY THE MCA. 

CALL DELSKIP (IA,A,NIN) 

CALL DELSKIP (IB,B,NIN) 
C THE FOLLOWING SUBROUTINE CALCULATES THE CENTRAL 
C CHANNEL LOCATION BASED ON THE CENTER OF THE 
C 11 CHANNEL REGION WHICH CONTAINS THE MAXIMUM 
C TOTAL NUMBER OF COUNTS 

CALL CENTER (A,IMAXA) 

CALL CENTER (B,IMAXB) 

C 

C THIS NEXT LARGE SECTION DEALS WITH THE 

C BACKGROUND CORRECTION. 

C 

PRINT*,' DO YOU WISH TO BACKGROUND CORRECT? * 

READ (*, ' (Al) •) NAME 

IF (NAME .EQ. *Y') GO TO 3 

CALL NOBACK (N,B,IMAXB) 

GO TO 200 

3 PRINT*,' DO YOU WISH TO BACKGROUND CORRECT USING' 

PRINT*, 'THE JACKMAN METHOD OR THE SPECTRAL FILTERING' 

PRINT*,' TECHNIQUE? (J OR S) ' 

READ (*,'(A1) ') NAME 

IF (NAME .EQ. »J») GO TO 5 

PRINT*, 'THE NEXT SEVERAL QUESTIONS REFER TO THE* 
PRINT*, 'RESOLUTION FUNCTION. ' 
CALL BACK (A,IMAXA,N) 

PRINT*, 'THE NEXT SEVERAL QUESTIONS REFER TO THE' 
PRINT*, 'EXPERIMENTAL FUNCTION. • 
CALL BACK (B,IMAXB,N) 
GO TO 377 

5 PRINT* , ' ENTER THE RESOLUTION OF THE MCA IN' 
PRINT*,' eV/ CHANNEL. (USE A DECIMAL POINT)' 
PRINT*,' (SUGGESTION: IF THE SYSTEM IS SET' 
PRINT*,' IN ITS NORMAL MODE USE 64. )' 
PRINT*,' ' 
READ (*,' (F6.2) ') EV 
X= ( 1 00 0 . D0/EV) +0 . 5 DO 
IK=IFIX (X) 
IK LA- ( IKAXA-20) /IK 
ILA=514-IKLA 
IKHA= (1 02 4- 1 MAX A) /IK 
IHA=514+IKHA 
IKLB= ( IMAXB-20) /IK 
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ILB-511-IKLB 

IKHB= ( 1 0 24- IMAXB ) /IK 

IHB=511+IKHB 

PRINT*, 1 THE EXPERIMENTAL SPECTRUM HAS BEEN COLLECTED ' 
PRINT*, 'OVER THE FOLLOWING KeV REGION:' 
WRITE (*,10) ILB, IHB 
10 FORMAT (15,15} 

PRINT*,' WE ARE NOW READY TO PROCEED WITH THE' 

PRINT* , ' BACKGROUND CORRECTION. OVER WHICH' 

PRINT*, 'LOW ENERGY REGION WOULD YOU LIKE THE' 

PRINT*,' LINEAR LEAST SQUARE FIT? IF POSSIBLE' 

PRINT*, 'SELECT 486 TO 497. THE PROPER FORM* 

PRINT*, 'OF ENTRY IS TWO 3 DIGIT INTEGERS' 

PRINT*, 'SEPARATED BY TWO BLANKS.' 

PRINT*,* ' 

READ (*,*) J LAB, J LB B 

IS* IMAXB- (511-JLAB) *IK 

IF= IMAXB- (511-JLBB) *IK 

N1=IF-IS 

DO 20 1*1, Nl 

J=IS-1+I 

D(I)=J 

E(I)=B(J) 
20 CONTINUE 

CALL FIT (N1,D,E,G) 

BBL=G(1) 

XMBL=G (2) 

PRINT* , ' OVER WHICH HIGH ENERGY REGION WOULD' 

PRINT* , * YOU LIKE THE LEAST SQUARE FIT?' 

PRINT* , * IF POSSIBLE SELECT 534 542.' 

PRINT*, * ' 

READ (*,*) JBAB,JBBB 

IS«IMAXB+ ( JHAB-511) *IK 

IF=IMAXB+ ( JHBB-511) *IK 

N1«IF-IS 

DO 30 I«1,N1 

J=IS-1+I 

D(I)=J 

E(I)«B(J) 
30 CONTINUE 

CALL FIT (N1,D,E,G) 

BBH=G(1) 

XMBB=G(2) 

PRINT*, 'TO WHICH ENERGIES WOULD YOU LIKE TO' 
PRINT* , ' EXTRAPOLATE THE LINEAR BACKGROUND 1 
PRINT* , * CORRECTION? (SUGGESTION: 503 518)' 
PRINT* , * • 
READ (*,*) IL,IH 

PRINT* , ' ON THE HIGH ENERGY SIDE OF THE PEAK, * 
PRINT*, 'WHAT IS THE HIGHEST ENERGY YOU WISH' 
PRINT*,' TO HAVE CONSIDERED FOR DECONVOLUTION? ' 
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PRINT*,' SUGGESTION: 540* 
PRINT*,' ' 
READ (*,*) I END 
IL=IMAXB- (511-IL) *IK 
IH-IMAXB+ ( IH-51 1 ) * IK 
IEND=IMAXB+ ( IEND-511) *IK 
IEB=IEND 
DO 35 I=IH,IEND 

B(I)«B(I)-(XMBB*I+BBH) 

IF (B(I) .LT. 0.) B(I)«0.D0 
35 CONTINUE 

XH»XMBH * IH+BBH 

XL-XMBL* IL+BBL 

XM= ( XL-XH ) / ( FLOAT (IL-IH) ) 

XB=XH- ( (XL-XH)/ (FLOAT ( IL-IH) ) ) *IH 

DO 40 I=IMAXB,IH 

B(I)=B(I)-(XM*I+XB) 

IF (B(I) .LT. 0.) B(I)=0.D0 
40 CONTINUE 

PRINT* , ' THE RESOLUTION SPECTRUM HAS BEEN COLLECTED' 
PRINT*, 'OVER THE FOLLOWING KeV REGION:' 
WRITE (*,10) ILA, IHA 

PRINT* , ' OVER WHICH LOW ENERGY REGION WOULD YOU LIKE' 

PRINT*,' THE LINEAR LEAST SQUARE FIT?* 

PRINT*, 'IF POSSIBLE SELECT 489 500. ■ 

PRINT*,' ' 

READ (*,*) JLAA, JLBA 

IS=IMAXA- (514-JLAA) *IK 

IF=IHAXA-(514-JLBA) *IK 

N1=IF-IS 

DO 45 I«1,N1 

J=IS-1+I 

D(I)=J 

E(I)=A(J) 
45 CONTINUE 

CALL PIT (Nl,D,E f G) 

BAL S G(1) 

XMAL=G(2) 

PRINT* , ' OVER WHICH HIGH ENERGY REGION WOULD YOU LIKE' 

PRINT*,' THE LINEAR LEAST SQUARE FIT?' 

PRINT*, 'IF POSSIBLE SELECT 537 545.' 

PRINT*,' ' 

READ (*,*) JHAA, JHBA 

IS=IMAXA+(JHAA-514) *IK 

IF=IMAXA+( JHBA-514) *IK 

N1=IF-IS 

DO 50 1=1, Nl 

J=IS-1+I 

D(I)=J 

E(I)=A(J) 
50 CONTINUE 



CALL PIT (N1,D,E,G) 

BAH«G(1) 

XMAH=G(2) 

PRINT*, 'TO WHICH ENERGIES WOULD YOU LIKE TO' 
PRINT*,* EXTRAPOLATE THE LINEAR BACKGROUND ' 
PRINT* , ' CORRECTION? (SUGGESTION: 506 521)' 
PRINT*,* ' 
READ (*,*) IL,IH 

PRINT*,' ON THE HIGH ENERGY SIDE OP THE PEAK, • 

PRINT*, 'WHAT IS THE HIGHEST ENERGY YOU WISH' 

PRINT* , ' TO HAVE CONSIDERED FOR DECONVOLUTION»• 

PRINT*, 'SUGGESTION: 521' 

PRINT*,' ' 

READ (*,*) I END 

IL»IMA7A- (514-IL) *IK 

IH=IMAJA+(IH-514)*IK 

IEND-IMAXA+ ( IEND-514) *IK 

IEA=IEND 

DO 55 I=IH,IEND 

A ( I ) «A ( I) - (XMAH*I+BBH) 

IP (A(I) .LT. 0.) A(I)=0.D0 
55 CONTINUE 

XB=XMAH*IB+BAH 

XL=XMAL*IL+BAL 

XM=(XL-XH)/ (FLOAT (IL-IH) ) 

XB*XH- ( (XL-XH) / ( FLOAT ( IL- IH) ) ) *IH 

DO 60 I=IMAXA, 10 

A(I)=A(I)-(XM*I+XB) 

IF (A(I) .LT. 0.) A(I)=0.D0 
60 CONTINUE 
GO TO 378 

377 CONTINUE 
IEA»IMAXA+109 
IEBsIMAXB+453 

378 CONTINUE 

C 

C THIS IS THE END OF THE BACKGROUND CORRECTION SEQUENCE 
C NEXT WE NEED TO REFLECT THE HIGH ENERGY HALF OF THE 
C BACKGROUND CORRECTED PEAK ABOUT THE CENTRAL POINT 
C TO FORM A SYMMETRIC PEAK. IN ADDITION A ZEROFILL 
C TO 2048 DATA POINTS IS REQUIRED 
C 

CALL ZEROFILL (A, IMAXA, IEA,N) 

CALL ZEROFILL (B, IMAXB, IEB,N) 

M1=IEA-IMAXA+10 

CALL TAPER (A,M1,N) 

M2=IEB-IMAXB+20 

CALL TAPER (B,M2,N) 

M3=K1 

CALL SMOOTH (A,M3,N) 
M4=M2 



CALL SMOOTH (B,M4,N) 
C THE NEXT TWO 00 LOOPS NORMALIZE THE AREA 
C OP THE RESOLUTION AND BROADENED SPECTRA. 
V=O.DO 
W-O.DO 
DO 62 1-1, N 
V-V+AU) 
W=W+B(I) 

62 CO NT IN DE 
DO 63 I«1,N 

A(I)=A(I)/V 
B(I)=B(IJ/W 

63 CONTINUE 
DO 68 1=1, N 

S(I)=B(I) 
68 CONTINUE 

DO 72 I=1,ND2 

J-ND2+I 

V=S(I) 

S(I)=S(J) 

S(J) «V 
72 CONTINUE 

PRINT*,' DO YOU WISH TO SAVE THE BACKGROUND ' 
PRINT*, •CORRECTED SPECTRUM? 1 
READ (*,»(A1) •) NAME 
IP (NAME .EQ. 'N*) GO TO 76 
PRINT* , 'ENTER THE FILENAME. ' 
READ (*,«(A7) ») NAME 
OPEN (90, FILE=NAME , PAD= * YES * ) 
WRITE (90,74) (S(I),I=1,N) 
74 FORMAT (5X t 5D22.10) 
CLOSE (90) 
CALL PARAMS (N,S) 

76 CONTINUE 
DO 77 1=1, N 

R(I)=DCMPLX(A(I) ,0.D0) 
H(I)=DCMPLX(B(I) ,0.D0) 

77 CONTINUE 

CALL REORDER (M,H) 
CALL FFT (M,H) 
CALL REORDER (M,R) 
CALL FFT (M,R) 
CALL DECON (M,F,R,H) 

DOPE{1)=150. 

DOPE (2) =150. 

DOPE{3)»950. 

DOPE (4) =650. 

DOPE(5)=0. 

DOPE(6)=0. 

PCSZ='NS' 
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PCSP-1 

c 

C THIS NEXT SMALL SECTION ALLOWS ONE TO VIEW THE POWER 
C SPECTRUM IF SO DESIRED. 

PRINT*,' WODLD YOO LIKE TO SEE AND PLOT THE POWER SPECTRUM ' 

PRINT*, 1 OF THE EXPERIMENTAL DATA AND THE DECONVOLDTED DATA. 1 

PRINT* , • • 

READ (*,'(A1) ') NAME 

IF (NAME .EQ. 'N' ) GO TO 179 

DO 78 1=1, ND2 

D(I)«H(I)*(CONJG(H(I))) 
78 CONTINUE 

PMNT*, | ENTER THE FILENAME FOR THE P.S. OF THE EX. DATA. ' 

READ (*,'(A7) •) NAME 
OPEN OCFILE^NAMEfPAD^'YES 1 ) 
WRITE (90,901) (I,D(I),I=1,ND2) 
901 FORMAT (5(14,022.12)) 
CLOSE (90) 
CALL REORDER (M,F) 
CALL FFT (M,F) 
T=0.DO 

DO 277 1=1, N 
T=T+F(I) 

277 CONTINUE 

DO 278 1=1, N 
F(I)=F(I)/T 

278 CONTINUE 

CALL REORDER (M,F) 
CALL FFT (M,F) 
DO 178 I=1,ND2 

E(I)*=F(I)*(CONJG(F(I))) 
178 CONTINUE 

PRINT* , * ENTER THE FILENAME FOR THE P.S. OF THE DECON. DATA • 

PRINT* , ' • 

READ (*,'(A7)«) NAME 

OPEN ( 90 , FILE=NAME, PAD= 1 YES ' ) 

WRITE (9 C, 901) (I,E(I) ,I=1,ND2) 

CLOSE (90) 

NPTS«74 

DO 910 X-1,10 

XP ( I) =ALOG10 ( FLOAT ( I) ) 

YP(I)=B(I) 

YP(I)=ALOG10(YP(I)) 
910 CONTINUE 

DO 912 1=15,100,5 

J=8+(I/5) 

XP ( J) =ALOGl 0 ( FLOAT (I) ) 
YP(J)=E(I) 
YP(J)=ALOG10(YP(J)) 
912 CONTINUE 
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DO 914 1=120,1000,20 
J=23+(I/20) 

XP ( J) =ALOGl 0 ( FLOAT ( I ) ) 
YP(J)=E(I) 
YP( J) =ALOG10 (YP( J) ) 
914 CONTINUE 

XP(74)«ALOG10(1020.) 
YP(74)=-10. 

XCH='LOG OP THE FOURIER COEFFICIENT NO • 
MX=34 

YCH='LOG OF THE POWER' 
MY=16 

PRINT*, •ENTER THE TITLE OF THIS PLOT * 

PRINT*, • ' 

READ (*,' (A80) ') TCH 

PRINT*, 'ENTER THE NO. OF CHARACTERS IN THE TITLE • 
READ (*,'{I2)') NT 
PCB= , 1 ' 

CALL CKSTUP ( 'TK15 • ,9600) 
CALL CKER 

CALL CKGRAF (DOPE,NPTS,XP, YP,XCH,MX, YCH,MY,TCH,NT, PCSZ 
C 9 PCSP r PCH) 
DO 930 1=1,10 
yP(I)=D(I) 
YP ( I ) ■ ALOG1 0 ( YP ( I ) ) 
930 CONTINUE 

DO 932 1=15,100,5 
J=8+(I/5) 
YP(J)=D(I) 
YP ( J) =ALOG10 ( YP ( J) ) 
932 CONTINUE 

DO 934 1=120,1020,20 
J=23+(I/20) 
YP{J)=D(I) 
YP ( J) =ALOGl 0 ( YP ( J) ) 
934 CONTINUE 
PCH='3 ' 

CALL CKPLT ( DOPE , NPTS, XP, YP.PCSZ r PCSP, PCH) 
CALL CKMA (0.,0.,1) 
READ (*,*) 
179 CONTINUE 

C 

C AT THIS STAGE THE RESOLUTION AND EXPERIMENTAL 
C SPECTRA HAVE BEEN TRANSFORMED INTO FOURIER SPACE. 
C THIS NEXT SMALL SECTION DETERMINES THE NUMBER OF 
C POSITIVE DECREASING FOURIER COEFFICIENTS IN THE 
C DECONVOLDTED RESULT. THIS INFORMATION CAN THEN 
C BE USED IN THE FILTERING ROUTINES IF SO DESIRED. 
C 

DO 80 1=2, N 
J=I-1 
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D(I)=REAL(P(I) ) 

D (J) 'REAL (P (J)) 

IF (D(I) .LT. 0.0) GO TO 81 

IF (D(I) .GT. D(J)) GO TO 81 

80 CONTINUE 

81 IPD»I-1 

WRITE (*,82) IPD 

82 PORMAT (SX,'THE NO. OF POSITIVE DECREASING CO EPS. IS ' 15) 
PRINT*,' • 

PRINT* , * DO YOO WISH TO USE A COSINE TAPER FILTER' 
PRINT*, 'TO MODIFY THE FOURIER COEFFICIENTS? 1 
READ (*, ' (Al) •) NAME 
IF (NAME ,EQ. *N') GO TO 189 

PRINT*, 'ENTER THE NUMBER OF THE FIRST COEFFICIENT* 
PRINT*, 'TO BE SET EQUAL TO ZERO. * 
READ (*,'(I3) ') IZ 

PRINT*, 'ENTER THE NUMBER OF COEFFICIENTS TO BE TAPERED. » 
READ (*,'(I3) ') IT 
IS=IZ-IT 
XS-REAL ( IS) 
XZsREAL(IZ) 
XY=1.0/ (XS-XZ) 
DO 181 I=IS,IZ 
X=REAL(I) 

F(I)=P(I)*0.5*(1-COS(PI*XY*(X-XZ))) 
181 CONTINUE 
J=IZ+1 
K=(N/2)+l 
DO 183 W,K 

F(I)=DCMPLX (O.DO,0.DO) 
183 CONTINUE 
K=K-1 

DO 185 1-2, K 
J-N+2-I 
F(J) =F(I) 
185 CONTINUE 

189 PRINT* , 'DO YOU WISH TO USE A BOXCAR FILTER TO' 
PRINT*, 'MODIFY THE FOURIER COEFFICIENTS IN' 
PRINT*, 'ORDER TO REMOVE THE HIGH FREQ. NOISE? ' 
READ (*,' (Al) ') NAME 
IF (NAME .EQ. «N') GO TO 84 

PRINT*, 'ENTER THE NUMBER OF COEFFICIENTS TO BE USED: • 

READ (*,*) J 

J=J+1 

K=N+2-J 

DO 83 I=J,K 

F(I)=DCMPLX(0.D0,0.D0) 

83 CONTINUE 

84 CONTINUE 

PRINT*, 'DO YOU WANT TO SEE THE FOURIER COEFFICIENTS? ' 
READ (*,' (Al) ') NAME 



IP (NAME .EQ. 'N') GO TO 85 
PRINT*, 'ENTER THE FILENAME FOR B(X). * 
PRINT* , ' ■ 
READ (*,'(A7) ') NAME 
OPEN (90,FILE»NAME,PAD= ' YES' ) 
WRITE (90,911) (I,B(I), 1=1, N) 
911 FORMAT (2(I4,2D22. 10)) 
CLOSE (90) 

OPEN OCFILE-'FCSF'jPADe'YES') 
WRITE (90,911) (I,P(I), 1=1, N) 
CLOSE (90) 

85 CONTINUE 

CALL REORDER (M,F) 

CALL FFT (M,F) 

CALL REORDER (M,H) 

CALL FFT (M,B) 
C THIS NEXT SECTION ALLOWS ONE TO OSE THE ITERATION 
C ROUTINE TO IMPROVE F(X). 

PRINT* , 1 WOULD YOU LIKE TO USE THE ITERATION' 

PRINT*, 'SUBROUTINE TO TRY TO IMPROVE THE' 

PRINT*, 'ACCURACY OP THE DECONVOLUTED SPECTRUM? 

READ (*,» (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 86 

CALL ITER (M,R,F,H) 

86 CONTINUE 

C THIS NEXT DO LOOP SHIFTS THE PHASE OF THE 
C EXPERIMENTAL FUNCTION B, AND THE DECONVOLUED 
C FUNCTION F, SO THAT THEIR MAXIMA ARE LOCATED 
C IN THE CENTRAL CHANNELS. 
ND2=N/2 
DO 88 1=1, ND2 
J=ND2+I 
T=F(I) 
T1«B(I) 
F(I)«F(J) 
H(I)=B(J) 
F(J)«T 
H(J)=T1 
88 CONTINUE 
DO 98 1=1, N 

F1(I)=REAL(F(I) ) 

98 CONTINUE 
W=0.D0 

DO 99 1=1, N 
W=W+F1(I) 

99 CONTINUE 

DO 100 1=1, N 

F1(I)=F1(I)/W 
100 CONTINUE 

PRINT*,'DO YOU WANT A PLOT OF h(X) AND f(X)? ' 
READ (*, ' (Al) ') NAME 
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IP (NAME .EQ. 'N*) GO TO 105 
NPTS-60 
W«F1(1025) 
DO 102 1=1,60 
J-1024+I 

XP(I)=.064*(FLOAT(I)) 
YP(I)=F1(J)/W 

102 CONTINUE 
XCH^KeV* 
MX=3 

YCB=»h(X) AND £(X) • 
MY=13 

PRINT* , * ENTER THE TITLE OP THIS PLOT. ' 

PRINT* ' ' 

READ (*,'(A80)») TCH 

PRINT* , 1 ENTER THE NO. OP CHARACTERS IN THE TITLE. » 
READ (*, , (I2)') NT 
PCH-'l ' 

CALL CKSTDP ( 'TK15 ' ,9600) 
CALL CKER 

CALL CKGRAP (DOPE,NPTS,XP, YP r XCH,MX, YCH,MY,TCH,NT,PCSZ 
C,PCSP,PCH) 
DO 103 1=1,60 

J=1024+I 

YP(I)=S(J)/W 

103 CONTINUE 
PCH= , 3 * 

CALL CKPLT ( DOPE , NPTS , XP , YP , PCS 2 , PCSP , PCH ) 

CALL CKMA (0.,0.,1) 

READ (*,*) 
105 PRINT*, 1 DO YOO WISH TO SAVE THE • 

PRINT* , 1 DECONVOLUTED SPECTRUM? ' 

READ (*,*(A1) ') NAME 

IF (NAME .EQ. •N') GO TO 101 

PRINT* , • ENTER THE FILENAME. ' 

READ (*,MA7)') NAME 

OPEN (90 , FILE=NAME, PAD= 1 YES" ) 

WRITE (90,74) (Fl (I) ,1=1 ,N) 

CLOSE (90) 
101 CONTINUE 

PRINT*, 1 DO YOU WISH TO CHECK YOUR RESULT BY' 

PRINT* , * CONVOLUTING f(X) WITH r (X) AND COMPARING ' 

PRINT* , 'THE RESULT WITH h(X)?' 

READ (*,' (Al) ') NAME 

IF (NAME .EQ. »N') GO TO 200 

CALL REORDER (M,F) 

CALL FFT (M,F) 

CALL CONV (M,F,R, HI) 

CALL REORDER (M,H1) 

CALL FFT (M,H1) 

DO 106 1=1, N 
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BC(I)»REAL(H1(I)) 

106 CONTINDE 
V«0.D0 

DO 107 1=1 ,N 
V=V+HC(I) 

107 CONTINUE 

DO 108 I«1,N 

BC(I)«BC(I)/V 

108 CONTINDE 

PRINT*,' ENTER THE FILENAME POR THE RESULT ' 

PRINT*,' OF THIS COMPARISON. 1 

PRINT*,' ' 

READ (*,' (A7) ') NAME 

OPEN (90, FILE=NAME, PAD= 1 YES ' ) 

WRITE (90,110) 

110 FORMAT (3X, , I , ,14X,'h(X)RAW',14X, , h(X)=(F*R) »,14X,'f (X) ') 
WRITE (90,120) (I, S(I) ,HC(I) ,F1 (I) ,1=1025,1300) 

120 FORMAT (1X,I4,2X,3D22.10) 
CLOSE (90) 

PRINT*,' DO YOO WANT A PLOT OF h(X) AND THE ESTIMATE' 
PRINT*, 'OF h(X) OBTAINED FROM f (X)*r(X)? • 
READ (*,'(A1)') NAME 
IF (NAME .EQ. 'N') GO TO 200 
NPTS=60 
W=HC{1025) 
DO 130 1=1,60 
J=1024+I 

XP(I)=.064*(FLOAT(I)) 

YP(I)=HC(J)/W 
130 CONTINUE 
XCB='KeV 
MX=3 

YCE='h(X) AND h*(X) ' 
MY=14 

PRINT*,' ENTER THE TITLE OF THIS PLOT.' 

PRINT*, ' * 

READ (*, ' (A80) ') TCH 

PRINT*,' ENTER THE NO. OF CHARACTERS IN THE TITLE. • 
READ (*, ' (12) ') NT 
PCB«'l ' 

CALL CKSTOP ( * TK15 ' ,9600; 
CALL CKER 

CALL CKGRAF (DOPE,NPTS,XP, YP,XCH,MX,YCH,MY,TCB,NT,PCSZ 
C,PCSP,PCH) 
DO 140 1=1,60 
J=1024+I 
YP(I)=S(J)/W 
140 CONTINUE 
PCH= ' 3 ' 

CALL CKPLT (DOPE,NPTS,XP, YP, PCSZ , PCSP, PCH) 
CALL CKMA (0. ,0. ,1) 
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READ (*,*) 
200 CONTINUE 

CALL P ARAMS (N,F1) 

STOP 

END 

cccccccccccccccccccccccccccccccccccccccccccccccc 
c 

SUBROUTINE DELSKIP (IC,C,NIN) 

DIMENSION IC(2048) 

DOUBLE PRECISION C(1024) 

J=0 

K=0 

DO 10 I«1,NIN 
J=J+1 
K=K+1 

IC(I)=IC(K) 

IP (J .EQ. 5) THEN 

J=0 

K»=K+5 

END IF 
10 CONTINUE 
NID2=NIN/2 
DO 20 I=1,NID2 

C(I)=FLOAT(IC(D) 
20 CONTINUE 

RETURN 
END 

C 

SUBROUTINE CENTER (A,IMAX) 

DOUBLE PRECISION A (1024) ,X f T(21) ,Y 

X=0.D0 

DO 10 1=300,800 

IP (A (I) .GT. X) THEN 

X=A(I) 

J-I-10 

END IF 
10 CONTINUE 
K=J+21 
Y=0.D0 
DO 20 I=J,K 

L=I-5 

M=I+5 

T(I)=0.D0 

DO 15 N«L,M 

T(I)=T(I)+A(N) 
15 CONTINUE 

IF (T{I) .GT. Y) THEN 

Y=T(I) 

IMAX=I 

ENDIF 
20 CONTINUE 



RETURN 
END 

C 

SUBROUTINE PIT (N,X,Y,C) 
C THIS SUBROUTINE IS A SIMPLIFICATION OF 
C A LARGER SUBROUTINE WHICH FITS A CURVE USING 
C UP TO A NINTH ORDER POLYNOMIAL. FOR A 
C REFERENCE SEE JPS DISSERTATION. 

DOUBLE PRECISION X(300) ,Y(300) ,C(2) ,A(3,3) ,XN(300) 
DO 10 1=1, N 
10 XN(I)=1.D0 
DO 30 1=1,2 
A(I,1)=0.D0 
A(I,3)=0.D0 
DO 20 J=1,N 

A(I,1)=A(I,1)+XN(J) 
A(I,3)=A(I,3)+Y(J)*XN(J) 
20 XN(J)=XN(J)*X(J) 
30 CONTINUE 
A(2,2)=0.D0 
DO 40 J=1,N 

A(2,2)*A(2,2)+XN(J) 
40 XN(J)=XN(J)*X(J) 
A(1,2)=A(2,1) 
DO 50 1=1,2 
SUM=0.D0 

IF (2 .GT. I) GO TO 45 
SUM=SUM+A ( 1 , 1 ) *A (1 , 2 ) 
A(I,2)=A(I,2)-SUH 
GO TO 50 

45 IF (ABS(A(I r I)) .LT. l.E-10) GO TO 99 

A(I,2)=(A(I,2)-SUM)/A(I,I) 
50 CONTINUE 

DO 90 J=l,2 
90 C(J)=A(J,3) 

C(1)«C(1)/A(1,1) 

SUM*A(2,1)*C(1) 

C(2)«=(C(2)-SUM)/A(2,2) 

SUM=A(1,2)*C(2) 

C(1)=C(1)-SDM 

RETURN 

99 PRINT* , 1 TROUBLE IN SUBROUTINE FIT' 
RETURN 
END 

C 

SUBROUTINE ZEROFILL <C,IM,IE,N) 
DOUBLE PRECISION C (2048) ,D (2048) , T 
M=IE-IM 
MP1=M+1 
DO 10 1=1, N 
D(I)=0.D0 



10 CONTINUE 

DO 20 1=1, MP1 
J-IM+I-1 

D(I)*C(J) 
20 CONTINUE 

DO 30 1=1, N 
C(I)=D(I) 
30 CONTINUE 
RETURN 
END 

SUBROUTINE TAPER (C,M,N) 
DOUBLE PRECISION C(2048),T 
T»l. D0/1000. DO 
DO 90 X-1,M 

IF (C(I) .NE. 0.D0) 60 TO 80 

J=I-1 

K«I+1 

C(I)«(C(J)+C(K))/2.D0 
IF (C(I) .LT. T) GO TO 95 
80 CONTINUE 
90 CONTINUE 
GO TO 97 

95 IK»I 

DO 96 J=IK,M 
K=J-1 

C(J)=0.95D0*C(K) 

96 CONTINUE 

97 CONTINUE 
ND2=N/2.D0 

DO 100 1=2, ND2 

J=ND2+I 

K=ND2+2-I 

C(J)*C(K) 
100 CONTINUE 
RETURN 
END 

SUBROUTINE SMOOTH (C.M,N) 
DOUBLE PRECISION C(204 8) 
DO 20 1=2, M 
J=I-1 

IP (C(J) .GT. C(I)) GO TO 10 
K=I+1 

C(I)=(C(J)+C(K))/2.D0 
5 IF (C(J) .GT. C(D) GO TO 10 
C(I)»0.999*C(I) 
GO TO 5 
10 CONTINUE 
20 CONTINUE 
ND2=N/2 
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DO 30 1=2, ND 2 

J-ND2+I 

K=ND2+2-I 

C(J)«C(K) 
30 CONTINUE 
RETURN 
END 

C 

SDBROOTINE REORDER (M,C) 
C IN ORDER TO USE TEE FOLLOWING PFT ONE MUST 
C FIRST* REORDER THE DATA. SEE JPS DISSERTATION 
C FOR DETAILS. 

COMPLEX *1 6 C(4096),T 

N«2**M 

MV2«N/2 

NMl=N-l 

J=l 

DO 7 1=1, NM1 

IF (I .6E. J) GO TO 5 
T=C(J) 
C(J)»C(I) 
C(I)=T 
5 K=MV2 

€ IF (K ,GE. J) GO TO 7 
J=J-K 
K=K/2 
GO TO 6 
7 J=J+K 
RETURN 
END 



SUBROUTINE FFT (K, A) 
C FOR DETAILS SEE JPS DISSERTATION 
DOUBLE PRECISION PI 
COMPLEX*16 A(4096) ,U,W,T 
PI=3.141592653589793D0 
N=2**M 
DO 30 L=1,M 

LE=2**L 

LEI=LE/2 

U=(l.D0,0.DO) 

W=DCMPLX ( DCOS (PI/LEI ) ,-DSIN (PI/LEI) ) 
DO 20 J=1,LEI 

DO 10 I=J,N,LE 
IP=I+LEI 
T=A(IP)*U 
A(IP)=A(I)-T 
A(I)=A(I)+T 
10 CONTINUE 
U=U*W 



20 CONTINUE 
30 CONTINUE 

RETURN 

END 

C 

SUBROUTINE CONV (M,F,R,H) 

DOUBLE PRECISION X 

COMPLEX *1 6 H(4096) ,F(4096) ,R(4096) 

N=2**M 

Wl«PLOAT(N) 

DO 10 1=1, N 

B(I)«P{I)*R(I) 
' 10 CONTINUE 

DO 20 1=1, N 

H(I)«H(I)/W1 
20 CONTINUE 
RETURN 
END 

C 

SUBROUTINE DECON (M,F,R,H) 
DOUBLE PRECISION X 

COMPLEX*16 F(4096) ,R(4 096) ,H{4096) 
N=2**M 

C IN POURIER SPACE DECONVOLUTION IS SIMPLY A 
C COMPLEX DIVISION. 
DO 10 1=1, N 

P(I)=E(I)/R(I) 
10 CONTINUE 
X= FLOAT (N) 
DO 30 1=1, N 

F(I)=F(I)/X 
30 CONTINUE 
RETURN 
END 

C 

SUBROUTINE ITER (M,R,F,B) 
DIMENSION A(2048) ,B(2048) 

COMPLEX*16 R(4096) ,F(4096) ,H(4096) ,H1(4096) ,T 
CHARACTER* 7 NAME 



PRINT* / 
PRINT* , 
PRINT*, 
PRINT*, 
PRINT* , 
PRINT* , 
PRINT* , 
READ (* 
PRINT* , 
PRINT*, 
PRINT*, 
READ (* 



ENTER TBE VALUE FOR THE PARAMETER WHICH' 

CONTROLS THE ALLOWABLE ERROR BETWEEN THE 

OBSERVED SPECTRUM AND F*R. • 

EPS=THE SORT OF THE SUM OF THE SQRS OF' 

THE DIFFERENCES BETWEEN THE H(I)S AND* 

THE F(I)*R(I)S. ' 

(A REAL NUM3ER PLEASE.) 1 

•(F13.7)') TPS 

ENTER THE MAXIMUM NUMBER OF ITERATIONS' 
YOU WISH TO PERMIT. • 
(AN INTEGER PLEASE) • 
' (13) •) I END 
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N«2**M 

ND2-N/2 

ND4-N/4 

DO 30 J>1,IEND 

CALL REORDER (M f F) 

CALL PPT (M,F) 

CALL CONV (M,P,R,B1) 

CALL REORDER (M,P) 

CALL PPT (M,F) 

CALL REORDER (M,B1) 

CALL PPT (M,B1) 

RSQ=0.0 

DO 2 IA=1,N 

A ( IA) =REAL (B ( IA) ) 

B(IA)«REAL(E1(IA)) 
2 CONTINDE 
V=0.0 
W=0.0 

DO 4 IAA=1,N 

V=V+A(IAA) 

W-W+B(IAA) 
4 CONTINUE 

DO 6 IAB=1,N 

A(IAB)=A(IAB)/V 

B(IAB)=B(IAB)/W 
6 CONTINDE 

DO 10 IAC=l f ND4 

RSQ=RSQ+ (A ( I AC) -B ( IAC) ) * *2 

10 CONTINUE 
RSQ=RSQ**0.5 
WRITE (*,11) RSQ 

11 FORMAT (SX^RSQ- ' ,E13.6) 

IP (RSQ .LT. EPS) GO TO 33 

P(1)=P(1)*(A(1)+2*A(2)+2*A(3))/(B(1)+2*B(2)+2*B(3)) 

P(2)=P(2)*(A(1)+2*A(2)+A(3)+A(4))/(B(1)+2*B(2)+B(3)+B(4)) 
DO 20 IB=3,ND4 

J=IB-2 

K=IB+2 

TB-0.0 

TH1-0.0 

DO 15 II=J,K 

TB=TB+A(II) 

TH1=TH1+B(II) 

15 CONTINUE 

IF (TBI .EQ. 0.0) GO TO 20 

F(IB)*F(IB) *TB/TB1 

20 CONTINUE 

DO 25 IC»2,ND2 

J=ND2+IC 

K=ND2+2-IC 

F(J)=F(K) 



c 



25 CONTINUE 
30 CONTINUE 
33 CONTINUE 
1=1-1 

WRITE (*,35) I 

<5X ' ,THE N °- 0P ite «a*ions is M; 

END 

SUBROUTINE P ARAMS (N,A) 

SSSSrKif'SSS A ' 2046) • X(12) •«<"> - V,K 
K=1024 

DO 20 1=1,12 
X(I)»0.D0 
DO 10 J=l,5 
K=K+1 

X(I)-X(I)+A(K) 
10 CONTINUE 
20 CONTINUE 
V=0.D0 

DO 30 1*1025,1084 
V=V+A(I) 
30 CONTINUE 

DO 40 1=1,12 
X(I)=X(I)/V 
40 CONTINUE 

DO 60 1=2,10,2 
M(I)=0.D0 
DO 50 J=1026,1084 
V=FLOAT(J-1025) 
W-FLOAT(I) 
V=V* *W 

en ^ M(I)*M(I)+A(J)*V 
50 CONTINUE 

M(I)*M(I)/59.D0 
60 CONTINUE 

M(2)=DSQRT(M(2) ) 
V=M(2) 

DO 65 1=4,10,2 

M(D=M(I)/(V**I) 
65 CONTINUE 1 

»S™!'! ENTER raE FILENAME FOR THIS SET OF- 
PRINT*, 'PARAMETERS. • OF 
READ (*,'(A7)') NAME 

OPEN (90,FILE=NAME,PAD= , YES I ) 
WRITE (90,70) 

70 FORMAT (5X, 'EVEN CENTRAL MOMENTS') 

WRITE (90,80) <I,M(I), 1=2 8 2 \ 
80 FORMAT (10X, '„ t \ 

WRITE (90,100) (M(10)) 



100 FORMAT (9X,»M(10) « a ,D17.10) 

WRITE (90,110) 
110 FORMAT (5X,'5 CHANNEL GROUPS (CENTER FIRST) 

WRITE (90,120) (I,X(I),I»1,9) 
120 FORMAT (10X, 'X( \I1, ') =',D17.10) 

WRITE (90,130) (I,X(I),I=10,12) 
130 FORMAT (9X, 'X( ' ,12, ■ ) = ',D17.10) 

CLOSE (90) 

RETURN 

END 

SUBROUTINE BACK (A,IMAX,N) 
DOUBLE PRECISION A (1024) ,C (2048) , V 
COMPLEX *1 6 B(4096) 
CHARACTER* 7 NAME 
ND2=N/2 
ND4=N/4 
DO 10 1-1, N 
C(I)=0.D0 
10 CONTINUE 
N1=ND2-IMAX 

IF (Nl .GT. ND4) N1=ND4 
DO 20 1=1, Nl 

J=IMAX+I-1 

C(I)=A(J) 
20 CONTINUE 
N2=IMAX-1 
N3=N-N2 
DO 30 J=1,N2 

I=N3+J 

C(I)»A(J) 
30 CONTINUE 
DO 40 1=1, N 

B(I)=DCMPLX(C(I) ,0.D0) 

40 CONTINUE 
M=ll 

CALL REORDER (M,B) 
CALL FFT (M,B) 

PRINT*, 'WOULD YOU LIKE TO SEE THE FIRST FEW' 

PRINT*, 'FOURIER COEFFICIENTS? ' 

READ (*,' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 43 

PRINT* , ' HOW MANY? • 

READ (*,*) IFC 

WRITE (*,41) (I,B(I) ,1=1, IFC) 

41 FORMAT (I4,2D22.10) 
43 CONTINUE 

PRINT*, 'WOULD YOU LIKE TO SEE THE FIRST FEW' 
PRINT* , ' VALUES OF THE POWER SPECTRUM? 1 
READ (*,' (Al) ') NAME 
IF (NAME .EQ. 'N') GO TO 46 



PRINT*, •HOW MANY? • 
READ (*,*) IPS 
DO 44 1=1 , IPS 

A(I)«B(I)*CONJG(B(I)) 

44 CONTINUE 

WRITE (*,45) (I, A (I) ,1=1, IPS) 

45 FORMAT (I4,D22.10) 

46 CONTINUE 

PRINT* i ' DO YOU WISH TO DELETE ALL OP THE SINE* 

PRINT* , • COEFFICIENTS? • 

READ (*, ' (Al) •) NAME 

IF (NAME . EQ. 'N') GO TO 51 

DO 50 1=1, ND 2 

B(I)«(B(I)+CONJG(B(I) ))/2.D0 

50 CONTINUE 

51 PRINT*, 'DO YOU WISH TO DELETE SOME OF THE COSINE' 
PRINT* , ' COEFFICIENTS? • 

READ (*, ' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 53 

PRINT* , ' HOW MANY? • 

READ (*,*) IDC 

DO 52 1=1, IDC 

B(I)=(B(I)-CONJG(B(I)))/2.D0 

52 CONTINUE 

DO 152 1=2, IDC 
J=2050-I 

B(J)=(B(J)-CONJG(B(J) ) )/2.D0 
152 CONTINUE 

53 PRINT*, 'INSTEAD OF DELETING SOME OF THE COSINE' 
PRINT*, ' COEFS. YOU MAY PREFER TO JUST TRIM THEM. ' 
PRINT*,' THIS IS DONE BY SETTING THE FIRST N TERMS 
PRINT*, 'EQUAL TO THE N+l TERM. YOU MAY INVOKE' 
PRINT*, 'THIS OPTION EVEN IF YOU HAVE PREVIOUSLY' 
PRINT*, 'DELETED THE FIRST N COSINE COEFFICIENTS.' 
PRINT*, 'DO YOU WISH TO DO SO? • 

READ (*, ' (Al) «) NAME 
IF (NAME . EQ. 'N') GO TO 200 
PRINT*,' HOW MANY? ' 
READ (*,*) IT 
DO 190 1=1, IT 
J=IT+1 

B(I)=DCMPLX(REAL(B(J) ) ,0.D0) 
190 CONTINUE 
200 CONTINUE 

PRINT*,' DO YOU WISH TO DELETE SOME OF THE SINE' 

PRINT*, 'COEFFICIENTS? ' 

READ (*,' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 55 

PRINT* , ' HOW MANY? » 

READ (*,*) IDS 

DO 5 4 1=1, IDS 



B(I)«=(B(I)+CONJG(B(I)))/2.D0 

54 CONTINDE 

DO 154 1=2, IDS 
J-2050-I 

B(J) = (B(J) +CONJG(B(J)) )/2.D0 
154 CONTINUE 

55 CONTINUE 

PRINT* , * DO YOO WISH TO SAVE THE FOURIER COEFS. 

PRINT*, 'B(X) PRIOR TO THE TAPER IN THE WINGS? 

READ (*,' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 57 

PRINT*, 'ENTER THE FILENAME. • 

PRINT*,' 1 

READ (*,'(A7)') NAME 

OPEN (90, FILE-NAME, PADs* YES') 

WRITE (90,56) (I, B (I) , 1=1,2048) 

56 FORMAT (2 (I4,2D22.10) ) 
CLOSE (90) 

57 CALL REORDER (M,B) 
CALL FFT (M,B) 

DO 60 1=1, ND2 

C(I)=REAL(B(I)) 
60 CONTINUE 
V-0.D0 

DO 70 1=359,484 
V=V+C(I) 
70 CONTINUE 
V«V/126. 

PRINT* , 1 DO YOU WANT A D.C. CORRECTION USING A' 
PRINT*,' LEAST SQUARE FITTING IN THE TAILS? ' 
READ (*,' (Al) ') NAME 
IF (NAME .EQ. 'N') V=0.D0 
DO 80 1=1, ND 2 
C(I)=C(I)-V 

IF (C(I) .LT. 0.D0) C(I)=0.D0 
80 CONTINUE 

DO 90 I=1,ND4 
J=(3*ND4)+I 
A(I)=C(J) 
90 CONTINUE 
N4=ND4+1 
DO 100 I=N4,ND2 
J=I-ND4 
A(I)«C(J) 
100 CONTINUE 
IMAX=513 
RETURN 
END 



SUBROUTINE NOBACK (N,B,IMAX) 
DOUBLE PRECISION B (1024) ,A(2048) ,V 



C0MFLEZ*16 C(4096) 
CHARACTER*7 NAME 
M»ll 
ND2=N/2 
ND4=N/4 
DO 10 1=1, N 
A(I)=0.D0 
10 CONTINUE 
N1=ND2-IMAX 

IF (Nl .GT. ND4) Nl=ND4 
DO 20 1=1, Nl 
J=IMAX+I-1 
A(I)«B(J) 
20 CONTINDE 
N2-IMAX-1 
N3=N-N2 
DO 30 J=1,N2 
I=N3+J 
A(I)=B(J) 
30 CONTINDE 
V=0.D0 
DO 33 1=1, N 
V=V+A(I) 
33 CONTINUE 
DO 36 1*1, N 
A(I)=A(I)/V 
36 CONTINUE 

OPEN (90,FILE=*RAWDATA' ,PAD=* YES' ) 
WRITE (90,70) (I,A(I) ,1=1, N) 
CLOSE (90) 
DO 40 1=1, N 

C(I)=DCMPLX(A(I) ,0.D0) 
40 CONTINUE 

CALL REORDER (M,C) 
CALL FFT (M,C) 
DO 50 1=1, N 

A(I)=CXl)*CONJG(C(I) ) 
50 CONTINUE 

PRINT* , ' ENTER THE FILENAME FOR E{X). • 
PRINT* $ ' ' 
READ (*, • (A7) ') NAME 
OPEN ( 90 # FILE-NAME. PAD*' 7ES') 
WRITE (90,60) (I ,C(I) ,1=1,2048) 
60 FORMAT (2 ( 14 ,2D22 .10) ) 
CLOSE (90) 

PRINT* , ' ENTER TEE FLENAME FOR THE P.S. 
PRINT* , ' ' 
READ (*, '(A7) ') NAME 
OPEN (90 ,FILE=NAME,PAD= ' YES' ) 
WRITE (90,70) (1,A(I), 1=1, 1024) 
70 FORMAT (5 ( 14 ,D22 .10) ) 



CLOSE (90) 

RETURN 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccc 



APPENDIX G 



In this Appendix the theoretical relationship between the vacancy 
concentration, C v , and the concentration of dislocation cores, C d , is 
demonstrated, (from Cotterlll et a]_. 1972). The starting equation is 

C v * A exp(-E vf /KT) (G.i) 
where C v 1s the vacancy concentration at thermal equilibrium 

E v f is the vacancy formation energy 
and KT has its usual meaning 

Now the concentration of vacancies after the quench, c', is related to C„ by 

v v 

C v s nC v (6.2) 
where n is the quench efficiency which is assumed to be nearly unity 

After the low temperature anneal the concentration of dislocation 
loops, Cg, is given by 



C i s < (6.3) 



where v is the inverse of the number of vacancies required to form an 
average loop 

The dislocation concentration, C d , in terms of an atomic fraction 1s then 

C d - C A *D/d 0 (G.4) 
where D is the average loop diameter 
and d 0 is the nearest neighbor distance. 
It can be shown that 
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v = 2fl d*/irD 2 (s.5) 
Combining Equations (G.3), (G.4), and (6.5) yields 

C d = 20 d Q C;/D (6.6) 



APPENDIX H 

In this Appendix the derivation of the defect specific lineshape 
parameter "R" is recreated from the work of Mantl and Trifthauser (1978). The 
following quantities require definition: 

S.W - are the values of the "S" and "W" parameters for the sample 
being analyzed. 

S t ,W t - are the values of the "S" and "W" parameters for a sample in 
which all of the positrons are annihilating from the trapped 
state. 

S B ,W B - are the values of the "S" and "W" parameters for a sample in 
which all of the positrons are annihilated from the bulk 
(defect free). 

P t - the probability of a positron annihilation from the trapped 
state. 

P B - the probability of a positron annihilating from the bulk. 
m - the trapping constant. 
C - the trap concentration. 
*B " tne annihilation rate from the bulk state. 
Now, it is known that 

S = S B P B + S t P t (H .i) 

and 

p B + p t ' 1 (H.2) 
(143) 



Also P t s »»C/(A B + W C) 

Substituting Eq. (H.3) Into Eq. (H.l) yields: 

wC 



S = 



Similarly 



W - W K (1 - 



mC 



b Vi ~ A. + jiC 



) + K ( 



PC 



t v A f + yC 



Rearranging Eq. (H.4) and Eq. (H.5) yields: 

S - % - (S t - S B )pC/(A f + pC) 
W - W B = (W t - W B )uC/(A f + V C) 

Taking the ratio of Eq. (H.6) to Eq. (H.7) yields: 



R = 



IS - S. 



W - W 



S t" S B 



W t ~ W B 
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(H.3) 

(H.4) 



(H.5) 

(H.6) 
(H.7) 

(H.8) 



Clearly the ratio on the far right has a constant value for a given defect in 
a yiven material. 
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Porous silicon is a highly complex system containing 
disorder that exists over a Tange of scales. Depend- 
ing on the anodization current and time, and the 
dopant level and type, the resulting porous layer can 
have a porosity of up to 80% with a unidirectional 
columnar-like structure for p~ silicon wafers [I]. 
The remaining silicon varies in size with character- 
istic lengths ranging between from 3-4 to 15- 
30 nm. Such small crystalline structures have an 
enlarged bandgap due to the quantum confinement 
effect. The role of confinement in the observed 
increase in radiative inter-band transition efficiency 
is however still unclear. How the broken silicon 
bonds are terminated can directly affect the nature of 
the inter-band transitions. For example, hydrogen 
passivation "sweeps" the dangling bond states out of 
the band gap region; reducing the available non- 
radiative recombination paths and thereby increasing 
the efficiency of the radiative transitions. This, 
however, is not the complete picture since the 
oxygen content of the porous region has been found 
in some cases to amplify and change the photo- 
luminescence properties [2]. The electro-chemical 
processes occurring on the surface during anodiza- 
tion are not yet clearly understood and further 
experimental characterization is needed. 

Positron annihilation lifetime spectroscopy (PALS) 
provides a unique non-destructive probe into the 
disorder inherent in porous silicon on the micro- 
scopic scale, through the identification of impurity 
concentrations and types; and on the mesoscopic 
scale at which the bulk of the porosity exists through 
the study of positronium lifetimes. Itoh et al. [3,4] 
and Dannefaer et a!. [5] have done such an 
investigation using a discrete lifetime analysis. Itoh 
et al. obtained a satis factory three-component fit, 
finding a bulk silicon component (at 230 ps), a 
defect component (at 600 ps), which they interpreted 
as arising from either silicon mono-vacancies and/or 
oxygen divacancies within an Si-02 layer [6], and a 
longer 25 ns lifetime due to ortho-positronium 
annihilating in the porous region. Dannefaer et ai y 
on the other hand, determined that a five-component 
fit was better, with additional lifetimes arising from 
para-positronium (at 1 25 ps) and another ortho- 
positronium component at 2.5 ns. They attributed 
the latter to ortho-postronium annihilating in larger 
vacancy defects either within the remaining silicon 
or the oxide layer. In this letter we report on a PALS 
investigation of porous silicon using a continuous 
lifetime fitting algorithm. 

0261-8028 y IM6 Chapman & Hail 



Our motivation lies in the underlying disadvantage 
in discrete lifetime fitting algorithms where the 
number of components must initially be assumed 
since in general a realistic spectrum does not 
uniquely determine this number. This becomes 
particularly apparent when looking at highly dis- 
ordered systems where, indeed, the notion of a 
discrete spectrum may be invalid. These problems 
can in part be overcome by using a continuous 
lifetime fitting algorithm; for example, performing 
an inverse Laplace transform on the lifetime 
spectrum, thereby obtaining a probability density 
function of characteristic lifetimes. The Fortran 
program CONTIN [7-9] is such an algorithm. 
Admittedly, use of the inverse Laplace transform 
does not take into account the instrument resolution 
function of the apparatus. For the determination of 
sub-nanosecond lifetimes this entails a distinct loss 
of accuracy. For well-separated lifetimes that are 
significantly larger than the typical full width at half 
maximum (FWHM) of the equipment, the effect is 
reduced and overshadowed by the loss of accuracy 
due to the worsening statistics typical for these 
components. To date, our work appears to be the 
only systematic study of porous silicon using a 
continuous lifetime fitting algorithm. This is the first 
step towards achieving the goal of determining 
intrinsic lifetime distribution widths for positrons 
annihilating within the porous silicon layer, and thus 
its electronic structure. 

For the measurement of the annihilation lifetime 
spectra of porous silicon we employed a standard 
PALS rig at 200 ps per channel. Such a large timing 
window is needed to ensure that all lifetime 
components decay completely into the background. 
The positron source was 22 NaC! evaporated onto a 
6 fim thick titanium film. Using our PALS rig at 
50 ps per channel and the discrete lifetime fitting 
program PFPOSFIT [10], we obtained for bulk 
silicon, 221 ± 1 ps at 97 ± 1% concentration, with 
a source component lifetime of approximately 
429±25ps at 3 ± 1% concentration (the iterated 
instrument resolution had an FWHM of 240 ps with 
long and short time tails of 50 and 20 ps r 
respectively. 

The porous silicon samples for this investigation 
were made from 200 urn p~ doped wafers, anodized 
for 30 min at 20 mAcm -2 . The electrolyte used was 
a 35% ethanoic solution of HF. Under illumination 
by an ultraviolet (UV) lamp this sample glowed a 
deep red. The initial thickness of the porous layer 
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was approximately 50-60 wm. A high resolution 
scanning electron microscopy (SEM) analysis per- 
formed approximately 6 months after the PALS 
experiment revealed a layer that was approximately 
10 fim thick, indicating that structural damage had 
occurred (Fig. 1). The sample nevertheless still 
glowed red at a reduced intensity. Fig. 1 reveals the 
large-scale columnar structure typical of p" doped 
porous silicon samples; the remaining (porous) 
silicon constituting the much thinner connecting 
membranes between the micrometre-sized pores. 

All PALS runs were taken at room temperature 
(24 °C) in an environment where the temperature 
could be held constant to within a degree. For each 
sample, a varying number of spectra were taken (see 
Table I and Fig. 2), each of which had a maximum 
peak height of 10 6 counts. Statistical analysis across 
the spectra for a particular sample indicated that to 
within the noise level, there was no detectable 
difference between individual runs, justifying the use 
of the mean spectrum in the analysis using CONTIN. 
The associated information content analyses [11, 12] 
are also given in Table I. Here, the Shannon number 
indicates that in principle up to eight distinct lifetime 
peaks are resolvable, and the resolution factor 
suggests that the minimum ratio between these 
resolvable lifetimes ranges between 2 and 3. 

Since CONTIN assumes a spectrum constructed 
from exponentials, we must choose an appropriate 
zero-time channel from which the analysis can 
begin. This must be just beyond the peak of the 
spectrum, yet not too close to include the "rounding 
off' effect of the finite instrument resolution. The 
resulting inverse Laplace transform found by CON- 
TIN is not very sensitive to the choice of the zero- 
time channel as far as determining distinct compo- 
nents. It is, however, sensitive to determining their 



intensity, since a shift 5t in the zero-time channel 
will affect the intensity of a component r by the 
approximate factor exp(-<3//r). 

Such an analysis for the 200 ps per channel bulk 
silicon wafer spectrum is displayed in Fig. 2 for two 
different zero-time channels; approximately 400 ps 
(solution I) and 800 ps (solution 2) beyond the 
instrumental zero-time channel. Solution 1 contains 
a dominant broad component centred at approxi- 
mately 220 ps, which can be attributed to* the bulk 
silicon component. Two long time components can 
be seen, occurring at approximately L5 and 45 ns, 
the shorter of which can be associated with the 
common ortho-positronium "pick off'* signal present 
in semiconductors [13]. The extra long 45 ns 
component does not seem to have been encountered 
in the literature. We do not know at this stage the 
origin of this lifetime peak; however, as we shall see, 
it also appears in the porous silicon sample at an 
increased intensity. In solution 2, the bulk silicon 
component becomes broader and less intense, the 
ratio of its peak height with that of solution 1 being 
approximately exp (it/220) where 67»400ps. The 
ortho-positronium peaks at 1.5 and 50 ns acquire 
more dominance, with their peak heights largely 
unchanged from solution 1. In addition, we observed 
a shorter, approximately 125 ps, component, which is 
most likely to be the accompanying para-positronium 
signal. That this was not resolved in solution J 
suggests that the "rounding off" effect was not 
completely removed at 400 ps beyond the approx- 
imate zero-time channel. 

Fig. 3 displays the CONTIN solution for the 
porous silicon spectra. Again, solution 1 and solution 
2 are for time shirts of approximately 400 and 
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Figure I High resolution SEM picture of the porous silicon sample 
used in the PALS investigajion. . 



Figure 2 PALS lifetime spectra of (3) bulk silicon and (b) porous 
silicon taken al 200 ps per channel. 



TABLE I The number of PALS spectra raken and ihc resulting resolution factor 2nd Shannon 
number 

Sam P ,e Number of Resolution Shannon 

spectra factor number 

Bulk silicon at 200 ps per channel 49 2 9 

Porous silicon at 200 ps per channel Z4 3 8 
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Figure 3 Characteristic lifetime probability density determined using 
CONTIN for the bulk silicon sample taken at 200 ps per channel, (a) 
Solution I: zero-lime approximately 400 ps beyond the instrumental 
zero-time channel, (b) Solution 2: zero-time approximately 800 ps 
beyond the instrumental zero-time channel. 



800 ps beyond the instrumental zero-time channel, 
respectively. Inspection of solution 1 indicates five 
distinct components. Component 1 is centred at 
approximately 200 ps and can be attributed to the 
positrons annihilating either within the silicon 
substrate or in the remaining silicon in the porous 
material. Component 2 is centred at approximately 
300 ps and has also been encountered by both hoh et 
al. and Dannefaer et ah Their measured lifetimes for 
component 2» however, ranged between 400 and 
600 ps. CONTIN analyses of simulated spectra 
obeying Poisson statistics indicate that this discre- 
pancy (and the reduced bulk lifetime) may be due to 
the 200 ps per channel resolution and the neglect of 
the instrument resolution function. Components 3, 4 
and 5, occurring respectively at approximately 2, 10 
and 50 ns, are the long lifetime components due to 
ortho-postronium. In solution 2, the sub-nanosecond 
components are reduced in intensity, and crystalline 
silicon and defect components (at 300 ps) are no 
longer separately resolvable, merging into a single 
broad peak centred at approximately 250 ps. The 
long time components now become broader and of a 
reduced height. This indicates that the widths we 
have observed are not intrinsic to the material but 
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to noise ratio. Marshall [14] has concluded that to 
determine accurately the intrinsic widths of peaks, 
this ratio must be at least 400. Inspection of Fig. 4 
indicates a signal to noise ratio of less than 40 for 
the positronium decay modes. 

The present method of analysis is unsuitable to 
obtain accurate absolute intensity information since 
the zero-time channel cannot be estimated beyond 
the 200 ps per channel limit Furthermore, modifica- 
tions to CONTIN suggested by Gregory [12, 15-17] 
that take into account the finite resolution (and 
therefore the zero-time channel), by obtaining the 
probability density profile with respect to a reference 
spectrum (in our case, bulk silicon), have proved 
unsuccessful. This is primarily due to the poor 




0 200 400 600 

Channel number (at 200 ps per channel) 

Figure 4 Characteristic lifetime probability density determined using 
CONTIN for the porous silicon sample taken at 200 ps per channel, (a) 
Solution 1 : zero-time approximately 400 ps beyond ihe instrumental 
zero-time channel, (b) Solution 2: zero-time approximately 800 ps 
beyond the instrumental zero-time channel. 

definition of the bulk silicon reference spectrum at 
200 ps per channel. 

Using the PALS data collected from sample A, we 
have confirmed and extended the findings of both 
Itoh et al. and Dannefaer et al. In particular, we 
resolved three rather than two ortho-positronium 
components, suggesting that there may be an 
additional intermediatory scale of porosity in which 
ortho-postronium annihilates. This has a lifetime 
indicative of hydrogenated amorphous silicon [18], 
suggesting ortho-positronium annihilation again 
within the amorphous layer. The integrated areas 
(corrected for zero-time shift) of the positronium 
peaks in solutions I and 2 show that components 3 
and 5 occur in approximately the same proportion; 
Dannefaer et al. % on the other hand, found that the 
—2.5 ns component occurred an order of magnitude 
less than the longest time ortho-positronium compo- 
nent The present analysis also establishes the 
existence of a very weak ortho-positronium compo- 
nent in the pre-anodized wafers at a time scale 
approximately equal to the longest time ortho- 
positronium component seen in porous silicon, using 
theory [19] that has traditionally been employed in 
the calculation of free volumes in polymer materials. 
In porous silicon Itoh et al. associated a void length 
scale of approximately I nm to this ortho-positro- 
nium component, suggesting that irregularities of 
this magnitude exist on the pre-anodized wafer and 
that these may. in part, be the catalyst for the initial 
pore formation process. 

The results presented here further demonstrate the 
rich variety of disorder inherent in porous silicon 
and indicate the potential of a continuous lifetime 
analysis of PALS spectra to provide a high resolution 
probe into such structures. If another order of 
magnitude increase in signal to noise ratio can- be 
achieved then there exists the possibility of resolving 
intrinsic peak widths for porous silicon samples, 
thereby extending the PALS method to resolving 
pore size distributions using the ortho-positronium 
components, and conceivably giving insight into the 
electronic structure of the confined electrons within 
the remaining silicon crystallites of the porous layer. 
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Abstract 

The study of deconvolution by Van Gtteit and Gold iterative algorithms and their use in the processing of experimental 
spectra of Poppler broadening of the annihilation line in positron annihilation measurement is described. By com pa ring 
results from both algorithms it was observed that- the Gold algorithm was able to eliminate linear instability of the 
measuring equipment if one uses the 1274 keV a Na peak, that was measured simultaneously with the amtihilation peak, 
for deconvolution of annihilation peak 311 keV. This permitted the measurement of small changes of the annihilation peak 
(eg. ^parameter) with high confidence. The dependence of /-ray-like peak parameters on the number of iterations and 
the ability of these algorithms to distinguish a y-ray doublet with different intensities and positions were also studied. 



1. Introduction 

Positron annihilation (PA) as one of the non-destructive nuclear methods is used for the study of structure 
of different materials and it is especially efficient in Determining free volumes or defects in bulk material. The 
main characteristics, that are determined by PA, are lifetimes of positrons in matter, Doppler broadening of 
annihilation line (DBAL), and angular correlation of annihilation /-rays. The last two characteristics give us 
information about the momentum distribution of electrons in matter. The more common angular correlation 
method is approximately one order more accurate, but it is more time consuming and requires a more intensive 
positron source. On the other hand DBAL is simple and not time consuming. It is mainly used for o>termining 
the relative changes in the annihilation peak (e*g* S-parameter, defined as the ratio between area of central part 
of the peak to the full area of the peak) depending on temperature, pressure etc After deconvolution of the 
annihilation peak using a resolution function of measuring system, one can obtain the momentum distribution 
of electrons in the sample; If we are only interested in relative changes of the ^-parameter, deconvolution is 
not necessary. But as we have seen in our work, deconvolution can partly eliminate the spurious effects due to 
instability of the measuring system. Another possibility to eliminate an instability of the system is the use of 
ss expensive spectnnss stabilizer unit 

Generally, deconvolution algorithms have recently found many applications in various domains of experi- 
mental science. They have been applied to such different problems as improving resolution in spectroscopy and 
measuring the thickness of multilayer structures. There is also the long-standing problem of determination of 
position and intensities in y-ray multiples [1]. 

2. Theory 

The relationship between a measured value x(f) and the raw result of measurement y(t) can be described 
by a convolution-type integral equation 

O168-90QZ/97/St7.0O Copyright © 1997 Bsevier Science B.V. Ail rig&u reserved 
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y(0 



oo 



x(T)h(t-r) dr, 



(I) 



where h(t) is an impulse response. Knowledge of the instrumental function h(t) is usually required. For a 
discrete system, Eq, ( 1 ) can be written as 



y(0 = *)*<*) » aB 0.I..;. t 2W-2 1 



(2) 



t=0 



where N is the number of samples of vectors A, x. 

The impulse response has a finite length. Therefore we will assume that h(i) = 0 for i < 0 and i ^ N. Then 
Eq. (2) can be written in a matrix form: 



r .v(0) 
.v(D 
>(2) 



or 



.y(2V-2) 



y = Hx. 



h(0) 0 0 
MD MO) 0 
h{2) A(I) h(0) 



h(N-\) A(/V-2) A(rY-3) ... 
0 A(/V-l) h(N-2) ... 
0 0 h(N-\)... 



x(0) 
x(l) 
x(2) 

x(*-l) 



(3) 



(4) 

It means that the columns of H are represented by vectora h mutually shifted by one position. Multiplying 
both sides of Eq. (4) by H T gives 



or 



H T y = H T Hx. 



jn^Hix, 



(5) 



(6) 



where Hi is Tbeplitz matrix [231- 

Solution of die linear equation system (6) (vector x), under the condition that the output vector y and the 
matrix of the impulse response H is known, is a problem of deconvolution. The output vector of system y is 
affected by noise that accompanies each measurement The existence of this noise strongly affects the process 
of deconvolution, and can lead to difficulties in solving the linear equation system (6). 

The above-formulated problem of input reconstruction is as a rule ill-conditioned, i.e. the estimates x(t) of 
x(f) satisfying Eq. (1) are extremely sensitive to errors in the measured data jr(r). It is expressed by the fact, 
that the matrix Hj is almost singular. The direct inversion of Hi for solving x cannot lead to a stable solution.. 

Therefore, in order to solve this problem, the method of regularization must be included. This means that 
original problem is replaced by an approximate one which solutions are significantly less sensitive to errors in 
the data /(f). 

The Van Cittert iterative method of deconvolution is widely applied in different areas, for example in 
spectroscopy or in image processing [2,4]. The Van Cittert algorithm of deconvolution is described in detail in 
Ref. [4], so we wilt describe it only very briefly. 

Its basic form for a general linear discrete system is 



jc* + ' =x k +rfy- Ax*), 



(7) 
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where A is system matrix, k represents the number of iterations and /* is the relaxation factor. The convergence 
condition of Eq. (7) is that the diagonal elements of the matrix A satisfy 

It is obvious that such a diagonal element dominance is rare in physical problems. However, the deconvolution 
algorithm in Eq. (7) can be modified in a such way that it will satisfy the conditions of convergence. Hence, 
Eq. (7) becomes 

x M ^M+iE-tiAlx'^My + Dx'. (9) 

where E is a unit matrix and 

D = E-/iA- (10) 
Under the condition that x° = jx, the successive substitutions give 

x*=^ + /*I>> + ...+/iD*" , J + ^ = /*(E + 0+... + D*)jr. (U) 

Assuming that Ao, A i% . . . , are eigenvalues of A. then (I - >Uo),(l -^Ai) (1 -/iA^-i) are 

eigenvalues of D. If 

lim (t-M>* = 0, i-0.l...-.Af-l. (12) 

then 

lim D* = [OJ . 

and 

lim x*«A" I j = x. (13) 

From equation (12) this implies that the necessary and sufficient conditions of convergence are 

|l-/Uf|<l. *' = <U (14) 

If we define A, and iu conjugate A,* as 

then the convergence condition (14) becomes 

A*[/i(a? + *?>-2*] <0, i=0J JV-1. (15) 

Inequality ( 15) gives two bounds for m 
M=0. ^-2-j^— t = 0,l. (16) 

conditions determine the bounds of die /t coefficient Unfortunately these conditions are not fulfilled for 
all cases. However, if the system matrix A is positive definite the convergent solution always exists. So we 
settle the algorithm in a such way that the eigenvalues A$ will be positive, real numbers. 

Let us return to Eq. (5). Matrix H T H is symmetric* so its eigenvalues are real. The eigenvalues of matrix 
(H T H)(H T H) are squares of eigenvalues of mainx H T H and therefore must be positive. Eq. (5) becomes 

(H T HH T ) v= (H T HH T H)x . (17) 
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and ihc iterative algorithm of dcconvolution becomes 
x 4+l = x* + Mf(H T HH T )^- (H T HH T H)x*J , 

or 

x i+, =i* + /i[/-HV]. 



509 



(18) 



(19) 



By this we ensure the existence of a common interval of solution for inequality ( 14) and convergence of the 
deconvolution algorithm. Eigenvalues A/ are real, positive numbers, so from Eq. (16) for /i we can write 



0</t<r— , 

where Amo is the greatest eigenvalue of H # 

Arm =max (Ao,A, A*_i). 

Now we determine the maximum eigenvalue Amu. For eigenvalues of system ( 17) we can write 

Ux = AiX. i = 0. 1 N-l. 

If xj is the biggest absolute element in x, then from Eq. (22) we get 

£H^X m = A,Xy 



(20) 



(21) 



(22) 



(23) 



or 



Then 



(24) 



(25) 



In practical cases we do not know the biggest element in x. We determine the value of Amu as the maximum 
value from A it determined by Eq. (25) eg., from the sum of absolute values of rows in matrix H\ 

This is the base of the Van Gttert algorithm of deconvolution. Now we introduce, in analogy with Eq. (24) , 
a local variable relaxation factor 



(26) 



s=0 



and we use it in Eq. ( 19). For the Ah element of vector x* +I we get 



' V-l 
m=0 



(27) 
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or 



^vZT 1 — 4- (28) 



m=0 



Eq. (28) is the Gold algorithm of deconvolution [5 J. It is an extension of Van Cittert 's iterative method 
The advantage of Gold's method is that the solution x is positive. 



3. Theoretical simulation 



To test an applicability of Van Cittert's and Gold's methods of deconvolution for y-rays spectra we performed 
a theoretical simulation. We convoluted the experimental resolution function (RF) obtained from a positron 
annihilation experiments with different kinds of Gaussians describing y-ray peaks. The simulated spectra were 
used as an input to both Gold's and Van Gttert's iterative deconvolution procedures. Simulated spectra included 
a noise component due to an experimental resolution function. 

3.1. The dependence of output peak parameters from the deconvolution on the number of iterations 

Because the deconvolutions are using the iterative method, it was interesting to study the dependence of the 
output from the deconvolution on the number of iterations {N). N was varied in the range 200-50000. 2000 
iteration steps lasted 40 min on a PC486, 33 MHz for 200 channels of the spectrum (a version with elimination 
of the multiplication of zero matrix terms took approx. 15 min). 

We analyzed the peak which was a convolution of a Gaussian with FWHM = 13 channels and a resolution 
function with a FWHM = 10.4. The resulting peak had a FWHM = 16.7, which was a value in the range of 
our experimental values for the FWHM of annihilation peaks measured in DBAL. The deconvoluted peak was 
fitted and the result should be one Gaussian with a FWHM = 13. 

Fig. 1 shows the dependences of the x* value and FWHM on the number of iterations for the algorithm 
of Vfcn Cittert and Gold. As we can see from this figure, the resulting deconvolution fit for the Van Cittert 
algorithm was much better for two Gaussians than for one until 14000 iterations, then the x 2 became practically 
the same. 

For the Gold algorithm the ^ of the deconvolution fit was very good for one Gaussian and reached a 
minimum value for 2000 iterations. The increase of x 2 with the number of iterations (for both Van Cittert and 
Gold methods) is a consequence of the round-off errors which arise during the computations with real numbers 
represented in standard single precision (REAL in Fortran or FLOAT in C) and a large number of these 




fig. I. Dependence of x 3 ^ FWHM from fit of decoavoluted peak om number of ucntioa stepi. 
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operations. It was verified by a calculation (for the Gold method) using the double precision floating-point 
representation of real numbers. The results of the double-precision calculation are shown in Fig 1. The fast 
convergence and the stability of the solution can be observed. It clearly shows the possible errors arising from 
using a finite computer word length of and a targe number of iteration steps. 

For the Van Cittert algorithm, it can be seen that by increasing the number of iterations the first Gaussian 
became more and more dominant and approached the input Gaussian in convolution. The same dependence can 
be observed in the intensity of Gaussians. 

For the Gold algorithm, the results were much better. From Fig. 1, it is clear that the FWHM from the fit 
of the deconvoluted peak reached the correct value of 13 for the number of iterations larger than 1000. The 
variance of the FWHM from this fit was < 0.03% and of the amplitude < 0.005% for ti > 1000. 

For the Van Cittert algorithm and N > 2000, the variance of the FWHM from the fit was < 0.3% and of the 
amplitude < 0.05%. 

In Fig. 2 the weighted differences of the original input Gaussian and fitted results of the deconvoluted peak, 
after N = 2000 iteration steps are shown. The difference is evident For die Gold deconvolution the variance 
was still less than 7.5% even for the ±25 displaced channels from the center of the peak. For the large number 
of iterations (Af > 30000) there was already defamation in deconvoluted peak, which could be seen as sharp 
peaks in the dependence of weighted differences. This deformation is a consequence of roundoff errors in the 
computer. 

Fig. 3 shows typical results of the deconvolution procedures. The peak, deconvoluted by the Gold's algorithm, 
is almost identical with the original Gaussian after carrying out 2000 iterations. Van Cittert's procedure showed 
negative values and there were oscillatiobs on both wings of the peak on the level of 2% of amplitude of 
peak. Because Fig. 3 is in a logarithm scale the negative values are not displayed. Gold's deconvolution did not 
produce negative values and oscillations were very small (< 0.01%). 

3.2. The deconvolution of Gaussians with different FWHM 

To determine the influence of the FWHM of the convoluted Gaussian on the result of deconvolution, we 
also made a convolution of the Gaussian with FWHM = 20 with RF. Then the difference in parameters of the 
resulting deconvoluted peak after 2000 iterations was compared with the difference for FWHM = 13. 

The Van Cittert's deconvolution of convolution of the original Gaussian with a FWHM = 20 was in very 
good agreement with the original Gaussian, x 2 fro™ fit 

was 0.04 and the FWHM was exactly the same 
The weighted difference of the fit of deconvoluted and original Gaussian was < 10% for counts > 0.1% 
of amplitude. The situation was worse for the convolution of the Gaussian with FWHM = 13. The FWHM 
of the deconvoluted peak was larger than about 1.8%. But the worst fact was that the x 2 from the fit was 
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very large (Fig. I) and the deconvolution was teller described by two Gaussians (x* = J. 8). That means that 
the deconvolution changed the shape of the peak in such a way that the deconvoluted peak could be better 
described by two Gaussians though there should be only one. The weighted difference of the deconvolution 
tit and the original Gaussian was ^ 10% for counts > 3% of the amplitude. The parameters from the fit of 
deconvoluted peak were in excellent agreement with the original Gaussian. Problem was that for a Gaussian 
with FWHM = 13 there was a better for two Gaussians as for one. 

For the Gold deconvolution agreement with the original Gaussian was very good also for Gaussian with 
FWHM = 13 as we can see from Fig. 2. What was important was that there are no deformations in deconvoluted 
peak because the x 1 for one Gaussian was very small (Fig. 2). 

J.J. The deconvolution of the sum of two Gaussians with different intensity and position 

Tb verify the possibility of distinguishing components of a doublet (sum of two peaks) by deconvolution* 
we made convolutions of the sum of two Gaussians with intensity ratios 1:1, 1:0.8, 1:0 J, 1:0.2 and mutual 
displacements of 1, 5, 10. IS, 20 channels. Both input Gaussians had FWHM = 13, the number of iteration 
steps N was 2000 and the deconvoluted peaks were fitted. We computed the weighted difference of these fitted 
parameters with the parameters of the original input Gaussians. 

For the \feh Cittert's deconvolution we got the following results: 

If the peak displacement was 15 and 20 channels there was no problem to fit the deconvolution results. The 
variation from the parameters of the input Gaussians was small. $ 2%. For a displacement of about 10 channels, 
the fitting procedure for the ratio of intensities 1:0.2 failed. The variance was < 5% for ratios 1:1,1:0.8 and 
< 19% for 1:0 J. For a displacement of about 5 channels we were not able to get a fit for 1:0.2 and 1.-0.5 
ratios. The variances were < 42% (the largest for the amplitude, for the remaining parameters < 10%). So for 
peak displacements larger than the FWHM of the peaks, the resolution was good also for the ratio of intensity 
of 20%. For a displacement smaller as the FWHM of the peak, the resolution was possible only for intensities 
close to each other. For a small displacement (< 10%), the fit was impossible. 

for the Gold deconvolution we also got a (it for a displacement of 10 channels and ratio 1:02 and a 
displacement of 5 channels and ratios 1:0-5. \ :02 with variances < 27%. For a displacement of 1 channel the 
fitting procedure tailed. So we can say that Gold's algorithm is much more able to distinguish a doublet and 
give us more precise results than the Van Chterfs algorithm. 

In Fig. 4 convolutions of Gaussians and the results of the deconvolution procedure described above can be 
seen. 
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4. Experimental spectra 



All these studies of the convolution of Gaussians and following deconvolution were aimed at finding out 
the possibility and limitations of the deconvolution methods and their possible use in the deconvolution of 
experimental y-ray spectra. Because the methods are quite reliable (especially Gold's method) we have used 
Gold's deconvolution method for our experimental spectra. 

The effect of DBAL is quite smallscale and one needs a very stable measuring system for its measurement 
We can get information about its stability from the 1274 keV peak* which is measured simultaneously with the 
511 keV annihilation peak if we use ^Na as a radioactive source. 

In our experiment we used an HPGe detector with a resolution of 1.56 keV at 570 keV with a ^Bi source. 
The energy calibration of the system was 163 eV/channeL The number of counts at the annihilation peak 
was ^630000 and the number of counts at the 1274 keV peak was «180000. In Fig. 5 we can see the 
^-parameter with physical information about the momentum distribution of the electrons and the ^parameter 
of the 1274 keV peak, which carries information about the stability of the measuring system. If there was no 
instability during the measurement it should be a constant As we can see, the ^parameter of the 1274 keV 
peak was not constant so the measuring system was not stable. So we do not know if the structure in the 
3-parameter of 511 keV peak was really a consequence of this instability or if it had a real physical meaning. 
Of course, there is a possibility to try to use the ^parameter of the 1274 keV peak as a correction on the 
instability, eg. by calculating the ratio R = 5(511 keV)/5(1274 keV). The result of such procedure is shown 
in Fig. 5. As we can see, this simple procedure does not suppress the effect of instability. Because we are 
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interested only in the relative changes of the S-parameter on temperature, we can make a deconvolution of the 
511 keV peak with the 1274 keV peak and try to suppress the effect of instability. 

Fig. 6 shows the result of these deconvolutions (the 5-parameter after deconvolution) along with the S- 
parameter from deconvolutions of the 511 keV peak with the 570 keV peak (^Bi), which was measured 
at the beginning of our experiment So we can see that the structure at die temperature of 175 K was only 
a consequence of instability of the measuring system, because it disappeared after the deconvolution of the 
511 keV peaks with the 1274 keV peaks. We can obtain physical information from the dependence of the 
.^parameter on the temperature by using a deconvolution procedure. 

This procedure can also be used for a reliable estimate of the electron momentum distribution. In such case 
we need a better resolution function (e.g. 514 keV from the decay of "Sr) measured simultaneously with the 
annihilation tine. 



5. Conclusion 



It turned out that the described deconvolution methods are in principle convenient for processing experimental 
spectra from DBAL measurements. They help to sup press influence of instability of the measuring system 
Particularly the Gold deconvolution method works very well and it does not deform the deconvolution result 
The reproducibility of the deconvolution results tested on simulated convolutions of Gaussians is very good. 
Vhn Cittert's algorithm gives the results with oscillations from positive to negative values on the wings of the 
deconvolved peaks. Gold's deconvolution gives only small non-negative oscillations. Unambiguously we can 
say that Gold's algorithm is more reliable with regard to the precision of the deconvolution. The time needed 
for the deconvolution is practically the same for both algorithms. 

So, the use of Gold deconvolution is a cheap and simple way to eliminate the electronic instabilities which 
blur physical effects. 
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